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GIVES AN 

ABSOLUTE DEPTH 
REFERENCE TO THE 
ELECTRICAL LOG FOR 
PERFORATING CONTROL 
IN WORKOVER 




















Many workover wells have the original casing 
set just above a sand in which the well is to be recompleted 
Recently, it has become standard practice to run a casing collar 
in combination with the regular electrical log. By so doing 
establish definite distances between the casing collars and the project 
When the well is perforated, the same collars can be relocated 
and the shots placed with absolute control 

You should consider running a Casing Collar-Electrical Log 


combination on your next workover well 





The Schlumberger Casing Collar Locator is another impor- 
tant development of Schlumberger’s progressive program 
of research and engineering to provide advanced services 


to the oil in 
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CHAMPIONS DON’T WIN BY PULLING PUNCHES 


To win consistently in the ring a 
man has to punch as hard and straight 
as he knows how. And to win consist- 
ently in the sometimes uncertain busi- 
ness of producing oil, you have to 
use casing perforating that punches 
straight through all casing, cement 
and far back into the formation. The 
record of results in thousands of 


wells perforated with the 


modern Welex Jet Process proves 

that JETS give 300% deeper penetra- 

tion than any other practical method. On 
job after job results have proven that 
JET perforating can make a profit- 
able oil well where other methods 
failed. Men who know will tell you 

that results cost less with Welex JETS. 

Call your nearest Welex Station for 


prompt service day or night. 


Welex 
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drillingand 
production techniques 


For a sure-fire’ shut-off 
use this Baker Casing 
Bridge Plug (Product No. 
400-D) which can be run 
and set either on tubing 
ora wireline. * For 
squeeze cementing, aci- 
dizing or testing, use the 
old reliable Baker Ce- 
ment Retainer, Product 
No. 400. 





The Baker Retainer Pro- 
duction Packer (Product 
No. 415-D) meets most 
production requirements 
with unequalled flexibil- 
ity, safety and economy. 
* Can be set on tubing 
or drill pipe—(or on a 
wire line by your choice 
of the leading electric 
wire line service organi- 
zations). 





Don’t waste time and 
money hauling kinked 
Kellys and bent drill pipe 
to the shop! * The Baker 
Portable Kelly and Pipe 
Straightener (Product No. 
550) does a better job while 
the Kelly or pipe is sus- 
pended vertically through 
the table; and avoids pos- 
sible damage in transport- 
ing and handling. 


Use Baker BAKWIK Drill 
Pipe Floats (Product No. 
480) to reduce wear and 
tear on wire lines and rig 
equipment. * Greatly re- 
duces possibility of wet 
strings due to clogged bits 
and drill pipe. 


When your casing 
is equipped with 
Baker Model ‘’G” 
Casing Centralizers 
(Product No. 
910-G) you know 
that is REALLY 
CENTERED with 
better chances for 
a ‘first-time’ 
water shut-off. * 
* In directional 
holes “sled-runner- 
action” of the 
springs helps to 
get pipe to bottom. 


No more torn 
swab rubbers or 
packing elements. 
Slips on testers or 
packers can always 
contact casing. 

* A Baker ROTO- 
VERT Casing 
Scraper (Product 
No. 620-C) run on 
drill pipe, or a wire 
line, scrapes casing 
free from gun- 
shot burrs and 
hardened cement. 


HOUSTON + LOS ANGELES 


For successful 
first-time” 
cementing use a 
Baker Cement 
WASH-DOWN 
WHIRLER Float 
Shoe (Product 
No. 120) cen- 
tered in the hole 
by a Baker Cas- 
ing Centralizer. 
* If an addi- 
tional ‘float’ 
valve is desired, 
run a Baker 
Cement Float 
Collar (Product 
No. 101) up the 
string. * If one 
“floating unit’ 
is sufficient for 
your casing 
string, use a 
Baker Float 
Collar a joint of 
two up the 
string anda 
Baker Cement 
WASH-DOWN 
WHIRLER Guide 
Shoe (Product 
No. 121) on 
bottom. 


BAKER OIL TOOLS, INC. 


NEW YORK 








DOWELL STRATALIFT using JEL-ACID provides an improved 
method of increasing production by fracturing formations 
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Editorial  Saeree Ya 
On Stockpiling Oil 


By R. W. French 


Sohio Petroleum Co 


I: YOU review our past editorial pieces, you may be struck as I was by their 
variety — despite our stratified professional training and immediate interests. 
From problems-and-progress of our AIME organization itself, they touch on many 
phases of development and use of engineer-citizens, winding up in such topics as 
greater participation in politics and estimation of undiscovered reserves. 
Still further reflection suggests a curious relationship to one vital common prob- 
lem now facing us — subsurface stockpiling of oil! Why a problem? In what con- 


nection? Here is one view, in brief: 


WHY A PROBLEM? 


Although our domestic supply and reserves (reserves are significant only as 
tomorrow's supply) may continue to increase as they have historically, we 
still face increasing imports or substitutes to balance climbing demand, even 
in peacetime 

Despite curtailment of “civilian” demand, and uncertainty of actual oil require- 

ments for fuel, lubrication, and chemical material, we have no proof that ade- 

quate domestic supply will be available in time of war. 

Unless such unlikely proof is forthcoming, we face the choice of stockpiling oil 

or substituting other energy sources such as coal, shale, gas, alcohol or even 

hydro-electric, solar, or nuclear energy, etc., through displacement or conversion. 
IN WHAT CONNECTION? 

If we agree that there is some problem, how do we think it should be met (1) as 
engineers or business men, (2) as citizens? Immediately we can ask knotty ques- 
tions, theoretical, technical, and “practical.” such as: 

(a) How much oil should be imported to defer current domestic production for 
stockpiling ? 

(b) Would it ever be feasible or necessary to actually inject oil into natural under- 
ground reservoirs for storage. with the obvious cost disadvantages and uncer- 
tainties of later recovery ? 

Who would own stockpile oil and what compensation or arrangements would 
be made with owners, operators, and royalty holders? 

Who would determine when the oil should be produced? At what price? For 
whose refineries? 

What attitude and action would be taken by state and federal tax and regula- 
tory authorities? And would this be impossible without federal control? 

But how does this all relate to our past editorial topics? Take “registration of 
engineers” by Geis — this says engineers must have certain competence, certainly 
have responsibility to the public. “Engineering education” by Calhoun, and Gay- 
lord’s excellent statement on “specialization” covered the academic and training 
phases basic to successful discharge of engineering work today. 

Then came views by Bullard and Rubel urging engineers to extend their efforts 
beyond engineering, to interpret and act in community and citizenship matters as 
aggressively and effectively as they do in their professional work. Returning to 
DeGolyer’s provocative and timely contribution on estimation of undiscovered 


a we see that it too bears directly on this topic as one important factor in 


xckpiling” to assure future supply. 

To me, these random topics bear on many aspects of today’s problems but with 
particular force on our nation’s future supply of crude oil. Those who heard Rail- 
road Commissioner William Murray's timely and forthright address at the AIME 
Annual Meeting in St. Louis last month will not forget his plea for realistic engi 
neering work as a requisite of national defense. 

Whether your primary work is in exploration, drilling, production, research or 
management, this “stockpiling” problem and its solution affect you professionally. 
Regardless of occupation, it concerns every American ultimately. Individually we 


must give it our most objective and far-sighted attention x * * 
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THE DEVELOPMENT 


and OPERATION of 


LPG STORAGE CAVITIES IN SALT STRATA 


By K. C. Howard 


Richardson and Bass, Fort Worth 


The process of storing liquefied petroleum gases in reservoir salt formations has proved 
successful in West Texas. Though untried, the process is believed to be applicable to forma- 


tions other than salt. 


ABSTRACT 


Underground storage of liquefied petroleum gases in reset 
voirs dissolved by circulating fresh water in Upper Permian 
salt beds in West Texas has been very successful. One system 
of 55,000 bbl capacity is in operation in Upton County, Tex.. 

another system in Winkler County, 
approximately 25,000 bbl capacity and is in the process of 


while Tex., is now at 
being enlarged. 

Development of the Winkler County project was delayed at 
about 8,000 bbl capacity for the 
allow demon trations of the system to the industry. In these 


purpo-e of testing and to 
tests, recovery of injected propane varied from 95.9 to 99.7 
per cent of injected volumes. The recovered LPG in all tests 
had increased in vapor pressure and had picked up a slight 
trace of moisture. 

The process, though untried, is believed to be applicable to 
formations other than salt such as shales using partially sat- 
urated salt water, and possibly to limestone formations using 
hydrochloric acid. 

Fall Meeting in 


was presented at the Petroleum Branch 


New O La., October 4-6, 1950 
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DESCRIPTION OF PROCESS 


It is believed that the practicality of artificially formed 
underground LPG storage space is established by the per- 
formance to date of two such systems in West Texas. In these 
systems, storage space was dissolved out of salt beds by cir- 
culating fresh water. 

The Upper Permian salt section of West Texas is generally 
a little over 1.000 ft thick, and in the well in the 
Keystone Field, the top of the salt is at 1.170 ft. The salt is 
overlain with a stratum of good solid anhydrite 100 ft in thick- 
ness in which 95,-in. casing was set at 1,162 ft. The hole was 
bottomed at 1,770 ft in salt and 2.5 tubing was run to 1,748 ft. 


storage 


ON EFFLUENT WATER 


RAT 


SAT 











FIG. 2 — SALT SATURATION OF EFFLUENT WATER 


\ conventional shallow well Christmas tree was used and 
hooked up with a positive displacement pump so that fresh 
could be injected either into the tubing or into the 


annulus 


water 


Fresh water was injected at 5.5 bbl per minute. The first 
returns were only 10 per cent saturated with salt, but as the 
hole was enlarged thereby increasing the area and time the 
salt was exposed to water, the salt saturat‘on of the effluent 
water continuously increased and leveled off at approximately 


97 per cent saturation after about 14 days circulation. The 


efluent water saturations for the first few days’ circulation 
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FIG. 3 — CALIPER LOG AFTER CIRCULATING 50 HOURS 


are shown in Fig. 2. At present the capacity of the Winkler 
County system is about 25,000 bbl and circulation will be 
continued until a capacity of 50,000 bbl is obtained. 

The salt beds consist mostly of sodium chloride with some 
chlorides of potassium and magnesium interspersed with an 
hydrite stringers and other insoluble matter. It was estimated 


that at 100 per cent saturation, 6.1 bbl of water will dissolve 
one barrel of formation and the computed volumes of space 


formed using this estimate checked very closely with the dis 
solved volume indicated by a caliper survey. The solubility of 
salt in water varies only to a small degree with temperature 
and the estimated solubility of the salt in the Winkler County 
project at an average circulating temperature of 80°F is shown 
in Fig. 1. It is believed that the volume of salt dissolved out 
of a formation may be determined with sufficient accuracy 
for all practical purposes from the quantity of fresh water 
injected and the salinity of the effluent water. The relative 
solubility of various strata in the salt formation is shown in 
Fig. 3, which is a representative portion of a caliper survey 
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run after approximately 50 hours circulation. The dissolved 


volume at the time of this survey was 1,115 bbl. 

No detectable loss of fluids was caused as surface pressures 
in excess of 500 psi were applied to the annulus on various 
occasions while circulating. This surface pressure plus the 
hydrostatic head of the circulating fluid caused a pressure 
against the salt formation equivalent to a gradient of 0.86 
psi per foot. Consequently, as no such pressures are required 
in normal operation, the possibility of loss of stored products 
to the formation is considered to be remote. 

As the Upper Permian salt section has shown no tendency 
to break down at excessive pressures, it has been found most 
convenient to store the LPG above water in the sterage space 
by injecting LPG into the annular space displacing the water 
down and up through the tubing te the surface. LPG is recov- 
ered by injecting water into the tubing thereby di«placing 
the LPG back through the annulus 

The surface pressure of the LPG if stored over water is 
proportional to the depth of the LPG-water contact as the 
system is merely a “U” tube with salt water inside the tubing 
and LPG in the annular space. The surface pressure of the 
LPG is equal to the difference in the columnar heads of the 
two fluids. The difference in columnar heads of saturated salt 
water and propane is about 300 psi per 1.000 ft. Storage of 
LPG obviously necessitates sufficient pressure to prevent va- 
porization of the product plus sufficient pressure to conven 
ently transfer the product, or at least furnish adequate pres- 
sure to the suction end of a transfer pump. A formation tha 
will support a column of saturated salt water and as much 
as 700 ft deep. may be used satisfactorily for the storage of 
200 psi vapor pressure propane, and of course, more shallow 
formations may be used for the storage of products with le-s 
vapor pressure 

The process is deemed to be practical in formations that 
will not support a column of water, provided the formation 
is at sufficient depth. Several schemes may be employed. For 


Continued on Page 3, Section 2 
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A Report by the 


SouTH 


AIME PRODUCTION REVIEW COMMITTEE 


PRODUCTION 


‘asim | 


D. V. Carter, Chairpran = \- 
R. B. Gilmore, Vige- Chairman 
“SEW. Berlin, Vigg-Ehairman 


\.. Milton Williams “ 
: 


Lawrence, W. Chasteen 
Gerald-€. Eddy = 
| Gordon Fisher qi. 

| “Raymond M. Lars¢én 
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John K. Werightur: 


s\ear Split Between Warand Peace Brought An Upsurge In Pettoleum Productian@ar-d Caused Fluc- 
tudtiens In Exploration Activities Throughout the Domestic.Area. Heré«ls a —— of 1950 Production 
Coordinated | by Ten Outstanding Industrymen. 


; abrupt.change in the world political situation in mid- 
1950..a wirtval Overnight ehange from {one of ignoring 
trends toward\ war to actual conditions of war, 
alterations in) praduction practices in the domestic petroleum 
industry. Productiom for the year was up all over the nation. 
Curtailment» were lifted and allowables ‘adjusted upwards. 
Such was not\ehe éase. however, in development and explord- 
tion drilling. Such ‘activity fluctuated across the coumMry with 
increases recorded in} sonfe sections and decreases in othegs. 
In areas of stro¥g petroledm reserves, there marked 
increases in activitié Of all fepes in exploration, drilling and 
production. 


brought many 


were 


There follows a summary of activity in the various sections 
of the United States. The data presented here were gathered by 
members of the AIME Production Review Committee and were 
reported upon in greater detail at the AIME Annual Meeting 
in St. Louis. Mo., Feb. 22 


ARKANSAS 
Production 


Cumulative oil and condensate production in Arkansas at 
the end of 1950 totalled 815.758.477 bbl. During the year the 
output of oil plus condensate totalled 30,840,007 bbl or 1,054, 
594 bbl more than the annual production of 1949. Runs to 
stills during the year amounted to 20,448,082 bbl or an in- 
crease of 2,721,829 bbl over 1949, and the intake volume of 
the gasoline plants increased from 59,752,686 cu ft in 1949 
to 65,292,690 cu ft in 1950. Dry gas production in Northwest 
Arkansas increased from 6,900,678 cu ft in 1949 to 7.802.827 
cu ft in 1950. 


Exploration and Drilling 


Geophysical exploration in Arkansas was confined almost 
entirely to the southern half of the state and decreased from 
610 crew weeks in 1949 to 467 crew weeks in 1950. About 
98 per cent of the drilling activity in Arkansas during 1950 
was confined to the southern half of the state and only 2 per 
cent to the northern half. A total of 425 wells were drilled 
for oil and gas during the year. representing an increase of 
66 wells over the number drilled in 1949. The 425 wells con- 
sisted of 103 wildeats and 322 field wells. 
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Twelve of the wild- 
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cats opened up new oil or gas cohkdéhsate producing areas 
in Séaitlk Arkansas and of these, teh discovered new fields 
and two were considered field extensions due chiefly to their 
proximity te eld fields Anestieér wildcat well opened a dry 
gas field in Nerthwest Arkansas, and the remaining 90 wild- 
cyts were abaftidened as dry. 

Of the 322 field wells drilled dufing thd year (which do not 
include the 13 discovery wells), 24%-were producers and 75 
were abandoned as dry. Also 82 old, producing wells, were 
recompleted as producers in the same of“other nes artd B8 
former producers were abandoned. 


Secondary Recovery Projects 

Excellent results were obtained by secondary methods of 
recovery in the oil and gas fields of South Arkansas. These 
methods proved successful in the following reservoirs: Haynes- 
ville, Lisbon, Magnolia, McKamie-Patton. Midway, Schuler. 
Stephens and Wesson. 


CALIFORNIA 


Crude oil operations for 1950 could best be described as 
reversing the trends of 1949 from a declining production rate 
to an increased production rate. At the beginning of 1949, 
the outlook was one of overproduction with the threat of 
while by the end of 1950 a period of expansion, 
production, was anticipated. The prin- 
was the threat of 


curtailment, 
both in drilling and 
cipal concern among crude oil producers 
material and manpower shortages which might result from 


defense 


Crude Oil Production 


Production at the beginning of 1950 was 866,969 B/D, 
which was approximately the rate in late 1945 and early 1946. 
After 1946, production steadily increased to a peak of 954.,- 
000 B/D realized in November, 1948. This rate was the all- 
time peak for the state until it was later broken by a produc- 
tion of 963,400 B/D in December, 1950. This new peak was 
especially significant because Elk Hills, which produced at the 
rate of 65,000 B/D in World War II, was practically shut in. 
In addition, there was also an estimated shut-in capacity of 
18,000 B/D in other fields. 

For the first four months of 1950, the state continued the 
started in early 1949 and reached a minimum of 


preparations. 


decline 
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857,439 B/D in April. In May. production began a rapid 
rise to the peak reached in December of 963.418 B/D. 


Development 

During 1950, some 1,262 oil wells and 46 gas wells were 
reported completed. Their total initial production rate was 
262,800 B/D. The new well completion rate reached a peak 
at the end of the year with 122 oil wells completed in Decem 
ber. The 1950 completion figure represents a decrease of 461 
under the 1949 development figure. However, initial produ 
tion increased some 25.800 B/D. Gas well completions in 
creased by 11 over the 1949 number. 


Exploration and Discoveries 

During 1950, a total of 489 exploratory wells were drilled 
as compared to 588 in 1949, representing a 17 per cent 
decrease. New field wildcats totaled 277 in 1950 as compared 
to 340 in 1949. A total of 11 new oil fields and two gas fields 
are credited as discoveries in California for 1950. In addition 
to the new fields, seven new pools and two important exten- 
sions were made. Although the list of discoveries is long 
there were no major additions made to reserves. Based on 
the data available, it is doubtful if California’s reserve posi 
tion will hold its own for 1950. 

Exploratory Activity in Other Pacific Coast Areas 

Exploration work on the Pacific Coast in areas outside 
California was highlighted by the shifting of activity from 
the Pacific Northwest to the Great Basin area of Nevada. 
Some 15 companies were actively engaged in geologic work 
in Nevada and by year’s end more than one million acres of 
Government land were under lease. One well was abandoned 
at 10,314 ft and two other wells were drilled without success 
in 1950. At year’s end, two more wells were being drilled, one 
in Clark County and the other in White Pine County. 

In Arizona, nine wells were drilled and were either aban 
doned or idle by December 31. Two wells were drilled in 
Oregon but neither were successful although reports of gas 
showings were made from a well drilled in Lake County 
In Washington, drilling activities ceased with the abandon- 
ment, after test, of Union Oil Co.’s State No. 3 in Grays 
Harbor County. One well spudded in Clallum County on the 
Olympic Peninsula in 1950 but was idle at the end of the 
vear. 

Petroleum Supply and Demand 

The Pacific Coast total petroleum supply and demand situ 
ation underwent a sharp change in 1950 compared with 1949. 
On the basis of data published by the Bureau of Mines for 
the first 11 months, total petroleum supply declined about 
2 per cent to around 990,000 B/D from 1,010,000 B/D in 
1949. Total demand, on the other hand, increased to 1,070,- 
000 B/D, a gain of about 11.5 per cent over the year earlier 
level of 960,000 B/D. As a result of these diverse trends, total 
stocks of all oils on the Pacific Coast declined about 28,500,000 
bbl for the year which was in marked contrast to the increase 
of 18,400,000 bbl in 1949, 

Contributing to the over-all improvement were gains of 9.7 
per cent in total demand for gasoline, 6.25 per cent in stove 
and diesel oil, and 15.1 per cent in fuel oil. An important 
factor in the gain in fuel oil demand was the shipment of more 
than 40,000 B/D to the East Coast. This eastward movement 
of fuel oil ended in the fourth quarter of 1950. 

The year proved to be one of sharp improvement over the 
unhealthy statistical position which had been built up previ 
ously. Through the curtailment of crude oil production and 
refinery activity, inventories were reduced to levels generally 
considered to be near the economic minimum in relation to 
the rate of demand. A reversal of both supply and demand 
trends occurred during the year, becoming evident in the 
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second quarter. In large measure, the improvement in the 

position of the Pacific Coast petroleum industry reflects the 

stimulus of the Korean War and the mobilization program. 

In this connection, it is interesting to note that the Pacific 

Coast petroleum industry is supplying more than one-half 

the total United States military petroleum requirements. 
ILLINOIS 

In Illinois 2.894 wells were drilled for oi] and gas in 1950 
as compared with 2,74! in 1949, an increase of 5.5 per cent. 
rhis number includes 830 wildcat wells which discovered 25 
new pools and 77 extensions to pools. There were no major 
pools among the 25 discovered in 1950. 

An outstanding event in the oil industry of Illinois was the 
beginning of water flooding in the Salem Pool, Marion County, 
under a voluntary unit operating agreement. Up to the end 
of 1950, the Salem Pool had produced approximately 216 
million bbl of oil in 12 years and it is expected to produce 
another 200 million bbl by flooding. Other Illinois water flood- 
ing operations are progressing favorably. Total oil production 
in Illinois in 1950 was 62,123,000 bbl (preliminary figure) 
valued at approximately $172,081,000. This represents a de- 
cline in production of approximately 4 per cent from the 
64,583,000 bbl produced in 1949. 

INDIANA 

4 total of 1,527 wells drilled for oil and gas were completed 
in Indiana in 1950 as compared with 1.288 in 1949, an increase 
of 18 per cent. Total oil production in 1950 was some 10,- 
702,683 bbl as compared with 9,734,000 bbl in 1949, an in- 
crease of 10 per cent. Wildcat drilling amounted to about 
one-third of the total drilling. Of these 504 wells, 28 discov- 
18 new plays; 26 extensions; and 432 were 
Most of the wells drilled were in Southwestern 


ered new pools; 

dry holes. 

Indiana. 
KANSAS 

Kansas operators posted new record highs in both the 
number of wildcat wells and total well completions during 
1950. Wildcat drilling increased by a phenomenal 56 per cent 
over 1949 to a total of 816 wells. Total well completions were 
up 20 per cent to an all-time high of 4,077. Oil production 
increased by 7 per cent to slightly over 107 million bbl, the 
second annual production recorded. Total footage 
drilled was nearly 12.75 million ft, an increase of 19 per cent 
over the previous record 1949 figure. 

The development of proven fields increased about 14 per 
cent over 1949. Drilling continued at a rapid rate in the 
Kraft-Prusa, Morel, Trapp, Silica, El Dorado, and Berland 
fields. Some 373 new gas wells were reported drilled in the 
Hugoton Gas Field. Permits were granted for 61 new water- ' 
flooding projects, chiefly in the East Kansas stripper area. 
It is estimated that more than 10 per cent of the current oil 
production in Kansas is from secondary recovery operations. 
In spite of the sharp increase in the number of exploratory 
wells drilled, the success ratio declined only slightly, from 
19.8 per cent to 17.4 per cent for 1950. One of the outstanding 
discoveries of the year was the Fall Creek Field in Sumner 
County. The discovery well was complete in August with a 
calculated potential of 6,226 B,D from the Simpson dolomite. 
Four wells had been completed in this field at year’s end. 
Exploratory activity in the Central Platform area accounted 
for the majority of the discoveries during 1950. 

KENTUCKY 

Kentucky produced 9.205.707 bbl of oil during the first 11 
months of 1950, and with December yet to be calculated at 
this writing, it is expected that the figure will exceed the 
10,000,000 bbl mark for the second time in the state’s history. 
In addition, 70 billion cu ft of gas was produced during the 
vear, of which 95 per cent came from the Big Sandy Gas Field 


largest 


Sowa ees oe 








of Eastern Kentucky. During the 1,226 wells 
drilled of which 524 were oil wells, 193 gas wells and 509 
irv holes. 

The most spectacular development in Western Kentucky 
during 1950 was in the Zion North Pool of Henderson County 
where 48 oil wells out of 59 wells drilled produced oil from 
the MeClosky limestone. The total initial production of the 
18 oil wells amounted to more than 12,000 B/D. The largest 
operations were in the Morganfield South in Union County 
where 66 wells were drilled, of which 57 were oil wells. In 
the Guffie Pool in McLean County, 48 oil wells and 12 dry 
holes were drilled. As of November, daily pipe line runs were 
1,700 B.D. The Griffith Pool in Davies County had the fourth 
largest number of operations with 35 oil wells, 


year, were 


three gas wells 
and 20 dry holes with pipe line runs of 705 B/D as of the 
same date. The only new oil development in Eastern Ken- 
tucky was small extensions of the Weir Oil Fields of Elliott 
County on the west flank of Paint Creek Uplift and extensions 
of the Corniferous White Creek Pool in Wolfe and Powell 
counties. 
LOUISIANA 

1950 witnessed an increase of 43 active wells over 1949 in 
Louisiana. The active wells include all wells drilled and drill- 
ing prior to January 1, 1951. There were 2,631 field wells and 
331 exploratory wells drilled resulting in a total of 2,962 
active wells for the year. At the end of the year. 174 field 
wells and 7] wildcat wells were drilling. 

In comparing the number of wells drilled in North Louisi- 
ana and South Louisiana it is found that there were more 
field and wildcat wells drilled in the northern part of the 
state. This resulted in a number of successful field 
wells in North Louisiana but not of successful wildceats. In 
Northern Louisiana’s districts there were 1.367 wells 
drilled and of this number 1,220 were field wells, 709 oil, 256 
gas, 255 dry holes and 147 wildcat wells. In South Louisiana 
1,083 wells drilled. Wells drilled in 
areas resulted in the discovery of 30 new fields, 23 vil and 
seven gas. Four of the 30 new discoveries may develop into 


greater 


two 


there were unproven 


major fields. They are Cox Bay, Romere Pass. Glenmore and 
Unionville. 
Production 
Production of oil 1950 in- 
Oil production (in- 
15.786.051 bbl. 
amounting to a 7.98 per cent increase. Gas production, both 


and gas in Louisiana during 
creased in comparison to 1949 figures. 
cluding condensate) showed a net gain of 
natural and that produced in conjunction with oil, increased 
by 93,946,807 cu ft. a 9.16 per cent 1949 gas 
production, 


Offshore Activity 


The number of successful wells completed in the Gulf of 


gain over 


Mexico during the year totaled 47, as compared with a total of 
4 in 1949. During 1950, 98 wells were drilled, 37 completed 
as oil wells, 10 as gas wells, 34 abandoned and 17 drilling at 
the end of the Marchand area showed the 
increase in 1949. A total of 25 
were completed, 18 as oil, one as gas and six abandoned 
MICHIGAN 

Estimated oil production for Michigan for 1950 was 15.- 
800,000 bbl. This figure compares disfavorably with the 1949 
total of 16.503.473 bbl. Estimated dry natural gas production 
billion 


from 


vear. The Bay 


greatest activity over wells 


is 10 billion cu ft, a decrease of approximately 4.5 
cu ft from 1949 production. this 
domestic production was more than offset by 


However. decrease 
increases in 
imported gas both to the markets and into underground stor- 
age. It is estimated that during the vear 90 billion cu ft were 


imported. 
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During the year, 901 oil and gas permits were issued. Of 
the 890 completions, it is estimated that 325 were oil wells, 
gas wells 190 were dry holes. Total 
drilled for the estimated at 2,150,000 ft of which 
probably 800,000 ft was wildcat drilling. 


75 were and footage 


year 1s 


Sixteen new fields, six new pays and seven extensions were 
discovered during the year. Of the new fields, Cannon Creek 
was probably the most important. Seventeen wells have been 
drilled, 


defined. 


proving up 2,720 acres, and the field is not as yet 
Ihe Howell-Salina gas field with 16 wells and an 
area of 2,600 acres during the year produced 3.75 billion 
cu ft 
MISSISSIPPI 

In 1950, there was a further decrease in the development 
drilling in Mississippi with only 177 development wells 
drilled as compared with 218 drilled in 1949. Of these 177 
wells, 127 were successfully completed as producers. However, 
on the upswing during 1950 with 
drilled compared to 117 in 1949. 
These wildcat wells discovered eight new producible fields. 
The year 305 wells drilled in Mississippi. During 
the period, the state produced approximately 38,500,000 bbl 
of oil and condensate from 44 oil and gas fields. This figure 
brings the cumulative total for the state since August, 1939, 
some 285,- 


exploratory drilling was 


128 wildcat wells being 


1950 saw 


when oil was first discovered in Mississippi, to 
300,000 bbl of oil and condensate. 

MISSOURI 
wells were completed in Missouri during 
of five over the 1949 figure. Of these, two 


A total of 27 
1950, an 
were completed as oil wells and 10 as gas wells. Seven ex- 
ploratory wells were reported drilled during 1950, all of which 
were dry. Oil production increased slightly to about 36,000 
bbl for 


increase 


the year 
NEBRASKA 
Drilling Nebraska has increased rapidly since 
the discovery of the Field in Cheyenne County in 
1949. Total completions were up from 27 to 171 wells while 


activity in 
Gurley 


wildcat wells increased from 18 to 82, of which 13 were suc- 
cessfully completed as oil or gas producers. Crude production 
for 1950 was approximately 1.25 million bbl as compared with 
330.000 bbl for 1949. 

Exploratory drilling resulted in the discovery of eight new 
oil and gas fields in Western Nebraska during 1950. Six of 
10 miles of Sidney, 

County. One new gas field is located near Big 
Duell County, about 
new oil field is located near 


Kansas border. 


these fields are within a radius of about 
Cheyenne 
Springs. 15 miles east of Sidney, and a 


Alma, Harlan County, near the 


NEW MEXICO 
New Mexico produced approximately 47,300,000 bbl of oil 
during 1950, an average of slightly less than 130,000 B/D 
There were approximately 6.100 wells at the end of the year. 
There were 607 wells drilled in the state during the year, of 
which Seventy- 
nine of the wells drilled were wildcats and 23 of these were 


84 were completed as oil or gas producers. 


completed as oil or gas wells. Cumulative production from 
New Mexico at the end of 1950 was approximately 681,000,000 
bbl with estimated reserves of approximately 610,000,000 bbl. 
There were 16 discoveries in the state during 1950, including 
new pays in previously producing areas. 
NEW YORK 
Oil production in New York for 1950 is estimated to be 
1,130,000 bbl. representing an approximate decrease of 6.1] 
per cent from 1949 production figures. During the year, 700 
field wells were drilled in all fields with 567 completions in 
the Allegany Field. This figure represents a decrease of 20 
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per cent under the drilling program of 1949. This decrease 
in drilling activity is probably due to the price of oil and to 
drilling only in proven areas with very little wildcat drilling 
being carried on. The total number of producing wells in 
New York, all stripper type, is in the neighborhood of 23,500 

The Medina and Oriskany gas production in New York for 
1950 is estimated at 4 billion cu ft. 


OHIO 


Total number of wells in Ohio during 1950 was 1,141. of 
which 284 were gas wells, 282 oil wells, 100 combination wells, 
and 475 dry holes. Total gas production amounted to 143,833 
cu ft and oil production totaled 10,328 bbl for the year. No 
new gas pools of size were discovered during the year and 
drilling in the gas areas dropped off. The principal oil de 
velopment during the year took place in New Castle town 
ship, Coshocton County, in the Clinton Sand. 


OKLAHOMA 

Activity in Oklahoma continued at an accelerated rate dur 
ing 1950. Well completions were up 25 per cent over 1949 
to 5,459, the largest number since 1923. The number of wild 
cat wells increased only about 5 per cent over 1949 to a total 
of 840, about 22 per cent of which were successfully com 
pleted as oil or gas wells. Total footage of drilled wells was 
up nearly 30 per cent to what probably is an all-time high of 
nearly 20 million ft. Crude production during the year in 
creased to a new 12-year peak of more than 163 million bbl. 
largely as a result of increased allowables and new develop 
ment. 

Most of the sharp increase in drilling during the year 
resulted from the accelerated development of such proven 
fields as Elk City and Ringwood in Western Oklahoma and 
the Golden Trend, Palacine and Velma-Sholem Alechem areas 
in South Central Oklahoma. 

Exploratory activity in the western part of the state was 
disappointing. Several deep tests were drilled on presumed 
seismograph highs in the western part of the Anadarko basin 
but all were dry. Probably the most important discoveries of 
the year were located in the south and south central parts 
of the state. Of extreme importance is the discovery and 
extension of deep Springer (Pennsylvania) sand production 
in the Ara, North Milroy, Fox-Graham, and Knox areas of 
Stephens, Carter, and Grady counties. Several narrow but rich 
bands of Springer production have been encountered. gen 
erally off the flanks of proven fields producing from shallower 
zones. One such band appears to lie southeast of Velma and 
along the north flank of the Fox-Milroy high. Two recent 
and probably very significant discoveries were made in Love 
County. in the Marietta basin. Production was encountered in 
one of the Deese sands (Pennsylvanian) in both wells, located 


about two miles apart, which suggests the possibility that 
one large field may have been discovered. Exploratory drilling 
continued at a high rate in the central and north-central 
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parts of the state. resulting in the discovery of a large num- 
ber of relatively small fields, chiefly in Pennsylvanian sands. 
The development of the Hydrafrae process, which has enabled 
made in many of the tighter 
were passed up. has stimulated the 


commercial completion to be 


sands which formerly 


development of this area. 
PENNSYLVANIA 
There was a slight increase in oil production in Pennsyl- 
vania during 1950 due to improved market conditions, but a 
decline of more than 10 per cent in drilling activities con- 
secondary operations, 
which acount for the bulk of the oil production in Pennsyl- 


ducted in connection with recovery 
vania. No new oil pool was discovered. In the shallowed sand 
gas territory, less drilling was carried out in 1950 than in 
1949. There were some extensions of oil pools and several 
new discoveries but none of these gave promise of developing 
into pools of significant size. Thirty-four new wells were com- 
pleted to the deeper horizons (into and below the Onondaga) 
in 1950, as compared with 27 in 1949. Of these, two completed 
for gas storage purposes. Of the remainder, 21 were gas wells 
and 11 were dry holes. 

The oustanding event of 1950 was the discovery of another 
major Oriskany sand gas pool in North Central Pennsylvania 
early in the year. By the end of 1950, twelve producing wells 
and three dry holes had been completed and 42 wells were 
drilling. The pool was producing at the rate of about 150 
MMef per day. 


ROCKY MOUNTAIN AREA 


Oil production from Colorado, Montana, Utah and Wyom- 
ing totaled more than 94,000,000 bbl in 1950, an overall high. 
Wyoming was responsible for most of the increase. Produc- 
tion in January was at the rate of 240,000 B/D and in De- 
cember it amounted to 298,000 B/D. 

Five wildcat wells in Idaho were dry and one deep test was 
still drilling. Nevada had four dry one at Meridian 
being 10,314 ft deep. Upon completion of this test late in the 
year, the drilling equipment was moved to Hayden Creek, 
White Pir County. Several unsuccessful wells 
drilled in ..orth and South Dakota. 

Colorado 


Oil production in 


holes, 


also were 


1950 was 23,140,000 bbl or 440,000 bbl 
less than in 1949. About 42 oil wells and 12 gas wells were 
drilled as compared with 55 plus four in 1949. 
Montana 

Oil production in 1950 in Montana totalled 8,279,000 bbl 
which was 845.000 bbl less than in 1949 and 1,000,000 bbl less 
than the 1949 maximum. Production varied from 20,000 B/D 
in January to 24,000 B/D in December. The only apparent 
new West Sumatra and Northwest 
Sumatra in Rosebud County. Notable extensions or new pools 


field discoveries were at 
were found east and southwest of the Pondera Field. 
Utah 

Oil production in Utah during 1950 totaled 1,208,000 bbl 
or 572.000 bbl more than in 1949. Three oil wells were com- 
pleted at Ashley Creek in 1950 and with 30 wells in the 
Weber sand, development is about complete unless other 
formations warrant separate drilling. A new oil discovery 
was made at Gusher but initial production at a depth of 
9.757 ft in the Green River formation was only 20-35 B/D. 
Wyoming 

Oil production in 1950 in Wyoming established a record 
high of 61,700,000 bbl which was 13.385.000 bbl more than 
in 1949 and 7.000.000 bbl more than in 1948, the previous 
New fields were discovered at Bonanza, Deep 


high year. 
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Creek (gas), Deer Creek. South Glenreck, Longs Creek (gas). 
Happy Springs, Meadow Creek, Nowood, Slick Creek, and 
Clareton. 
SOUTHERN STATES 

Alabama 

During 1950, there were 46 wells drilled in various parts 
of Alabama, 18 of which produced oil, one produced dry gas, 
and 27 were dry holes. There were 22 wildcat wells completed 
in 12 counties during the year, two of which resulted in new 
Most of the wildcat 
drilling was in the southwestern part of the state. Twenty-four 
field wells were completed with 17 producing oil and seven 
were dry holes. Alabama produced approximately 743,000 bbl 
of oil during the year, bringing the cumulative total for the 
state to approximately 2,733,000 bbl. 
Florida 

There were 14 wells drilled in Florida during 1950, eleven 
wildeats which were all dry field 
which was completed asa producer. The 1] wildcat wells were 


discoveries of reservoirs for the state. 


and three wells, one of 
drilled in nine counties and the three development wells were 
all in the Sunniland Field. There are now 12 wells producing 
in the lone Sunniland Field which produced approximately 
176,000 bbl of oil during 1950, bringing the cumulative total 


for the state and the field to approximately 1,553,000 bbl 
Georgia 

Georgia had five wells drilled during 1950, all of which 
holes. There two drilled in 


proved to be dry wells 


Seminole County and one each in Calhoun, Decatur and Heard 


were 


counties. These five wells bring Georgia’s total wildcats since 
1938 to 46 and all these wells have been dry. 


TENNESSEE 
An encouraging upsurge in drilling activity in Tennessee 
during the early months of 1950 was followed by a lull dur- 
ing the summer and early fall but as the year closed the deep 
tests started at the first of the year continued to drill and 
new tests are drilling or slated to start in the immediate 
future. Gas production during 1950 was slightly more than 
88 MMcf, a decided drop from production of 137 MMcf 
during 1949. 
TEXAS 
East and East Central Texas 
For the most part, new discoveries in the East Texas area 
were small and disappointing. Most 
in the Woodbine and Sub-Clarkville sands. Of the new sand 
the Woodbine Reservoir in the Pine Mills Field which 
previously has produced only from the Sub-Clarkville sand 


of the discoveries were 


is of some interest. Exploratory wells drilled in 1950 amounted 
to 19 o:l. six gas or gas condensate wells and 198 dry holes 
for a success percentage of 11.2 

North Texas 


Approximately 3,020 wells were completed in North Texas 
during 1950, of which about 1.660 were completed as oil or 
gas wells. Six hundred 


which 108 were completed 


seven wildcats were drilled, of 
as oil or gas wells. Production 
iately 53,300,000 bbl from some 23,000 
producing wells. Cumulative production in the area at the end 
of 1950 was in excess of 1,151,000.090 bbl. Some 60 new pools 
were reported during the vear, nor 


le of very 


and 


amounted to approx: 


of which appear at this 
time t great importance 
North Central Texas 

This area produced sligit!y more than 25.5 million bbl of 
oil during 1950 or a daily average of 70,000 bbl which 
approximatel, 3,000 B/D greater 1949. 


was 


At the end of 


than 
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there 
total cumulative production to January 1, 1951, was more than 
502,000,000 bbl. During the year some 2,040 wells were drilled 
of which about 1,000 were completed as oil or gas wells. 


the year were more than 8,600 producing wells and the 


Panhandle Area 

During 1950, the Texas Panhandle area saw some 890 wells 
drilled, of which 814 were completed as oil or gas wells. 
Iwenty-four of the wells were drilled as wildeats, resulting 
in the completion of one oil well and two gas wells. At the 
end of the year there were some 7,900 wells producing in the 
area. Production during 1950 amounted to some 33,120,000 
bbl or a little more less than 91,000 B/D and very nearly the 
same as during 1949. Cumulative production to the end of 
1950 approximated 696,000,000 bbl. 


South Texas 
Che Wilcox and Frio-Vicksburg trends received about equal 
n exploratory activities during 1950 in the South 


lexas area. The Portilla Field, San Patricio County, appears 
of major importance. Several Frio and 


attention 
to be a discovery 
Vicksburg sands produce oil in this field. Possibly the largest 
Wilcox discovery in this area in recent years is the Helen 
Gohlke Field in DeWitt and Victoria counties. Exploratory 
wells drilled during the year amounted to 114 oil, 56 gas or 
and 455 dry holes for a success percentage 


s condensate 


of 27.1. 
South Central Texas 

No discovery of apparent major importance was made in 
1950 in the South Central Texas area although the successful 
completion of a deep Edwards Limestone oil well in Atascosa 
County may bring renewed interest to the Charlotte-Imogene 
trend in that area. Wildcatting was scattered throughout 19 
of the 30 counties in the area. Of the exploratory wells drilled 
14 were oil four were gas or gas condensate and 189 
were dry holes for a success percentage of 8.7. 


wells, 


Upper Gulf Coast Area 


Although some fields of potential importance were discov- 
ered during 1950, the results of exploratory drilling in the 
Upper Gulf Coast were generally disappointing. Possibly the 
best of the year’s discoveries was the Nona Mills Field in 
Hardin County. Other discoveries of interest were the Rose 
City Field, Orange County: West Hamshire Field, Jefferson 
County (producing gas); and the Trinity Bay Field, Cham- 
bers County. Of the exploratory wells drilled 75 were oil 
wells, 11 were gas or gas condensate and 252 were dry holes. 
West Texas 

Slightly more than 5,000 wells were drilled in the West 
Texas area during 1950, of which approximately 4,130 were 
completed as oil or gas wells. About 670 were drilled as 
wildeats, of which some 116 were completed as oil or gas 
wells. The number of producing wells at the end of the year 
was approximately 30,300. Production amounted to slightly 
more than 284,500,000 bbl or an average of more than 779,- 
000 B/D. Cumulative production for the area to Jan. 1, 1951, 
was slightly more than 2,915,000,000 bbl. 

Activity Reef fields in Scurry, 
Kent and 1.300 wells had been 
completed by the end of the year. An important development 
was the Sweetie Peck discovery in Midland County. Some 19 
wells were producing at year’s end. Approximately 130 fields 
were discovered in West Texas during 1950 including the new 
production on 


the Canyon 
Borden counties. More than 


continued in 


brought into previously productive 


~*~ * * 


zones 
structures 
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THE CALCULATION OF WATER RESISTIVITIES FROM CHEMICAL ANALYSES 


H. F. DUNLAP, MEMBER AIME, AND R. R. HAWTHORNE, JUNIOR MEMSER AIME, ATLANTIC REFINING CO., DALLAS, 
TEX. 


ABSTRACT 


A method of calculating formation water resistivities from 
chemical analyses is presented which is somewhat faster and 
more accurate than previously described methods. For 26 
formation water samples taken from wells in the Texas Gulf 
Coast, the average deviation of calculated from measured 
water resistivities was 3.3 per cent, with an average bias of 
+1.1 per cent in the calculated values. This compares with 
average deviations of 4.6 per cent and 7.8 per cent, with 
corresponding average biases of +3.2 per cent and -6.3 per 
cent, for the two methods previously described in the litera- 
ture. All three methods are sufficiently accurate for most 
purposes; the main advantage of the method presented is its 
speed and convenience. 


DISCUSSION 


It is frequently desirable to calculate the formation water 


resistivity from the chemical analysis of the water, as for 


example, in electric log interpretation. Two methods have 
previously been given in the literature’* fog accomplishing 
this. Both methods require that the chemical analysis be given 
in the Palmer® system, as primary and secondary alkalinity 


Table I- 


Resistivity 


NaCl. 
Val. 
Vall. 
NaCl. 
NaCl. 
NaCl. 
NaCl. 


Total Equivalent ppm NaCl. 


Va ppm x 1.00 = N, Equivalent ppm 
ppmx .95 = N, Equivalent ppm 
ppm x 2.00 = NV, Equivalent ppm 
ppm x 1.00 = NV, Equivalent ppm 
ppmx .50 = ', Equivalent ppm 
ppmx .27 = N, Equivalent ppm 
ppm x 1.26 = N, Equivalent ppm 

= Ni 
and salinity. This conversion of the values of the ppm of 
the various ions to primary and secondary salinity and alka- 
linity is time consuming, and many laboratories report only 
the concentrations of the various ions in parts per million 
However, all laboratories report the individual ionic concen 
trations, whether or not the Palmer system of reporting is 
also used. Hence it was thought that a method that would 
permit resistivity calculations directly from these data would 
be useful. 

In the method developed, the various ions reported in the 
chemical analysis are converted to an equivalent amount of 
VaCl. The sum of the various contributions to the NaCl (ac 
tual or equivalent) can then be used in conjunction with a 
plot of NaCl concentration versus resistivity to obtain the 
resistivity of the solution in question. The conversion factors 
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Conversion Factors for Calculating Water 
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FIG. 1 — NaC] CONCENTRATION VS. RESISTIVITY 


u-ed in this process are shown in Table I and the resistivity 
vs. NaCl saturation plot (1.C.T. data) in Fig. 1. An example 
of a typical calculation is shown in Table I, together with 
the measured resistivity. 

These conversion factors for the various ions were obtained 
by measuring the resistivities of a number of simple solutions 
containing only two types of ions, one of which was always 
either Na+ or Cl-. Although for dilute solutions of most ions 


Continued on Page 7, Section 2 


Table Il — Example of Method Applied to a Frio 
Water from the Aransas Field, San Patricio County, 
Texas 

Equivalent 


Concentration 
of NaCl (PPM) 


21,899 
2.325 


Conversion 
Factor 
Va (calculated) 21,899 1.00 
C1 2.449 2 
Ve 382 2.00 764 
cl 38.292 1.00 38.292 
SO, 381 50 19] 
HCO 1.098 27 296 
CO 0 1.26 0 


Ionic 


Concentration 
Constituent PPM 


63,767 
.1260M 
1310M 


Total Equivalent NaCl concentration = 
Calculated Resistivity (using Fig. 1) 
Measured Resistivity 

Error 3.8 per cent 
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F. E. Hatfield, manager of exploita- 
tion of the Fort Worth Division of the 
Gulf Oil Corp., outlined last month 
highlights of a recent trip through the 
Middle East oil fields. His remarks 
were directed to members of the North 
+exas Section in meeting at Fort Worth. 

Known reserves of the Middle East 
are, Hatfield pointed out, enormous. 
Estimates are in the order of 40 billion 
or more bbl proven reserves and poten- 
tial reserves of 100 to 150 billion bbl. 
or about one-half the world’s supply. 
For comparison, United States reserves 
are estimated at 25 to 30 billion bbl 
proven and 50 billion bbl potential. 
Current Persian Gulf da‘ly production, 
averaging 5,000 bbl per well daily. is 
more than 1,800,000 bbl, or 18 per cent 
of the world’s total of 10,000,000 B/D. 
Some time next year. if 500,000 B/D 
are taken from Kuwait, as it is rumored, 
production may reach 2 million B/D 
and in 1952 and 1953, with Middle East 
pipe line completed, production could 
reach 2.5 million B/D. 


The largest and oldest producing 
company is Anglo-Iranian. Although oil 
seeps in the Middle East have been 
known since Biblical times, the fore- 
runner of Anglo-Iranian. the D’arcy 
Petroleum Co., discovered the first oil 
in the Persian Gulf in the Sulaiman 
Field in 1908. The concession taken by 
D’arey in 1902 cost $20,000. Anglo- 
Iranian’s daily production in Iran is 
about 625,000 bbl and it has a refining 
capacity at Abadan, Iran. of 550,000 
B/D. 

Commercial oil was next found in 
Iraq by the Iraq Oil Co. in 1927, which 
resulted in building the first trans- 
desert pipe line in 1934, a double 12-in.. 
to the Mediterranean Sea. In 1948, one 
of these lines which went to Haifa in 
Palestine was cui off by the Arab Iraq 
government and this line is still out of 
service. Later a 16-in. line was laid 
from Kirkuk to Tripoli in the Lebanon 
and Bechtel is now beg nning another 
30-in. line from Kirkuk to Banias. 
Syria. a distance of about 560 miles. 
which will handle around 275,000 B/D 
and capacity of all Iraq Petroleum Co. 
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Oil in Kuwait and the Persian Gulf 


Reported By 
John Clayton Jacobs, Jr. 
North Texas Section 


lines would then be more than 500,000 
B/D in 1952 and 1953. Present produc 
tion is just under 170,000 B/D. 

The second largest producing area is 
Aramco’s series of fields in Saudi Ara- 
ba near the Gulf discovered, just 
across a narrow strait from Bahrein 
Island, in 1937. Since then, seven fields 
have been found in all, the largest. 
Abquaiq. producing around 525,000 
B/D from 50 wells. Within the last 
year and a half, a 140,000-bbl refinery 
has been built by Aramco at Ras Tan- 
ura on the Gulf and the new 30-in., 
1.000-mile pipe line from Abquaigq 
Field to the Mediterranean at the an- 
cient city of Sidon, Lebanon, was com- 
pleted December 1 at a cost of $200.- 
000.000. These facilities will handle 
315,000 B/D and save about $.18 per 
barrel toll through the Suez Canal as 
well as tanker costs of $.65 to $.70 per 
barrel for the 3,200-mile trip. Aramco’s 
production is now more than 640,000 
B/D. 

The Bahrein Island Field, about 25 
miles from the Arabian mainland. was 
discovered in 1932 and produces some- 
thing less than 30,000 B/D. 

Now something about Kuwait and 
the Burgan Field, the world’s largest 
oil field. After a geological and geo- 
physical survey by the joint owners of 
the concession — Gulf and Anglo-Iran- 
ian— which disclosed some asphalt 
seeps and evidence of structure, the 
first well begun in May, 1936. was dry 
at about 8.000 ft. Oil was encountered 
in Burgan No. 1 in February, 1938. 
Drilling went on slowly until nine wells 
were completed in 1942 when, due to 
the war, it was thought advisable to 
plug them off and shut everything down. 
In late 1945, the wells were cleaned 
out, drilling of additional wells began 
and the first tanker was loaded in 
June, 1946. At present about 100 wells 
are completed, but only about 65 are 
used as producers because of the fact 
that the highest gravity oil, around 31 
comes from the highest part of the 
structure. 

The pay is sand. Middle Cretaceous. 
lying in three to four zones between 
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ENGINEERING DEVELOPMENTS IN LOCAL SECTIONS 


3,500 ft and 4,800 ft with as much as 
550 ft of net thickness. The wells are 
practically all dual zone. completed in 
the conventional way. The reservoir 
contains approximately 60,000 acres 
and is about seven by 14 miles in area. 
Spacing is one well to 600 acres. More 
than 240,000,000 bbl of oil have been 
withdrawn with no discernible drop in 
bottom hole pressure. 

A brief look at the surface installa- 
tions which are typical of the large 
Middle East fields reveals that there are 
five gathering stations which cost, 
erected, about $1,600,000 each, han- 
dling from 40,000 to 175,000 B/D. 
Horizontal separators are used with 
four to seven stages of separation. 
Tkere are two 50,000-bbl receiving 
tanks and a gas engine-driven pumping 
installation to deliver oil to the tank 
farm about 11 miles away. At the tank 
farm there are 21 big storage tanks 
where the oil is gauged and fiscalized 
for government royalty purposes. 

The oil is then gravitated five miles 
to the oil loading dock on the Persian 
Gulf. This installation can load, under 
ideal circumstances, as much as 700,000 
B/D from the new oil and cargo jetty, 
finished late in 1949. On the coast near 
the jetty there is also a 25,000-bbl re- 
finery which was placed in operation in 
October, 1949, producing kerosene, gas- 
oline, diesel and bunker fuel for use in 
the country and fuel for ships calling 
at Mena El Ahmadi. A steam electric 
station and sea water vacuum distilla- 
tion plant, to produce 15,000 kilowatts 
and turn out 10,000 B/D fresh water 
from sea water, is now completed. 

The Ahmadi headquarters camp is 
up on the hill near the tank farm. 
There are about 260 permanent houses 
for British and Americans besides sev- 
eral hundred temporary houses, gen- 
eral office, shops, warehouse, cooling 
and heating plant, mess halls, barracks, 
ete. 

Some 60 members and guests of the 
North Texas Section heard Hatfield's 
description of Middle East operations 
R. S. Christie introduced the speaker. 
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Operating Conditions in South 


William H. Sargent, recently re- 
turned from an extended tour through 
Seuth America in his capacity as head 
of foreign sales for Sargent Engineer- 
ing Corp., spoke to members of the 
Junior Section of the Pacific Petroleum 
Chapter in Downey, Calif., last month. 

Sargent prefaced his discussion, 
which depicted his travels through 
South America, with some general com- 
ments on what the “gringo” can expect. 
He de:cribed the Latin American coun- 
tries as a land of many extremes. On 
the one hand is found the most modern 
garbage disposal system in the world, 
while, at the same time, there is no 
sewage system, and sewage is allowed 
to run down the streets. At Lake Mara- 
caibo, Venezuela, intensive modern oil 
drilling and production proceeds di- 
rectly over the water, while the natives 
live in primitive huts along the shore. 

This Maracaibo development is tre- 
mendously impressive. In the beginning. 
one company leased the shore rights 
and began drilling. A second company 
came in on the offshore rights, and a 
third leased the center of the lake. To- 
day all three operators have producing 
wells, some as far as 22 miles off shore. 
The problem of establishing and main- 
taining a location on the lake is severe 
due to the high waves which can engulf 
and capsize boats, and because of the 
silty bottom of the lake which seem- 
ingly swallows up foundations like 
quick sand. It is common practice to 
sink heavy cement pillars for corners 
of the derrick floor some 300 ft into 
the silty lake bottom without ever 
reaching bed rock. These pillars are 
allowed to settle in by placing weights 
upon them, as it has proven impossible 
to jet or drive them. In some cases, the 
wells are drilled in 200 ft of water 
which necessitates a total depth of 500 
ft for the cement pillars. 

Sargent noted that a majority of oil 
company technical personnel in South 
America come from the Gulf Coast or 
Europe. The opportunity for capable 
technical men is excellent, but Sargent 
warned against the idea of getting rich 
quick. He noted that the South Ameri- 


March, 1951 


Reported By 


W. G. Paul, Jr. 


Pacific Junior Group 


can may appear ignorant, but that ac- 
tually ke is far from it. The gen- 
eral tendency is to welcome anyone who 
has something constructive to offer, but 
to eliminate the man who is only inter- 
ested in making a dollar at the expense 
of the native. 

For those contemplating South Amer- 
ican visits, Sargent advises a_ large 
wardrobe of nylon clothing. While it 
has the disadvantage of retaining heat, 
it has a major advantage in that it can 
be washed by hand. Native washings 
return clothes either green or brown, 
depending upon whe.ner there is moss 
or mud in the local river. In many 
cases, laundry may not be returned at 
all due to the native affinity for white 
clothing. Since the native South Amer- 
ican is always presumed right in any 
dispute with a foreigner there is little 
recourse. Khaki clothing is always ac- 
ceptable in the field and it is very sel 
dom that a business suit would be 
necessary outside the major cities. 

To avoid the tropical diseases which 
are prevalent, Sargent advised against 
eating anything touched by the natives. 
Eat only fruits or vegetables which can 
be peeled and avoid water and fresh 
milk. Every liquid should be boiled and 
ice should never be added. Sargent 


America 


noted that in the camps of the major 
companies, food is flown in to avoid 
contamination. If a man will eat only 
in the company camps or in the best 
known eating places, he should be able 
to survive a South American visit with- 
out difficulty. 

Sargent devoted some time to the 
city of Caracas, the capitol of Vene- 
zuela. This city is the center of foreign 
oil operation in South America and 
virtually every oil company has a head 
office there. The problem of a sales or- 
ganization at Caracas is two-fold; the 
high cost of maintaining a man and the 
beauty of the city itself, which deters 
quite a few from traveling back to 
where the actual producing operation 
is going on. Sargent estimated that it 
would cost between $40,000 and $50,000 
per year to station a representative in 
Caracas. 

Sargent concluded his talk with a 
series of pictures on Trinidad. He de- 
scribed in detail the pitch lake from 
which the highest grade roofing tar 
in the world is produced. The pitch is 
pried in chunks from the lake’s surface 
and trucked to a nearby refinery. After 
six months, the hole will have refilled 
itself and the same operation can be 
repeated again and again. x * * 


PRODUCING PROBLEMS in South America are compared with domestic ones 
by W. G. Paul, Jr., left, and William H. Sargent of Sargent Engineering. 
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Utilization of Connate Water Data 


Che Study Group program announced 
by the North Texas Section in the 
January issue of Petroteum TecuNo.- 
ocy has been inaugurated and is being 
enthusiastically received by those par- 
ticipating in the meetings. The meetings 
of the Dallas Group, which are held 
on the second and fourth Wednesdays 
of each month at 7:30 p.m., are at- 
tended by 40 to 60 engineers and geol- 
ogists. Marshall Fagin of the Petroleum 
Engineer is in charge of the Dallas 
Group. The Fort Worth Group, which 
meets on the first and third Wedne:- 
days of each month at 7:30 p.m. at the 
Fort Worth Public Library, has an 
attendance of 20 to 30. Paul Wiley of 
The Texas Co. is in charge of the Fort 
Worth Group. 

The first meeting of the Fort Worth 
Group was built around a discussion 
of the “Derivation of the Material Bal- 
ance Equation.” This discussion was 
led by R. M. Liebrock of the Stanolind 
Oil and Gas Co. At the second meeting 
Liebrock led a discussion on “The Ap- 
plication of the Material Balance Equa- 
tion to a Field Problem.” At the third 
meeting discussions of “Field Data for 
Reservoir Engineering Studies” and 
“PVT Analyses” will be held under the 


A Report From 
The North Texas Study Group 


direction of Dave Meisenheimer of The 
exas Co. and Ray Ard of the Gulf Co. 
The third meeting was held February 
21. Those attending the meetings have 
been polled as to their choice of pro 
gram and a schedule of future pro 
grams is to be announced shortly. 

The Dallas Group chose “Connate 
Water” as an initial discussion topic. 
At the first meeting the discussion, 
which was led by Bill McGuire of the 
Atlantic Refining Co., covered the “Cen- 
trifuge Method of Determining Connate 
Water.” At the second meeting the dis- 
cussion dealt with the definition and 
determination of connate water, and the 
utilization of connate water data. 

Atlantic Refining’s R. L. Slobod was 
the discussion leader. A summary of 
the discussion follows. 

While it was generally agreed that 
there is a great need for data which 
would enable the engineer to arrive at 
a figure fer connate water for any 
given reservoir, there was considerable 
divergence of opinion as to the exact 
definition of the terms connate water, 
interstitial water, residual water and 
irreducible water. Not only does this 
divergence exist within the petroleum 
engineering field, but also apparenth 


among geologists. Whatever term is 
applied to the result of a laboratory 
test or series of tests, there must be 
a clear understanding on the part of 
everyone involved, from the laboratory 
technician making the tests on up to 
the reservoir engineer utilizing the data. 
of the entire scope and aim of the en- 
deavor. The limitations of sampling, of 
method, of _ statistical 
must be considered in evaluating the 
final result. The blind acceptance of a 
“per cent connate water” or “average 
connate water” figure should be avoided. 


analysis — all 


The various laboratory methods of 
obtaining capillary pressure curves and 
from them the irreducible water satu- 
ration of core samples were discussed. 
The capillary disc method, with its in- 
herent disadvantages, seems to be the 
most common method used in the indus- 
try. The concensus seems to be that it 
is the basic method to rely upon. Other 
techniques, such as the centrifuge 
method and the evaporation method 
(which gives not the full capillary pres- 
sure curve but only its end point) have 
not been generally accepted throughout 
the industry. 

The use of electric logs for the deter- 
mination of connate water was dis- 
cussed and, while opinions varied wide- 
ly as to the reliability of data so ob- 
tained. it was generally agreed that the 
most valuable feature of these deter- 
minations lies in the fact that the meas- 
urements are made on the reservoir 
itself. 

In regard to the upper and lower 
limits of reliable connate water deter- 
minations it was brought out that fig- 
ures as low as 10 per cent and as high 
as 50 per cent have been generally 
accepted, but when the percentages go 
much higher or lower than these limits 
the tendency is to doubt their reli- 
ability. 

Correlation of laboratory results with 
direct water determinations made on 
cores drilled with oil base mud is con- 
sidered generally to be the most re- 


liable method of checking the various 
laboratory techniques. Regrettably, only 
a relatively small amount of such .cor- 
relation has been presented in the lit- 
erature. x * * 


FEAST TEXAS SECTION OFFICERS — The new officers and directors of the 
ast Texas section for 1951 are, front row, L to R, W. R. Mays, vice-chairman; 
W. E. LaRoche, secretary-treasurer; W. W. Leonard, chairman; and P. J. Lehn- 
hard, director. In the back row, L to R, are Directors H. L. Long, J. A. Walker, 
Ww. * Morris, P. E. DesJardins, W. M. Jones, Bert E. Crowder, and B. F. Patter- 
son, Jr. 
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Elements of Rock Bit Design 


Reported By 


J. E. Kastrop 
Gulf Coast Section 


One of the most stimulating talks in 
several months was presented at the 
February meeting of the Gulf Coast 
Section at Houston, Tex. Henry Woods, 
assistant chief research engineer for 
Hughes Tool Co., Houston, gave a prac- 
tical discussion of the “Elements of 
Rock Bit Design.” Through the use of 
slides and demonstrations, he showed 
the possibilities and limitations in the 
design of rock bits. One of these limi- 
tations, he said, was the size of the 
hole which the bit was intended to drill. 
To enlarge one part means reducing 
other parts. 

For emphasis, a rock bit was de- 
signed for three different formations. 
Based on knowledge of bit design and 
formation characteristics, a bit that will 
produce maximum results can be made. 
In quartzite with a compressive strength 
of 68,000 psi, which is considerably 
harder than thove generally encountered 
in oil well drilling, it is essential that 
a crushing action be obtained from the 
teeth. A drag or scraping action would 
result in breaking the teeth from the 


LEON SINGLETON 


Leon Singleton, Houston Oil Field 
Material Co., is the Gulf Coast Sec- 
tion’s new JOURNAL Secretary. He has 
been closely associated with the in- 
dustry, particularly in oil field equip- 
ment and services since his discharge 
from the U. S. Army. 
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cone, or they will wear more than the 
rock. In drilling quartzite, a pressure of 
at least 68,000 psi must be exerted by 
the teeth in contact with the formation 
if any fracturing is accomplished. The 
number, size and shape of the teeth on 
such a bit are determined by the load 
which they must carry. 

In a formation of intermediate hard- 
ness, Woods showed that it was best to 
use a combination gouging and scrap 
ing action. The teeth must be designed 
to penetrate the formation and give a 
twist or jerk that will break the chip 
loose. This twist action is accomplished 
by swelling the apex of the cone so 
that it is not true rolling, and by off- 
setting the cone somewhat so that the 
bearing about which the cone turns is 
slightly ahead of center. 

The teeth are 
apart so that the formation 1s not too 
likely to pack between them. To break 
up the ridges made by the teeth on the 
heel row, different numbers of teeth 


spaced far enough 


are used on each cone. 

A third and perhaps extreme type of 
formation drilled is that of soft clay. 
In removing this type material, not 


HENRY WOODS 


only must the teeth penetrate the for- 
mation, but they must dig it out with 
a shovel action. To drill such forma- 
tions, the teeth need be widely spaced 
to prevent balling up the bit. The use 
of interrupted heel teeth is desirable. 
Teeth should be as slim as possible 
without giving too much breakage. The 
maximum amount of tungsten carbide 
is put on the teeth to make it drill in 
zones streaked with intermittent soft 
and harder formations. 

In designing a bit, Woods pointed 
out, it is necessary to know whether it 
is better to concentrate on increasing 
the rate of penetration or making more 

(Continued on Page 6, Section 2) 


Thin Oil Column Reservoirs 
Reported By 
J. W. Walker 
Gulf Coast Study Group 


Problems of producing reservoirs 
with thin oil columns were discussed 
by R. W. Holman of The Texas Co., 
following a brief introduction to the 
subject by Norman Clark of the Hum- 
ble Oil and Refining Co., before a meet- 
ing last month of the Gulf Coast Sec- 
tion Study Group. Clark showed that 
the thin-oil column reservoir is a spe- 
cial case of the combination drive res- 
ervoir. Reservoir mechanics and analy- 
sis are the same in both cases; however. 
individual well problems assume more 
importance in exploiting thin-oil col- 
umn reservoirs than in the normal com- 
bination reservoir 

Individual well problems were out- 
lined by Holman who presented a de- 
tailed discussion of the theory of con 
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ing of oil and gas. The effect of sand 
penetration, completion methods, den- 
sity differences, and rate of production 
were discussed. Methods of reducing 
gas/oil ratios by mechanical means 
were reviewed. The economical aspects 
of exploiting thin oil columns encour- 
aged considerable group discussion. 
The problems involved in obtaining 
accurate PVT data were discussed by 
J. W. Walker of Shell Oil Co. Discus- 
sion revealed that in some cases the 
gas cap of a thin-oil column reservoir 
has been produced for some time before 
the presence of a thin oil column is 
detected. Methods of obtaining repre- 
sentative analyses both of the reservoir 
liquids and of the gas-cap gas were 
reveiwed. xe 
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University of Tulsa Student Activities 


The Student Chapter of the AIME 
at the University of Tulsa has been 
very active this past semester. In addi- 
tion to the sponsorship of lectures by 
prominent men of the industry, a num- 
ber of field trips were held to plants 
in and around Tulsa. A freshman ori- 
entation program was also held in No- 
vember to acquaint the entering petro- 
leum engineering majors with the aims 
and activities of the AIME. 

On November 14, Warren B. Davis 
of Gulf Oil Co. spoke to the chapter 
on the Dornick Hills-Springer Sand Res- 
ervoir, Velma Pool, Okla. He supple- 
mented his talk with slides and the 
talk was well received by a large group 
of students. 

The following week Lloyd Elkins, 
director of production 
Stanolind Oil and Gas Co. gave an in- 
structive lecture on “Research; Past, 
Present and Future.” In this talk. El- 
kins defined research first, then went 
on to a discussion of the development 
of research methods. 

Following the membership drive in 
September and October, in which 50 
new memberships were received, a num- 


research for 


ber of field trips were taken to various 
places of interest. These trips included 
tours of the National Tank Co., Mid- 
Continent Refinery, Bethlehem Supply 
Co., and an all-day tour through South- 
ern Oklahoma oil operations in the 
Golden Trend Area. These trips were 





Meetings of Interest 

Spring Meeting of the Interstate 
Oil Compact Commission will be 
teld in the La Fonda Hotel, Santa 
Fe, N. M., April 26-28. 

It will be impossible to place all 
those who plan to attend the meet- 
ing in the headquarters hotel, but 
the Commission’s headquarters of- 
fice, Box 3127, State Capitol, Okla- 
homa City 5, Okla., through which 
all reservations are being handled. 
has arranged for accommodations 
elsewhere. Committee meetings will 
be held in the headquarters hotel 
luring the Santa Fe meeting; how- 
ever, general sessions will be held in 
the Capitol Auditorium on the Capi- 
tol grounds. 
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By Jack Combes 


University Student Chapter 


all under the sponsorship of Professor 
A. W. Walker, faculty advisor. 

Near the close of the semester, a 
joint dinner meeting of the Mid-Conti- 
nent Section and the University of 
Tulsa was held in the new Student Ac- 
tivities building on the campus. Speak- 
er of the evening was Thomas West of 


A 


Blanco Oil Co., San Antonio, Tex., who 
spoke on “A Gravel Pack Completion 
for Exclusion of Gas and Water.” 
New officers elected for the 1950-51 
school year were Jack Combes, presi- 
dent; Jerry Stafford, vice-president; 
Tom Carlson, secretary: and Dale Rub- 
inson, treasurer. x * * 
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UNIVERSITY OF TULSA STUDENT CHAPTER — 1951 officers and faculty 
advisers of the Student Chapter at the University of Tulsa are, front row, L to R, 
E. E. Barberii, instructor; Thomas Carlson, secretary; Jerry Stafford, vice-presi- 
dent; Warren Davis, Gulf Oil Corp., speaker; Dale Robinson, treasurer; Jack 
Combes, president; and A. W. Walker, faculty advisor. 


Bottom-Hole Pressure Build-Up 


Southwest Texas Section 
Report 


William Hurst, petroleum consultant, 
Houston, Tex., presented to the South- 
west Texas Section a talk titled “Sig- 
nificance of Bottom-Hole Pressure 
Build-up with Respect to Well Com- 
pletion” before the February meeting 
of the section in Corpus Christi, Tex. 

The speaker dealt with the subject 
pertaining to the physical concepts as- 
sociated with bottom-hole 
build-up. In this respect the effects of 
well unloading and loading in its rela- 
tionship to pressure draw-down and 
build-up were discussed. This likewise 
included the interpretation of produc- 
tivity index with respect to well com- 
pletion, and its significance as estab- 
lished from bottom-hole pressure analy- 
sis. The talk was heizhtened by an in- 


pressure 
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teresting discussion which followed, 
bearing upon the observations of fac- 
tual field examples. xe 





The 1951 officers of the Petroleum 
Engineers Club (AIME Student 
Chapter) at the University of Okla- 
homa are: 

E. J. Dickinson, president; How- 
ard M. Perdue, vice-president; S. D. 
Carnahan, secretary; F. A. Harris, 
treasurer; A. J. Hankinson, chair- 
man of entertainment; W. T. Mas- 
sey, manager of publicity; Glen and 
Omar Emerick, co-chairmen of the 
open house comittee; and P. G. Von 
Tungeln, St. Pat’s Council repre- 
sentative. 
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1951 Petroleum Branch Chairman Richard W. 
French, left, beams broadly at John E. Sherborne, 
1950 Chairman, just before the latter relinquished 
the chairmanship at the Petroleum Branch Ban- 
quet in St. Louis. Chairman French had at the time 
already conducted sessions of the Branch Executive 
Committee, but the formal changeover did not take 
place until the night of the banquet. 
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. Louis, Mo. 


; 


it, 





i | 
A portion of a technical session audience listen intently to Lloyd Elkins of Stanolind Oil and Gas Co., Tulsa, Okla., i 
the presentation of one of the 10 papers presented at the right, has something light to say at the Petroleum Branch 7% 
four-day meeting. John C. Calhoun of Pennsylvania State Col banquet. Some 150 members and guests turned out for the 
lege, center middleground, thumbs through his preprint dur informal affair held in the Statler Hotel. Texas Railroad Com- 
ing the delivery. missioner William J. Murray was the speaker for the evening. 


HE 171st General Meeting of the AIME, held 
T last month in St. Louis, Mo., has been ranked 
with the most successful staged by the Institute. The 
practice, established in recent years, of holding the 
Annual Meeting on alternate years outside New 
York City has again been proved a wise move. All 
Branches and Divisions occupied session rooms in 
four hotels in the Mid-west city. The Statler Hotel. 
with but four elevators, shuttled Metals and Petro- 
leum people to and from the 16th and 17th floors, 
respectively, without excessive delay. Meetings were 
also conducted in the Jefferson, Lennox and DeSoto 
hetels. The loca] Lions Club interrupted one of the 
Petroleum Branch sessions to claim the room for 
the’r weekly luncheon site. Most sessions were held 
before overflow attendance. At right is a view of 
the regi-tration desk in the Statler lobby. 
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DISCUSSION OF THIS AND ALL FOLLOWING TECHNICAL PAPERS IS INVITED 
Discussion in writing (3 copies) may be sent to the Editor, Journal of Petroleum Technology, 601 Continental Building, 
Dallas 1, Texas, and will be considered for publication in the Transactions volume Petroleum Development and Technology. 
Discussion will close December 1, 1951. Any discussion offered thereafter should be in the form of a new paper. 











CAPILLARY PRESSURE INVESTIGATIONS 
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METHOD 
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ABSTRACT 


his paper presents results of static capillary pressure meas- 
urements made both by the restored-state and by the mercury- 
injection methods, and of dynamic capillary pressure meas- 
urements made by the Hassler technique. 

[he results indicate a variation from sample to sample of 
the factor used to convert capillary pressure obtained by mer- 
cury injection to those applicable to a water-gas system. Static 
capillary the restored-state method 
were found to be in close agreement with dynamic capillary 
pressures, indicating that data obtained from static measure- 
ments may be utilized with confidence in the solution of dy- 
namic problems of fluid flow. 


pressures measured by 


Experimental data are presented on the use of “J-curves” 
for the correlation of capillary pressure-saturation relations 
for a number of cores from a particular geologic formation. 
It is shown that the correlation was improved in some cases 
by restricting it to certain lithologic classifications. 


INTRODUCTION 


The Meaning and Importance of Capillary Pressure 

The coexistence of two or more immiscible fluids within the 
voids of a porous medium, such as reservoir rock, gives rise 
to capillary forces. Because of the interfacial tension existing 
it the boundary between two immiscible fluids in a pore space, 
the interface is curved, and there is a pressure difference 
across the interface. This pressure difference is termed the 
capillary pressure. 

The magnitude of the capillary pressure between two immis- 
cible fluids which together fill a porous medium depends upon 
several factors: the textural properties of the medium, the 
wettability of the medium, the interfacial tensicn between the 
fluids, the respective saturations of the fluids, and the manner 
in which these saturations are attained. For a given pair of 
immiscible fluids in a particular sample of reservoir rock, the 
capillary pressure is a unique function of the fluid saturation, 
provided that the saturation of the fluid which wets the rock 
has previously been decreased unidirectionally from an ini- 
tially complete saturation. Similarly, for a given pair of fluids 
it a particular saturation which has been achieved by a 
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50. Paper presented at the Petroleum Branch Fall Meeting in New 
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CAPILLARY PRESSURE INVESTIGATIONS 


GAS (nother technique for determination of capillary pressure 
NLET has been proposed by Purcell.’ In this technique, mercury, the 
INL non-wetting phase, is injected under pressure into an evacu- 
J ated core sample. The mercury vapor, together with any resid- 

ual gas present, corresponds to the wetting phase. The mer- 
cury-injection method has the advantage of rapidity of opera- 
tion, since equilibrium is attained in a short time after the 
application of each pressure increment. In contrast, equilib- 
rium at each step is attained in the restored-state technique 
only after a matter of days. 








. 











Both the restored-state and mercury-injection techniques are 
static methods. Dynamic capillary pressures may be deter- 
mined by a modification of the technique developed by Hass- 
ler‘ for measuring relative permeability. In this method both 
phases are permitted to flow through the core sample, with 
control being exerted over the pressure difference existing 
between them. At equilibrium the two phases continue to flow. 
but their respective saturations do not change. 


























Leverett’® and Rose and Bruce* have pointed out the hetero- 
geneity of reservoir rock and the dependence of capillary pres- 
sure on textural properties of the rock. These authors have 

















capillary pressure to saturation with textural properties of 
rock, in order to extend the applicability of capillary pressure 
TO measurements. This method of correlation involves the use of 
the “J-curve.” which is discussed in later sections. 
ATMOSPHERE PORCELAIN 
PLATE 





~ L also proposed a method of correlating the relationship of 
OIL INLET 


PROBLEMS INVESTIGATED IN THIS WORK 


OIL Of the numerous problems existing in the measurement and 


BURETTE ©) application of capillary pressures, three were investigated 


during the course of this work. 





FIG. 2—DYNAMIC CAPILLARY PRESSURE APPARATUS (HASSLER’S 
PRINCIPLE). 
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previous decrease of saturation of the wetting fluid, the capil- 
lary pressure depends upon the textural properties and wet- © Restored _ state 


tability of the reservoir rock. 

Although the absolute magnitude of the capillary pressures 
in most petroleum reservoirs is usually not large, the effects 
of these capillary pressures are extremely important. Together 
with gravity they control the original distribution of the fluid 
saturations within the reservoir, particularly the distribution 
of connate water. By virtue of their effect on the shapes of the 
fluid interfaces within the pore spaces, they control in large 
measure the relative freedom of movement of fluids present in 
the reservoir, and are important factors in influencing the 
behavior and distribution of fluids in the production processes. 


e Mercury injection 


— oe 


SANDSTONE CORE 
Porosity —28.1 % 
Permeability — 1.43 darcys 
Factor — 7.5 
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Previous Investigations on Capillary Pressure 


Methods of measuring capillary pressures of fluids in small 
core samples, and of determining the relation between these 
capillary pressures and fluid saturations have been given by 
Leverett’ and by Bruce and Welge.* These authors have pro- 
posed the use of the “restored-state” technique, in which the 
core sample, saturated with a wetting fluid, is placed on a 
porous diaphragm permeable to the wetting phase. By the 100 80 60 40 
application of a known pressure to the non-wetting phase 
which is confined above the diaphragm, a portion of the wet- 
ting phase is expelled from the core sample. By the use of LIQUID SATURATION PER CENT 
successively higher pressure increments in this procedure a 
relation between the capillary pressure and fluid saturations FIG. 3—CAPILLARY PRESSURES BY RESTORED-STATE AND MERCURY. 
may be obtained. INJECTION METHODS 
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FIG. 4— CAPILLARY PRESSURES BY RESTORED-STATE AND MERCURY 
INJECTION METHODS. 


Correlation of Results of Restored-State 
and Mercury-Injection Methods 


\lthough the restored-state technique is probably the most 
commonly used method and is generally considered to be the 
most reliable, the mercury-injection technique has the practi- 
cal advantage of being much more rapid. Since fluids similar 
to those encountered in the reservoirs may be used in the 
restored-state determination, the results are usually more di 
rectly applicable to reservoir problems; the results of the 
mercury-injection method, on the other hand, must be con- 
verted by means of a conversion factor to yield those capillary 
pressures which would be obtained with reservoir fluids. To 
convert capillary pressures obtained by the mercury-injection 
method to those obtained by the restored-state technique using 
water and gas, Purcell’ has suggested a constant conversion 
factor which is independent of the character of the reservoir 
rock. 

In view of the complicated character of the material com 
prising most reservoir rock, and in view of the limited amount 
of work which has been reported on capillary pressure meas 
urements by the mercury-injection technique, it was deemed 
advisable to investigate the two methods on a comparative 
basis. In particular, it was thought pertinent to examine the 
constancy of the conversion factor for cores from various types 
of formations. 


Comparison of Static and Dynamic 
Capillary Pressures 


Most of the determinations of capillary pressure-saturation 
relations have been made under static conditions. Although it 
has been presumed that capillary pressures have an important 
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effect on oil recovery by causing fluid movements to take place 
when a reservoir is flushed by gas or water, it has not been 
known whether, when the fluids are in motion, the capillary 
pressures are the same as when the fluids are at rest. Several 
papers’”** presented in recent years have proposed methods 
of applying static capillary pressure data to the solution of 
dynamic problems of flow. Because of the fundamental nature 
of the differences between the static and dynamic systems, the 
validity of such application is subject to some question. 

In investigating the various problems connected with capil- 
lary pressure determinations, it was particularly desired to 
compare capillary pressures measured under conditions of 
static equilibrium with those measured under conditions of 
dynamic equilibrium during flow. 


Correlation of Capillary Properties 
of Rock with Rock Texture 


The method of correlating capillary pressure-saturation 
relations with permeability, porosity, and wettability of the 
rock which has been proposed by Leverett’® and by Rose and 
Bruce’ requires the use of a dimensionless “capillary pressure- 
function” group. A plot of this dimensionless group against 
wetting phase saturation gives the curve generally known as 
the capillary pressure-function curve or “J-curve.” It has been 
suggested that this curve is unique for core samples from a 
particular geologic formation. If this be true, it is apparent 
that the J-curve, once derived for a particular formation, is 
of signal importance either in permitting the computation of 
capillary pressure-saturation relations from known values of 
permeability and porosity, or in permitting the calculation of 
the permeability from known values of capillary pressure and 
porosity. Adequate information on the subject has been lack- 
ing, however, and the need for additional data prompted an 
attempt to confirm the correlation. 


EXPERIMENTAL PROCEDURE 


Static Capillary Pressure Measurements 


Static capillary pressure-liquid saturation relations for core 
samples were obtained by use of both the restored-state tech- 
nique and the mercury-injection method. The core samples 
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FIG. 5— CAPILLARY PRESSURES BY STATIC AND DYNAMIC METHODS 
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used in these measurements consisted of both limestone and 
sandstone samples from producing zones. 


Restored-State Technique 


The displacement cell employed in measurements made by 
the restored-state technique was of the conventional type, the 
essential feature of which is the porous porcelain diaphragm. 
This diaphragm is so chosen that when it is saturated with 
the fluid to be displaced from the core, it is permeable only 
to that fluid. When capillary contact is established between 
the diaphragm and the saturated core, the application of a 
suitably chosen constant pressure differential across the dia- 
phragm will cause fluid to be displaced from the core, through 
the diaphragm, to the extent that sufficient fluid remains to 
satisfy capillary forces. The pressure differential applied 
across the porous diaphragm is the capillary pressure, be 
cause, by definition, capillary pressure is the difference be 
tween the pressures of the wetting and non-wetting phases. The 
quantity of fluid remaining in the core is determined by 
weighing. 

All of the restored-state results presented in this paper were 
obtained under conditions of decreasing wetting phase satu- 
ration. 


Mercury-Injection Method 


The apparatus used to determine capillary pressure curves 
by the mercury-injection method was patterned after that 
described by Purcell’ and is shown in Fig. 1. The essential 
components of the apparatus were a mercury displacement 
pump, a gas burette, a sample holder (shown in detail), and 
a manifold system wherein the gas pressure could be varied 
from a high vacuum to the pressure in a commercial nitrogen 
cylinder. 

The gas burette was an addition to the original apparatus, 
and was added so that the porosity, as well as the capillary 
pressure curve, could be determined with the same apparatus. 
Prior to this modification only one of the two quantities needed 
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FIG. 6 — CAPILLARY PRESSURES BY STATIC AND DYNAMIC METHODS 
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FIG. 7 — J-CURVE FOR ALL CORES, EDWARDS FORMATION, JOURDAN 
TON FIELD 


to determine porosity could be measured with the mercury 


apparatus, ¢.e., bulk volume. This made it necessary to deter- 
mine the solid volume in a separate apparatus. The addition 
of the gas burette, connected to the chamber of the mercury 
apparatus, made it possible to determine the solid volume of 
the sample by the conventional gas-expansion technique. A 
steel blank of known volume was used to calibrate the system. 
thus simplifying the calculation of the solid volume of the 
core sample. This change obviously increases the utility of 
the apparatus, and should make possible the design of a more 
compact and portable unit for field use. 

In operation, the core sample was enclosed in the pressure- 
tight chamber, the system was pressurized to a predetermined 
pressure, and the solid volume of the sample was determined 
by the conventional gas-expansion technique. Next, the system 
was evacuated, and the volume between two reference marks 
on the Lucite below the chamber was 
determined by introducing mercury into the system with a 
mercury pump. This volume was subtracted from the known 
volume of the sample chamber to obtain the bulk volume of 
the sample. The vacuum pump was then isolated from the 
system and gas was admitted in increments, thereby increasing 
the pressure on the mercury surrounding the sample. The 
penetration of mercury into the pores of the core was indi- 
cated by a recession of the mercury-gas interface from the ref- 
erence mark on the upper Lucite window, and the amount of 


windows above and 


penetration was determined by advancing the piston of the 
mercury pump until the mercury meniscus was again at the 
reference mark. The procedure of alternately building up the 
pressure and advancing the pump piston was continued until 


the desired maximum pressure was reached. 


DYNAMIC CAPILLARY PRESSURE 
MEASUREMENTS 


During the course of research in fluid mechanics conducted 
to determine relative permeabilities of reservoir rock to oil 
and gas, an apparatus designed in principle after that devel- 
oped by Hassler’ was constructed which makes it possible to 
fix directly the capillary pressure between gas and oil, and 
to maintain that at a fixed and known level while 
the fluids are flowing 


pressure 
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FIG. 8—J-CURVE FOR LIMESTONE CORES, EDWARDS FORMATION 
JOURDANTON FIELD 


The essential features of the Hassler-type apparatus. used 
to measure dynamic capillary pressures, are shown in Fig. 2 
Each of the two end plates consists of a central dise dia 
phragm and an outer ring diaphragm. The two diaphragms 
of porous porcelain, are isolated from each other by means 
of a metal cylinder. The face of each end plate is grooved to 
facilitate the passage of gas. The oil is flowed through the 
outer rings, and the center discs are utilized for measuring 
the pressure in the oil phase. 

In operation, gas and oil were admitted to the core at two 
different pressures, this pressure difference being the capillary 
pressure. Since it was possible to measure the pressure drop 
in the oil phase from one core face to the other by the use of 
the center section of the porous plate, the gas pressure drop 
was adjusted at the outlet end until its pressure drop just 
equaled the pressure drop in the oil phase. When this condi 
tion had been attained, equal capillary pressures existed at 
each end of the core and the saturation was essentially the 
same throughout the core, if the sample was homogeneous 
Any saturation variations due to expansion of the gas were 
very small because of the low pressure differences employed 
While approaching this state of equilibrium, the pressure drop 
in the gas phase was maintained at a higher value than that 
in the oil phase. In consequence, prior to equilibrium th¢ 
saturation at the inlet end was lower than that at the outlet 
end because of the decreasing capillary pressure gradient in 
the direction of flow. The liquid in each portion of the cor 
was, therefore, being displaced to a region of higher satura 
tion. The capillary pressure data obtained under these condi 
tions should correspond, as nearly as is physically possible 
to the drainage curve arrived at by static methods. 

After attaining equilibrium in the manner described above 
the core sample was removed and the saturation was deter 
mined by weighing. 


Determination of Lithologic Properties 
of Core Samples 


A petrographic study was conducted to determine the lith 
ologic character of the core samples which were used in corre 
lating capillary pressure curves. These samples were all from 
the Edwards formation (Cretaceous Age) of the Jourdanton 
Field in southwest Texas. 
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lhese samples were all carbonate rock, but some were 
classified as limestone and the others as dolomite or dolomitic 


limestone. 
These samples were predominantly of the intercrystalline 
porosity type, though some samples had _ solution-enlarged 


pores as well. 

The samples were divided by visual methods into the follow- 
ing textural groups, classified according to grain size: 

1. Microgranular — characterized by not being visibly crys- 
talline in reflected light under 20-power magnification, but 
individual grains being barely discernible in thin section with 
50-power magnification. 

2. Fine-grained — visibly crystalline under 20-power mag- 
nification, but individual grains not recognizable. Average 
grain size about 0.02 mm. 

3. Coarse-grained —- individual grains easily recognizable 
under 20-power magnification. Dimensions of individual grains 


exceed 0.04 mm 


RESULTS 


Comparison of the Restored-State and 
Mercury-Injection Methods 


Capillary pressure-liquid saturation curves obtained by both 
the displacement cell and the mercury-injection apparatus 
were found to be nearly identical when plotted with appro- 
priate scale relationships between the water-gas capillary pres- 
sure and the mercury capillary pressure. This identity of 
results indicates that the capillary pressure-saturation curves 
obtained by the two methods are geometrically similar, one 
curve being obtainable from the other for any given core by 
multiplication of the capillary pressure by a constant. In Figs. 
} and 4 are presented capillary pressure curves obtained by 
the two methods of capillary pressure measurement for sand- 
stone and limestone cores, respectively. 

The results of these studies confirm in general the results 
published by Purcell,’ with the exception that Purcell’s use 
of a constant conversion factor for all cores is questionable 
in the light of the data shown in Figs. 3 and 4. The constant 
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FIG. 10—J-CURVE FOR MICROGRANULAR LIMESTONE CORES, ED 
WARDS FORMATION, JOURDANTON FIELD 


conversion factor used by Purcell was arrived at by assuming 
the following values: 

Surface tension of water: 70 dynes per cm 

Surface tension of mercury: 480 dynes per cm 

Contact angle of water against the solid: 0 

Contact angle of mercury against the solid: 140 
The conversion factor is then, 


(Mercury Capillary Pressure ) (480) (cos 140°) 


(70) (cos 0°) 


_ 


<0 


(Water-Air Capillary Pressure) 
Purcell, in presenting his results. used as a first approxima 
tion for the conversion factor a value of 5.0. and found fair 


agreement between the results of the two methods. 

In contrast, it was found in this study that the agreement 
could be made almost exact if a different factor. ranging from 
5.4 to 8.3, were used for each core sample. The fact that the 
lowest factor encountered experimentally was 5.4 suggests that 
the theoretical value of 5.25 may be a limiting value. This 
supposition is substantiated by the known sensitivity of water 
surface tension to contamination, making it quite possible that 
the figure of 70 dynes per cm is a maximum. Therefore, the 
use of a higher constant factor than that used by Purcell seems 
justified, and is supported by experimental experience. Data 
obtained to date indicate that the value of such a factor should 
be set tentatively as 6.4 for limestone and 7.2 for sandstone 
It is expected that future experience will dictate revisions, but 
for the present these proposed values will 
better correlations. 

As stated, it was found that agreement 
methods for measuring static capillary pressure could be made 
almost exact, provided the appropriate conversion factors were 
employed. The determination of one capillary pressure point 
hy the restored-state technique permits the calculation of thi- 
factor for a core, if the point is properly located. Fig. 3 may 
be used to exemplify what is meant by proper location of the 
point. It may be observed in this figure that the wetting phase 
saturation is most sensitive to changes in capillary pressure 
along the horizontal portion of the curve. It is here that the 


in general. give 


between the two 


correlating factor must cause the two curves practically to 
coincide, if large errors in saturation are to be avoided. The 
known point, therefore, must be in this region. The determina- 
tion by this procedure of the correct correlating factor makes 
it possible to obtain a complete curve in detail by the rapid 
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mercury-injection method, and then to convert the mercury 


curve to a restored-state curve 


Comparison of Dynamic and Static Capillary 
Pressure Measurements 


Comparisons are presented in Figs. 5 and 6 between static 
and dynamic capillary pressure data for a sandstone and a 
limestone sample. The two methods show almost exact agree- 
ment for each core, and confirm the validity of the assumption 
that static capillary pressure data apply to flowing systems 
as well as to static systems. and hence can be utilized in the 
solutions of dynamic problems of fluid flow. It must be remem- 
bered, however, that in each of the two methods only two 
phases, oil and gas, were present. Neither the two-phase sys- 
tem of oil and water nor the three-phase system of water-oil- 


gas was investigated 


Correlation of Capillary Properties 
of Rock with Rock Texture 


Leverett’s’’ capillary pressure function relationship, pro- 
posed by Rose and Bruce" for correlating capillary pressure 
curves for samples from the same formation, is 

T(S.) = (P./y cos 0) (K/¢)°" 
where 
P. = capillary pressure 
K permeability 
porosity 
interfacial tension between the two fluids, and 
the contact angle between the interface separating the 
two fluids and the surface of the rock. 
y cos 0) (K/@)°" is dimensionless, and, as 


The group (P 
the S, referring to the satu- 


shown, is designated as J(S,). 
ration of the wetting fluid. A plot of J(S.) versus wetting 
saturation results in the capillary-pressure-function 


phase 
Chis curve is more familiarly known as the J-curve. 


curve. 

The J-curve was employed to correlate capillary pressure 
liquid saturation data obtained from core samples of the Ed- 
wards formation in the Jourdanton Field. The capillary pres- 
sure-liquid saturation data were obtained by the restored-state 
method on cores from various depths in seven wells of this 
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field. The porosities of these cores varied from 9.5 to 28.2 per 
cent, and the permeabilities varied from 0.346 to 186 md 

On plotting the capillary-pressure function, J(S,), against 
wetting phase saturation it was found that all the data points 
fell roughly along the same trend. The result. presented in 
Fig. 7, gives support to the proposal of Rose and Bruce.’ It 
was felt, however, that some refinement of the curve would 
be desirable. 

Subsequent studies of the data showed that the curve might 
be better for cores in the same lithologic classification. Conse 
quently, the cores were divided into groups and plotted as 
separate J-curves according to whether they were limesten¢ 
Fig. 8, or dolomite, Fig. 9. 

The only improvement noted by this division was for the 
dolomite samples. As shown by a comparison of Fig. 9 with 


Comparison Between Measured 
Indica 


Measured 
Permeability 
(m 


186.0 
70.3 
137.0 
154.0 
1.03 
0.346 
20.6 
25 
1.97 
7.17 
34.0 
15.0 
1.99 
£.55 


»>90 


Porosity 
(fraction) 


.279 
225 
.254 
.246 
117 
104 
.206 
.159 
143 
127 
20] 
237 


134 


Well No. 1 


Well No. 2 


160.0 
113.0 
15.9 
2.85 
30.6 
10.7 
53.5 
8.9 
1.07 
2.79 
14.3 
93.0 
28.5 
2.62 
5.07 
1.26 
6.13 
10.8 
67.5 
84.5 


2.5 


Well No. 


Well No. 4 


Well No. 5 


) 


53.0 
13.8 
39.5 
74.0 
Average Deviatic 
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Fig. 7, the scattering of points was reduced materially, giving 
much less spread of the data points. 

To investigate quantitatively this improvement, the perme- 
abilities of 40 dolomite samples were calculated using the J- 
curve for all cores, Fig. 7, and the J-curve for dolomites only; 
Fig. 9. These calculated permeabilities are compared with the 
measured permeabilities in Table I. The per cent error, or 
deviation from the measured value, was determined for each 
core. The average error for all samples was found to be 71.5 
per cent when using the J-curve for all cores, Fig. 7. The use 
of the dolomite J-curve, Fig. 9, reduced the average deviation 
to 48.0 per cent. The variations in porosity and permeability 
for this group of cores covered such a wide range, threefold 
for porosity and 500-fold for permeability, that these average 
per cent errors should be fairly representative for any similar 


and Caleulated Permeabilities for Dolomite Cores Using J-Curves 


ted 


/-Curve for All Edwards Samples /-Curve for Dolomite Only 


Calculated 
Permeability 
K, (md.) 


Calculated 
Permeability 


K fans Deviation 


Deviation 
% 


22.6 
53.6 
2.9 
77.9 
22.3 
386.0 
1.0 
31.4 
61.9 
9.5 
11.8 
15.8 
24.6 
18.7 
61.5 
50.7 
29.4 
54.7 
69.1 
56.5 
76.6 
7.3 
4.5 
27.8 
42.3 
25.2 
11.9 
272.0 
115.0 
16.5 
55.5 
29.0 
16.5 
8.0 
0.3 
21.0 
18.1 
12.4 
39.7 
0.1 
48.0 


144.0 
108.0 
141.0 
274.0 
1.26 
1.68 
20.8 
3.30 
3.19 
7.92 
30.0 
52.1 
1.50 


10.2 
121.0 
38.0 
110.3 
72.8 
599.0 
33.0 
65.0 
126.4 
54.7 
23.6 
62.0 
6.0 
66.1 
126.7 
30.6 
0.9 


7 


Die 
140.3 
13.1 
137.0 
14.8 
33.7 
5.2 
19.4 
9.8 
24.7 
105.0 
200.0 
25.0 
37.7 
0.7 
33.6 
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189.0 
324.0 
1.78 
2.42 


97 


Zi 


4.14 
1.16 
11.1 
12.0 
72.9 
2.11 
7.56 
5.19 
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225 
15.7 
116.0 
144.0 
7.79 
3.80 
0.785 
6.09 


7 
al 


81.2 
39.7 
50.5 
82.5 
50.0 
x15 


97.7 


) 100 


ACTIONS, AIME 


m. “¢ 
K.-K.,, 
K,, 





T.P. 3025 


group of cores from this formation. In this case, therefore, 
the breakdown of Edwards core samples into two groups, lime- 
stone and dolomite, improved the correlation for the dolomite 
group. 


The limestone samples of Fig. 8, were subdivided into two 
groups according to their texture. One group was the micro- 
granular type, and the other was coarse-grained. As shown in 
Figs. 10 and 11, respectively, this breakdown into lithologic 
groups according to texture appears to improve the correla 


tion for both groups. 

The dolomites divided 
grained dolomites as one group, and coarse-grained dolomites 
as the other. The improvement in the individual J-curves over 
the J-curve for all dolomites was negligible. 


were into fine and fine-to-coarse 


The two remaining groups, fine-grained and _fine-to-coar=« 
grained limestone were neglected, since the number of cores 
available was not judged sufficient. 

It is apparent from these data that the J-curve may have 
utility in correlating capillary pressure data obtained on sam 
ples from a single formation, but that the results derived from 
its use may be of limited precision. It is also apparent that 
the correlation may be improved by restricting it to specific 
lithologic types of samples 


CONCLUSIONS 


1. A comparison of the mercury-injection method with the 
restored-state method of measuring capillary pressure showed 
that the two yielded nearly identical results if the proper con 
version factor were chosen for each core sample, but that the 
conversion factor was not constant for all samples. This factor 
was found to average 6.4 for limestone samples and 7.2 for 
sandstone samples. 

2. Static and dynamic methods of obtaining drainage capil 
lary pressure curves gave nearly identical results. This identity 
confirms the validity of the use of static capillary pressure 
data in dynamic problems of fluid flow. 

3. Data on a number of samples of Edwards limestone sup 
port the use of the J-curve in correlating capillary pressure 
data obtained on core samples from a single formation. The 
correlation may be improved by restricting its use to specifi 
lithologic types from the same formation. 
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By W. R. Purce Houston, Tex., Junior Mem 
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Shell Oil Co... 


Of the various interesting experimental results reported in 
Brown’s paper, those concerning the conversion of water /ni- 
trogen capillary pressures appear to warrant additional com- 
factors. 
ranging from 5.4 to 8.3, are higher than the factor of 5.25. 


ment. Brown suggests that his observed conversion 
previously proposed by the writer. because of a reduction by 
contamination of the surface tension of water to a value below 
the 70 dynes per cm which was used in calculating the factor 
would however, that the mercury would 
be as susceptible to contamination as the water, if not more 
so, and it seems difficult. therefore. to resolve all of the dif- 
ferences merely on the basis of surface tension lowering alone. 
It should be pointed out that the factor 5.25 as originally pro- 
posed was calculated on the assumption that the pores of a 
rock sample behave as a group of circular capillary tubes. It 
would seem reasonable therefore, to conclude that the depar- 
from the idealized 


5.25. One suspect, 


ture ot the observed conversion factors 
theoretical value of 5.25 reflects to a large degree the effect 
of departure of pore geometry from the hypothetical cylindri- 
cal tube. It is my opinion that. in the consideration of capil- 
lary pressure phenomena, the effect of pore geometry is equally 
is great as that of surface and interfacial forces,’ and it is 
this effect of pore geometry which hes been neglected by 
Brown, and previously by the writer, in considering conversion 
factors for changing capillary pressures as determined by one 
system to equivalent values for another system 
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DISCUSSION 


By {. Bruce, Carter Oil Co. Tulsa, Okla., Member AIME 


It is of interest to me to note that Brown suggests refine- 
ment of the J-function technique to the extent of subdivi- 
porous material from a formation into two or more 
distinct groups. This seems to be a desirable practice if the 
rock is divisible into distinctly different classifications such as 
dolomite or calcite as in Figs. 8 and 9 above. However. 
petrographic classification is not the key to the 
J-function, probably because a number of petrographic fac- 
grain size and standard deviation are already 
accounted for macroscopically by the permeability 
that go to make up the J-function. * * * 
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VARIABLE CHARACTERISTICS OF 
THE OIL IN THE TENSLEEP SANDSTONE RESERVOIR, 
ELK BASIN FIELD, WYOMING AND MONTANA 


RALPH H. ESPACH, MEMBER AIME, AND JOSEPH FRY, U. S. BUREAU OF MINES, LARAMIE, WYO. 


ABSTRACT 


In the spring of 1943, when it was evident that the Tensleep 
sandstone in the Elk Basin Field, Wyoming and Montana 
held a large reserve of petroleum, Bureau of Mines engineers 
obtained samples of oil from the bottom of nine wells and 
analyzed them for such physical characteristics as the volumes 
of gas in solution, saturation pressures or bubble 
shrinkage in volume caused by the release of gas from solu 
tion, expansion of the oil with decrease in pressure, and othe: 
related properties. The composition of the gas in solution in 
the oil was studied. The pressures and temperatures existing 
in the reservoir and the productivity characteristics of the 


points 


oil wells were determined. 

The data obtained indicate that the oil in the 
Reservoir of the Elk Basin Field has unusually varying physi 
cal characteristics, such as a saturation pressure of 1,250 psia 
and 490 cu ft of gas in solution in a barrel of oil at the crest 
of the structure and a saturation pressure of 530 psia and 134 
cu ft of gas in solution in a barrel of oil low on the flanks 
The hydrogen sulfide content of the gas in solution in the oil 
varies from 18 per cent for oil on the crest to 5 per cent for 
oil low on the flanks of the structure. Of even greater signifi 
cance is the fact that these and other variable characteristics 
of the reservoir oil are related to the position of the oil in 
the structure. Many geologists and petroleum engineers hav 
considered that all the oil in a petroleum reservoir has rathe: 
uniform physical characteristics and that equilibrium condi 
tions prevailed in all underground accumulations of oil and 
gas; that such is not alwavs so is berne out by the results of 
the study by the writers 


lenslee P 


INTRODUCTION 


The Rocky Mountain region is one in which may be found 
striking examples of the unusual! in oil and gas accumulations 
as is evident from the following: The high helium content 
(7.6 per cent) of the gas in the Ouray-Leadville limestone 
sequence in the Rattlesnake Field, New Mexico, and gases of 
similar helium content in other fields; 50° to 55° APTI gravity 
distillate in solution in carbon dioxide gas and recoverable 
through retrograde condensation, in the North McCallum 
Field, Colorado; the occurrence of gas. oil, or both in closely 
related structures contrary to the usual concepts of gravi 
metric segregation; the accumulation of gas and/or oil in 
structures closely related to other structures that apparently 
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are more favorable but do not contain oil or gas accumula- 
tions; the high hydrogen sulfide content (as high as 42 per 
cent) of the gas associated with oil in some fields in the Big 
Horn Basin, Wyoming; and the wide range of fluid charac 
teristics found in the Elk Basin reservoir. 

Elk Basin, an interesting old oil field that has been pro- 
ducing oil from the Frontier formation since 1915, is situated 
in a highly eroded basin resulting from the erosion of the 
crest of an anticline and some of the underlying softer shales. 
The field came back into national prominence during 1943 
when it became known that it was the largest single reserve 
of new oil discovered in the United States that year. The 
Tensleep sandstone was found to contain oil in November. 
1942, when a well drilled to a depth of 4,538 ft (44 ft into 
the Tensleep sandstone) flowed oil at the rate of 2,500 B/D. 
By the end of 1949, 137 oil-producing wells and five dr’ holes 
had been drilled, and approximately 32 million bbl of oil had 
been produced. Approximately 6,000 acres may be consid- 
ered productive of oil in the Tensleep Reservoir, and estimates 
of the oil that will be produced average 200 million bbl. 

The Tensleep Reservoir has further interest because it ha- 
greater closure than any oil field in the Rocky Mountain 
region; the closure of the Elk Basin anticline is variously 
estimated at 5,000 to 10,000 ft. of which the top 2,000 ft of 


the structure contained oil 


SUBSURFACE OIL SAMPLING 


Fig. | is a structural map of the Elk Basin Tensleep Reser- 
voir, on which the nine wells used in this study and the num- 
bers corresponding to the well designations hereafter referred 
to are shown. Wells 1, 2, 3, 4, and 8 were tested and sampled 
during October and November. 1943, and Wells 5, 6, 7, and 
9 during June and July, 1944. 

An electromagnetic type sampler developed by the Bureau 
of Mines and described by Grandone and Cook’ was used in 
obtaining the subsurface oil samples. 

(s the wells were tubed nearly to bottom, the sampler wa- 
run as far as possible in the tubing but never below the top 
perforations. 

The following procedure was used in testing and sampling 
the wells: A well was shut in for at least three days, after 
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which s‘iut-in subsurface pressures and temperatures were 
measured. Generally, the well then was allowed to flow at a 10 
to 20-bbl-per-hour rate for approximately one day, when a 
flowing temperature and pressure traverse was run. One day 
later, at the same production rate, a subsurface sample was 
obtained and transferred to a storage container. A second sub- 
surface sample of oil then was procured; and, if in transfer- 
ring, the pressure and volume data duplicated those of the first 
sample, both were considered satisfactory for analyses. After 
the sampling was completed the flow rate of the well was 
increased to two higher rates and subsurface pressure tra- 
verses were obtained at each rate. A 12-hour subsurface pres- 
sure build-up was obtained when the well was shut in follow- 
ing the highest rate of flow. A pressure gauge recorded the 
tubing and casing pressures at the wellhead during the tests 
The oil produced by the well during the tests was gauged in 
lease tanks at regular intervals and the separator gas was 
measured with a recording orifice meter. The essential data 
pertaining to the wells are given in Table 1. 


Wells 7 and 9 were not producing when sampled. Well 7 
had flowed most of the indicated production; but the flow 
gradually lessened, and the well finally stopped flowing. It 
was sampled several days later. Well 9 was swabbed unsuc- 
cessfully to cause it to flow: it had been standing for about 
30 days before sampling. and as tubing had not been run to 
bottom, it was necessary to sample the oil at a relatively high 
point in the well. 


It will be noted that Wells 2 and 3 have essentially the 
same structural position, and for the purposes of this study 
it might seem that one well would suffice. Well 3, however 
produced from the upper 44 ft of the Tensleep formation, and 
Well 2 was completed so that only the zone between 66 and 
192 ft below the top of the Tensleep formation produced oil 
through the tubing. Sampling both wells provided an excellert 
opportunity to study the possibility of the existence of more 
than one reservoir in the Tensleep formation, a postulate 
that had merit because of the rather extensive dolomitic zones 
extending throughout the sandstone body at depths greater 
than 70 to 80 ft in the formation. The almost identical results 
obtained in the analyses of the oil samples from both wells 
indicated that the oils were the same and that the Ten leep 
formation contained only one reservoir 


RESERVOIR-OIL CHARACTERISTICS 


The subsurface oil samples were analyzed by the differential 
gas-liberation method outlined by Grandone and Cook.’ The 
laboratory analyses of the oil samples were made at the Petro 
leum Experiment Station (now Petroleum and Oil-Shale Ex 
periment Station), Bureau of Mines, Laramie, Wyo., during 
several months following sampling of the wells. The general 
results are given in Table II. Individual gas-liberation data 
for the nine samples are not presented, but Figs. 2, 3, and 4 
show how the data changes for the oil at wells 1, 6, and 8 


In these figures, the graphs of “expansion of sample” ar 
self-explanatory. The values for the gas in the “gas in solu 
tion” graphs are figured at conditions of 60°F and 14.4 psia 
pressure. The “relative oil-volume” graphs show the volume 
relationship between the oil in the reservoir under reservoii 
conditions and the residual oil at 60°F. The “density of sam- 
ple” graphs indicate the change in density as changes in pres- 
sure and release of solution gas occur. The “specific gravity 
of gas in solution” graphs show the gravity of the gas that i- 
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in solution in the oil at the different pressures and in volume 
as shown by the “gas in solution” graphs. The values for the 
specific gravity of the gas in solution were obtained by using 
1 220-ml gas-density balloon, weighing samples of gas liber- 
ated between pressure points and calculating the composite 
to any pressure by using the volumes and weights of the sev- 
eral increments of liberated gas. 


The “hydrogen sulfide content of gas in solution” graphs 
show the percentage of hydrogen sulfide in the gas that is in 
solution in the oil at the different pressures and in volume as 
shown by the “gas in solution” graphs. The percentages of 
H.S were obtained by running both Tutwiler and Orsat anal- 
yses of the increments of liberated gas during analyses of 
samples from wells 5, 6. 7, and 9 and Orsat analyses for acid 
gases of the increments of liberated gas during analyses of 
samples from Wells 1, 2, 3, 4, and 8. Where Orsat values 
alone were determined. 5 to 61 per cent was subtracted from 
them as a correction for the carbon dioxide content of the 
gas. The “hydrogen sulfide content of gas in solution” graphs 
represent values obtained from: The determined percentages 
of hydrogen sulfide in the several increments of gas liberated 
between any two pressures; the volumes of the several incre 
ments; and the calculated composites of these to various 
pressures. 


Of interest in connection with the “gas in solution” graph- 
are the data on the H.S contents of the several increments of 
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FIG. 1— ELK BASIN, 





Date of 
well 
completion 


9. 9-43 
3- 9-43 
4-18-43 


Date of 
sampling 


11- 4-43 
10-23-43 
10-15-43 
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Table | 


Oil produced 
to sampling 
date, 


20.000 
114,000 
130,000 


Reservoir 


r, Elk Basin Field 


Top 


Gravity of 
produced 
oil, “API 


31.8 
31.4 
31.1 


Data Relative to Wells Sampled, Tensleep 


of Tensleep formation 


Original 
pressure 
psia’ 


1853 
1917 
1934 


Temper- 
ature, 
112 
115 
116 


Sampling point 


Eleva- 

tion, 
ft 
488 
283 

+413 


Flow rate 


Pressure, while sampling 
i B/D 


psia 
1564 
1734 
1644 


922 
454 
536 


31.0 
30.6 
30.0 


10- 4-43 
6-11-44 
6-21-44 
6-15-44 
10- 9-43 
5-31-44 


34,800 
6,200 
9,200 
13,800 29.1 
52,700 28.4 

200 27.1 


7- 6-43 
1-28-44 


‘Calculated from estimated original reservoir pressure of 2,080 psia a 
From temperature gradient chart for the Elk Basin structure 


Production by swabbing in an attempt to cause the well to flow. Wel 


gas liberated during the differential gas liberation of a sub 
surface sample of oil. Such data from the sample of oil from 
Well 1 are given in Table III. The great change in H.S con 
tent of the gas as it is released from solution when the pres 
sure is lowered is significant. It reflects the solubility of H.S 
in crude oil. The “hydrogen sulfide content of gas in solution 
graphs indicate that at Elk Basin the H.S content of the gas 
from the oil in the reservoir would increase as reservoir pres 
sures decreased in accordance with the graphs provided the 
gas/oil ratio during production would not exceed that in solu 
tion in the oil. However, gas/oil ratios will increase in solution 
gas-drive reservoirs as the pressure declines. The percentag: 
of hydrogen sulfide in the gas evolving from solution also will! 
increase as reservoir pressure declines; this relationship i- 
shown in Table III for an oil sample from Well 1: The rela 
tionship is similar for oil samples from other wells. Thes« 
two factors will cause a decrease in the percentage of hydro 
gen sulfide in the produced gas as reservoir pressures decline 
for the volume of gas produced in excess of the volume of 
gas in solution in the oil will have a smaller percentage con 
tent of hydrogen sulfide and will act as a diluent for the 
solution gas. Later in the life of the reservoir when declining 
gas/oil ratios become effective. an increase in the H.S con 
tent of the produced gas can be expected. 


Owing to the method by which the data were obtained 


those data relating to the specific gravity of the increments 


Table LI 


sea-level 


was pumped after 


General Results of 


420 
486 
240 
shut-in 
440 


shut-in 


1905 
1570 
1720 
1882 
2023 
1770 


12) 30 
125 265 
129 557 
134 773 
137 1044 
14] 643 


2049 
2120 
2218 
2319 
2392 


2472 


683 
887 
1107 


1 datum and density of subsurface samples 


sampling 


of gas liberated from the subsurface oil samples and the 
hydrogen sulfide content of the gas samples are the least 
accurate of all results reported. In fact, the data show changes 
in one or two instances that are inconsistent with the data 
from other samples. 


The writers are indebted to two of the operators for in- 
formation on the viscosity of the reservoir oil from Well 2. 
Oil from Well 2 is representative of the oil in the upper 500 
or more feet of closure. The viscosity in centipoises of the oil 
at a temperature of 123°F and at several gauge pressures 
and the estimated viscosities of the oil near the edge of the 


reservoir are given in Table I\ 


DISCUSSION 


\ review of the data given in Table IJ and shown in Figs 
2, 3, and 4, reveals a very interesting conclusion that should 
be considered in all studies of petroleum reservoirs; namely, 
that the oil in a reservoir may have rather widely varying 
physical characteristics 

In general, many engineers and geologists have reasoned 
that, considering the length of time involved in the formation 
of an oil field, equilibrium must have been attained between 


Analyses of Subsurface 


Gravity, 


At 
original 
reservoir 
pressure 


45.4 
14.9 
15.4 
14.3 
14.6 
412.3 
39.2 
37.3 


y 34.9 


API 
Residual 
oil, after 
liberation 
of gas 
30.7 
30.2 
30.6 
30.6 
30.2 
29.5 
29.) 


27.7 


27.0 


Oil Samples, Tensleep Reservoir, Elk Basin Field 


Gas- 
liberation 
tempera- 

ture, °F 


121 
121 
121 
123 
129 
133 
134 
140 
141 


Saturatior 
pressure 
psia 


1250 
1255 
1250 
1235 
1177 

990 


99 


772 
695 
530 


(ras In 
olutior 
u ft/bbi 


esidual oil 
190 
178 
169 
159 
137 
$50 
225 
205 


134 


Based upon differential gas-liberation method. Base, 14.4 psia and 60°! 
*Ratio of reservoir-oil volume (at original reservoir conditions at top of Tensleep formation) to 
Volumes per volume per psi pressure change, reservoir oil at pressures above saturation pressure 
'This figure probably is in error. Corresponding data from other samples indicate that it should have been approximately 75 x 10 


Relative 
oil 
olume 

292 
1.220 
1.221 
l A | » 
1.210 
1.178 
1.122 
1.116 
1.086 


68 x 10 
64x 10 
63x 10 


residual-oil volume at 60°F 


Specific 
gravity 


of gas in 


solution 
0.989 
996 
1.001 
1.011 
1.028 
1.048 
1.021 
1.026 
1.041 


Hydrogen su! 

fide content of 

gas in solution 
per cent 
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the fluids in a reservoir and if any of the oil was saturated 
with gas under existing pressures and temperatures, then all 
of it was saturated or very nearly so, or, conversely, if some 
of the oil was undersaturated, then all of it would be under 
saturated to about the same degree. Some recorded observa- 
tions in other fields seem to indicate that equilibrium condi 
tions may not exist in all fields, but incontrovertible proof of 
that has been lacking. Estabrook and Rader’ noted the pres- 
ence of a “zone of dead oil” in the Second Wall Creek sand 
in the Salt Creek Field and reported that “the distribution 
of dissolved gas in the oil seems to have been fairly uniform 
over the field except toward the edges.” 


Study of the magnitude of change of properties of the sul 
surface oil samples from Elk Basin definitely indicates that 
the fluids in a reservoir are not necessarily in equilibrium. It 
is concluded that either extremely long periods of time are 
necessary to effect equilibrium or that accumulation can be 
quite young, or is still taking place in some petroleum reser 


voirs 


RELATIONSHIP OF OIL CHARACTERISTICS 
TO RESERVOIR STRUCTURE 


The relationship of the characteristics of the oil to the loca 
tion of the oil in the Elk Basin, Tensleep Reservoir, is shown 
in Fig. 5. In Fig. 5 values of the saturation pressure, gas in 
solution, relative oil volume, density of the reservoir oil, spe 
cific gravity of the gas in solution, and hydrogen sulfide con 
tent of the gas in solution for the nine subsurface oil samples 
are plotted as ordinates against the common abscissa of sea 
level elevations of the top of the Ten-leep formation over the 
structure. The nine wells shown in Fig. 1 are identified by 
number along the abscissa of Fig. 5 at positions corresponding 
to the elevation of the top of the Tensleep in each well 
Although the wells are plotted according to formation “tops 
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FIG. 2—RESULTS OF DIFFERENTIAL GAS-LIBERATION ANALYSIS AT 
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Table Lil — Hydrogen Sulfide Content of Liberated 
Gas, Differential Gas Liberation, Oil Sample from 
Well 1 


Pressure range in 
which gas increment 
liberated, psia 


Hydrogen sulfide 
content of gas 
increment, per cent 
1250 to 835 7.6 
835 to 725 5.1 
725 to 625 5.9 
625 to 535 6.8 
335 to 450 6.1 
450 to 365 10.2 
565 to 300 13.0 
300 to 225 16.0 

225 to 165 
165 to 120 26.7 
120to 85 
85 to 60 
H0to 45 
to 11 


h wells came from greater 
depth, depending on the penetration (160 to 200 ft) of the 
formation in each well 


e oil from each of these some 


The composite sample data for the oil from each of the 
wells sampled, as derived from the individual analyses, are 
shown in Fig. 5 to have a very definite relationship to the 
from which the oil came. To the 


writers’ knowledge. this is the first time that results of detailed 


structure of the reservoir 
reservoir oil sampling have been published showing such a 
relationship to exist. 
\ study of Fig. 
surface sample of oil from Well 7 was not truly representa- 
oil at Well 7 and that the sample had 
been conditioned somewhat. This was referred to in the sec- 
The data for the specific 


5 would indicate that possibly the sub 


tive of the reservoir 


tion on subsurface oil sampling. 
gravity of the gas in solution show the least consistent rela- 
all the data. 


this report to the manner in which these data were acquired 


tionship of Attention has been called earlier in 
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Table IV — Viscosity of Reservoir Oil 
Oil from edge 


wells, centipoises 


Pressure, Oil from Well 2, 
i (estimated) 


psi centipoises 

2500 1.59 1.6 
2000 1.50 
1500 1.42 
1150 1.35 
1000 1.40 
800 1.50 
600 1.65 
100 1.82 
200 2.05 
0 4.50 


The physical characteristics of the oil throughout the crestal 
500 to 600 ft of the reservoir are very much alike — all oil 
samples from the crestal part of the reservoir showed a satu 
ration pressure of 1,235 to 1,250 psia and contained 460 to 
490 cu ft of gas in solution based on 1 bbl of residual oil. As 
a consequence of this gas in solution, 1.22 bbl of oil in the 
reservoir are required to produce 1 bbl of residual, or approxi 
mately stock-tank, oil. At a depth of about 100 ft above sea 
level, a decided change begins to occur in the character of 
the reservoir oil—the saturation pressure becomes lower and 
the volume of gas in solution decreases. The other characte: 
istics begin to change also. These changes continue with depth 
until the oil from a well in which the top of the producing 
formation is at an elevation of 1,107 ft below sea level. con 
tains only 134 cu ft of gas in solution per barrel of oil (on 
a residual oil basis) and its saturation pressure is only 530 
psia. Edge water is present in the reservoir at an elevation of 
about 1,500 ft below sea level. The writers have no informa 
tion on the characteristics of the oil coming from wells in 
which the top of the Tensleep sandstone is below a depth of 
1,107 ft below sea level. Characteristics may be inferred from 
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extension of the several curves to elevations of 1,200 or 1,500 
ft below sea level. 


The hydrogen sulfide content of the gas in solution in the 
oil as compared to the volume of gas in solution in the dif- 
ferent oil samples is a rather interesting relationship; the gas 
in solution in the oil on top of the structure contains 18 per 
cent hydrogen sulfide. that in solution in the oil on the edge 
of the structure contains 5 per cent, and the variation is almost 
directly proportional to the volume of gas in the different sam- 
ples. This suggests a plausible hypothesis for the movement 
and characteristics of the fluids in the reservoir which is con- 
sistent with the supposition that equilibrium has not yet been 
attained in the whole mass of reservoir oil and which has been 
included in another article’ concerning the Elk Basin Field. 
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DISCUSSION 
By Henry J. Welge. Carter Oil Co., Tulsa, Okla.. 
Vember AIME 
For some time geologists and engineers have tacitly assumed 
that the elapsed time during which petroleum reservoirs have 
been relatively undisturbed is sufficiently large so that their 


fluids are in complete equilibrium with one another. Hence 
the authors are to be commended for directing attention to 
the possibility of the occurrence of non-equilibrium conditions 
in certain petroleum reservoirs. Work recently performed at 


The Carter Oil Co. Research Laboratory substantiates this 
point. In the interest of adding to the picture. a portion of 


the work done recently is herewith presented. 


In order to test the possibility of an incomplete approach 
toward equilibrium by diffusion of gas downstructure in the 
Elk Basin Pool, calculations have been made of the minimum 
time required for the attainment of the present gas distribu 
tion at the downstructure wells (6, 7, 8, and 9 on Fig. 1) by 
diffusion. The calculations presume a constant dissolved gas 
content of 467 cu ft per bbl at all times at levels higher than 
the present sea level contour (Fig. 1). They also assume that 
no, or relatively little, gas was in solution elsewhere in the 
reservoir at the time of trap formation. In other words, the 
picture assumed is that of gas (or gas-rich oil) present ini- 
tially in the upper part of the reservoir, overlying relatively 
“dead” oil. It is well known that rather long intervals of time 
are required for a solute to diffuse through a liquid if long 
distances have to be traversed. The calculations were made in 
order to see how long a time interval is required in the in 
stance of the Elk Basin Pool. and to compare it with the geo 
logic time available. 

The relation governing unsteady state diffusion is mathe- 
matically identical with that treating heat diffusion. For con 
ductors of certain shapes. the solutions have been given 
graphically.’ 

The special case of a constant high concentration of gas at 
the head of a petroleum column, maintained by the presence 
of a gas cap diffusing downward into the oil, is analogous to 
heat flow through a slab of uniform cross-sectional area. from 
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a constant temperature source at one face toward an opposite 
insulated face 

The solution for the slab shape assumes a uniform cross 
sectional area normal to the direction of diffusion. This con- 
dition would prevail approximately on the elongated southwest 
and northeast flanks of the pool, as indicated in Fig. 1 of the 
paper under discussion. With the aid of this graphical solu- 
tion, it is possible to compute the time required for the present 
gas distribution to be attained. The value used in this caleu- 
lation for the diffusion constant, D, is 0.43 sq ft per year, and 
it compares favorably with values given in the literature.’ The 
lengths used for the diffusion paths toward the downstructure 
wells have been estimated with the aid of the field map, Fig. 1. 

Che computations are summarized in the table below 

1. Well No 6 7 8 

2. Gas-oil ratio, cu ft/bbl = G 350 225 205 
3. Gas concentration ratio, 

Y = (467-G) /467 

Horizontal distance (Fig. 1) 
contact to 


» 
0.56 
from oil-water 
well, ft a 
Horizontal distance from vil- 
contact to sea level 

b 1970 
0.692 


3440 1780 


water 
4970 
0.242 


7040 
0.430 


4970 
0.359 


contour 
Position ratio n ach 
Angle of dip (Fig. 1) 
degrees 20 12 20 20 
Inclined 
level contour to oil-water con- 

tact, ft = R 5280 
R’*, ft. x 10 27.9 
XY (from Y and n 0.33 
Diffusional time 


years x 10 


distance from sea 

7190 
51.9 
0.30 


5280 
27.9 


0.32 


5280 
27.9 
0.22 

RD 


21.4 36.2 20.8 14.3 


It will be noted from the last line in the table that the times 
required for diffusion are long, of the order of 15-40 million 
vears, and are not incompatible with the hypothesis that gas 
diffused downstructure from the crestal region of high gas 
concentration throughout the geologic time elapsed since trap 
formation, without reaching complete equilibrium in the time 
available. 


Several other possible explanations for the observed varia- 
tion in composition of the petroleum in the Elk Basin Field 
have also been considered. In all cases they can explain only 
relatively small variations in composition. For example. Sage 
and Lacey* have computed for a similar reservoir the variation 
expected because of gravitational segregation. They found 
that, over as much as 2,000 ft of closure. the expected varia- 
tion for methane was only about 1/20 of the average amount 
of methane in solution. This explanation obviously would not 
account for the differences observed at Elk Basin. 


Considering, finally, the possibility of segregation due to 
temperature difference between the top and bottom of the oil 
column, Sage and Lacey® have found that the change in solu- 
bility of natural gas expected with a temperature change of 
from 120°F to 140°F was also about 1/20 of the total gas 
-olubility. 

Apparently considerations other than the 
incomplete attainment of equilibrium lead to the conclusion 
that comparatively little variation in petroleum composition 
-hould be observed in a reservoir at equilibrium, even though 
its closure extends over several thousand feet. This conclusion 
is demonstrated by the existence of a number of oil reservoirs 
a large closure. but whose petroleum is of substan- 
considerable 


assumption of 


having 


tially uniform composition. Therefore. when a 
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variation in composition is encountered, as at Elk Basin, a 
satisfactory explanation may be incomplete diffusional ap 
proach toward equilibrium. 
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DISCUSSION 


By Michael J. Rzasa, Stanolind Oil and Gas Co., Tulsa, Okla 
Vember AIME 


Espach and Fry are to be commended on publishing mat 
rial which shows the change in physical characteristics of 
crude oil in a unit reservoir with elevation. However, there 
are several statements and conclusions in their paper which 
are worthy of attention. 


In Table ILI, the authors have shown the hydrogen sulfide 
content of liberated gas from oil samples from Well No. | 
These data show that during the first pressure increment of 
1.250 to 835 psia, the hydrogen sulfide content of the gas de 
creased from 7.6 to 5.1 per cent. Inspection of equilibrium 
vaporization constants for hydrogen sulfide would indicate 
that these results are to be expected in the particular pressure 
range under consideration. As further expected, the hydrogen 
sulfide content of further gas increments increases with pres 
sure reduction. | am not sure what the authors mean by stat- 
ing that the hydrogen sulfide content of the gas from the oil in 
the upper part of the reservoir would increase with pressure 
decrease, provided the gas-oil ratio during production did not 
exceed that in solution in the oil, and that increasing gas-oil 
ratios would result in decreasing the hydrogen sulfide content 
of the produced gas, because any volume of gas in excess of 
that in solution in the oil will contain smaller percentages of 
hydrogen sulfide than the solution gas. The amount of hydro 
gen sulfide which will come out of solution with pressure de 
crease is determined not only by the absolute amount of oil 
which releases this gas, but also by the vaporization charac- 
teristics of hydrogen sulfide. Since all the other components 
in the released gas will increase in absolute amount as the 
producing gas-oil ratio increases, it would seem to me that 
the hydrogen sulfide content would not necessarily decrease, 
but would appear in approximately the same proportion as 


given in Table ITI. 


One point of great significance is the statement of the 
authors that their work proves that equilibrium in the Elk 
Basin Sand has not been attained. This conclusion is based 
on the fact that the physical characteristics of the oil from 
the deeper wells differ considerably from those upstructure 
The change in elevation between Well No. 1 and Well No. 9 
is approximately 1,750 ft. The very fact that there is this 
change in physical characteristics of the oil with change in 
elevation shows that the conclusion reached by the authors 
is incorrect. Thermodynamic equilibrium demands that thi- 
change in properties occur with change in elevation. 


In ordinary engineering work it is customary to assume 
that pressure, temperature. and the composition of a phas« 
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are the independent variables sufhcient to define the state of 

the system. This conclusion ignores gravitational and other 

types of energy associated with the system and it assumes 

that the pressure is uniform on the system. Stated in other 

words, when gravity is taken into account, the phase rule is 
2—p, but becomes f = c + 3-p. 


no longer f = c + 2 


Under equilibrium conditions, the chemical potential or 
molal free energy of all components must be equal to zero. 
It can be shown from thermodynamics that relationship exists 
between the change in composition of a fluid and change in 
elevation under isothermal conditions. Previous publication 
of calculations on a Dominguez Field crude oil show that for 
a change in elevation of 2,000 ft, the bubble point of the oil 
at 160°F changes approximately 885 psi. Furthermore, the 
isothermal effects indicating change in composition with ele- 
vation would be expected to be enhanced by a gradual increase 
in temperature as depth is increased. 


With deeper depth, the crude oil will contain higher pro- 
portions of the heavier components and consequently a lower 
bubble point. With greater proportion of heavier components 
in a crude oil, the gas in solution and the relative oil volume 
will also decrease as depth from some datum is increased. 
Obviously, with increase in pressure the density of the oil will 
increase and the specific gravity of the gas in solution will 
also increase. These qualitative conclusions are borne out by 
the data published by Espach and Fry, and further substan- 
tiate the statement that thermodynamic equilibrium does exist 


at Elk Basin. 


It would be of significance and value if the chemical com- 
position of the crude oil at various elevations in Elk Basin 
were published as a part of this paper in order that such 
thermodynamic calculations can be made. 


AUTHOR'S REPLY TO MR. RZASA 


M. J. Rzasa’s comment that change in composition of reser- 
voir oil with elevation is to be expected at conditions of ther- 
modynamic equilibrium is correct. This does not, however. 
make incorrect the conclusions reached by the authors. The 
magnitude of change reported in the Elk Basin study between 
the properties of the reservoir oil from Wells 1 and 9, with 
some 1,750 ft difference in elevation, is many times too great 
to be consistent with changes expected at conditions of ther- 
modynamic equilibrium. Although changes in liquid compo- 
sition of the reservoir oil were not available to be compared 
with the change of composition as defined by thermodynamic 
equilibrium, it is apparent upon examination of Sage and 
Lacey’s' sample calculations, showing such small changes in 
composition, that equilibrium is not the explanation for the 
changes in properties of the original oil in the Tensleep Res- 
ervoir of the Elk Basin Field. This is substantiated further by 
a study of the properties versus structure of the original reser- 
voir oil in the Weber Reservoir of the Rangely Field, Colo. 
The Weber sandstone at Rangely and the Tensleep sandstone 
at Elk Basin are geologically equivalent formations. Labora- 


Table \ 


Composition as a Function of Depth, Based 
on Subsurface Oil Samples from the Weber Sandstone 
Reservoir. Rangely Oil Field, Colo. 


Deptt B 


ft 


itanes and lighter 
ol per cent liquid 


0 14.29 
840 30.57 
840 13.72 


interface 


Pentanes and heavie 
mol per cent liqui 
55.71 
69.43 
Change + 13.72 


Measured from gas-oil 
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VARIABLE CHARACTERISTICS OF THE OIL IN THE TENSLEEP SANDSTONE 


RESERVOIR, 


ELK BASIN FIELD, WYOMING AND MONTANA 


tory analyses of reservoir oil sampled at Rangely provided 
nearly identical change in oil properties, solution gas/oil ratio 
saturation pressure, etc., per foot change of reservoir eleva 
tion with the change reported for the Elk Basin study. Consid 
erably more data on composition of the reservoir oil were 
obtained for the Rangely samples, so that the change in com 
position of the reservoir oil samples might be compared with 
calculated change at conditions of thermodynamic equilibrium 
in the reservoir. 

The change in mol per cent liquid composition of the bu 
tanes and lighter fraction of the reservoir oil calculated by 
use of equations defining thermodynamic equilibrium is shown 
in Table VI. This change will account for only one fourteenth 
of the actual change found by sampling, as shown in Table \ 
For the single component methane in the liquid phase. the 
calculated change was only one thirtieth of the actual change 
found by sampling. It is difficult to reconcile the magnitude 


lable VI—Composition as a Function of Depth, Based 

on Caleulated Effect of Gravitational Concentration of 

Components for the Weber Sandstone Reservoir Oil, 
Rangely Oil Field. Colo. 


and lighter 
ent liquid 


Deptt Butanes 
ft « 


0 14.29 
840 13.35 
Change: 840 0.94 
Measured from gas 
of change of reservoir oil properties with elevation found at 
both Rangely and Elk Basin with conditions of the thermo 
dynamic equilibrium. It is thought that the small changes in 
composition of the reservoir oil defined by thermodynamik 
equilibrium fix the limit toward which the Elk Basin Reservoir 
oil was approaching. but had certainly not reached at the 
time of discovery. 

We believe that the change of 885 psi in bubble point 
calculated by Rzasa from Dominguez crude oil data and due 
to 2,000 ft change in elevation at thermodynamic equilibrium 
is in error and that these data indicate a change of about 70 Ib 

A few the 
hydrogen sulfide content of gas produced from wells as reset 


comments will be made to clarify percentage 


voir pressure declines. In Fig. 6 the upper curve is the per 


centage hydrogen sulfide content of the solution versus 


gas 
ga 


pressure as taken from Fig. 2 


of the subject paper. A reservoi: 
oil sample taken when the reservoir pressure is at any value 
between 1.880 psia (original pressure) and 1.250 psia (satu 
ration pressure), when differentially liberated to atmospheric 
pressure, will yield solution gas of quantity equal to the 
solution gas 
contain 18.2 per cent hydrogen sulfide. A reservoir oil sample 
taken at any lower pressure. for instance, 600 psia will yield 


oil ratio at 1.250 psia pressure. This gas will 


solution gas of quantity equal to the solution gas /oil ratio at 
600 psia, and this gas will contain 23.7 per cent hydrogen sul 
fide. If the produced gas only the 


g solution 
ratio. the percentage hydrogen sulfide would follow 


oil ratio were 
fas oil 
the upper curve of Fig. 6 as the reservoir pressure declines 

The lower curve of Fig. 6 is the hydrogen sulfide content 
of liberated gas. It is the information of Table II of the sub 


ject paper plotted versus mid-point pressures. This curve pro 
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FIG. 6 — ESTIMATE OF PERCENTAGE HYDROGEN SULFIDE CONTENT OF 
PRODUCED GAS FROM WELL 1, ELK BASIN FIELD 


ides the percentage hydrogen sulfide content of a differential 
quantity of gas liberated from the reservoir oil sample for a 
small pressure decline. For example. as the pressure declines 
from 725 psia to 625 psia. gas will be liberated in a quantity 
equal to the difference in the solution gas/oil ratios at the 
two pressures and of approximately 5.9 per cent hydrogen 


-ulfide content. This gas will be liberated in the reservoir. 


If some of this free gas is produced from the reservoir in 
addition to the reservoir oil, the produced gas at the wellhead 
will be a mixture of solution gas from the reservoir oil and 
ree gas from the reservoir. It is obvious that the percentage 
ivdrogen sulfide in the produced gas will be reduced from 
the values of the upper curve of Fig. 6 by the free gas pro- 
duced the 


ince of percentage hydrogen sulfide content of the produced 


from reservoir. A qualitative estimate of the vari- 
gas with reservoir pressure is made in Fig. 6. The values trend 
iway from the upper solution gas curve as free gas produced 
from the reservoir accounts for increasing percentages of the 
produced gas oil ratios. This will occur as reservoir pressures 
are reduced from saturation pressure and will continue until 
the produced gas /oil ratios reach a maximum value. As the 
pressure is further reduced and the produced gas/oil ratios 
decline, the values will trend toward the upper “solution gas” 
curve, It is true that for the period of production character- 
ized by extremely high produced gas/oil ratios the values 
will be nearly equal to the lower or “liberated gas” curve or 
the values of Table III of the subject paper. For the initial 
portion of the pressure decline the values will be more nearly 
equal the values of per cent hydrogen sulfide of the “solution 


gas” curve 
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RELATIVE PERMEABILITY TO LIQUID 
IN LIQUID-GAS SYSTEMS 


L. A. RAPOPORT AND W. J. LEAS, MEMBER AIME, THE CARTER OIL CO., TULSA, OKLA. 


ABSTRACT 


As a preliminary, consideration is given to the conventional 
definition of relative permeability and to the conditions gov 
erning the simultaneous flow of oil and gas through porous 
media. 

For the conditions of flow prevailing throughout most of a 
gas drive reservoir, the oil and gas can reasonably be sup 
posed to be in capillary equilibrium with each other. Under 
these conditions, and these conditions only, the relative per- 
meability to liquid can be expressed as a function of satura 
tion. The relative permeability to liquid in that case is de 
pendent upon the distribution of fluids which itself is shown 
to be related to the capillary pressure, and, in turn, to the 
saturation. As a consequence, relative permeability to liquid 
can be expressed in terms of the volume and surface area of 
a network of liquid channels bounded by the rock and the 
gas phase. While the volume of this network can be evaluated 
accurately, the surface area cannot. However, for any such 
volume, maximum and minimum values of the corresponding 
surface area can be calculated from capillary pressure data. 
It is then possible to establish for any saturation the limits 
within which the value of the relative permeability to liquid 
must lie. 

As a consequence of the theoretical development, the valid 
ity of an experimental method for measuring relative perme 
ability to liquid which utilizes a stationary gas phase is dem 
onstrated. In this method capillary barriers are cemented to 
the ends of the core sample to permit the maintenance of 
capillary equilibrium between the two phases. At the same 
time, this procedure eliminates undesirable secondary phe 
nomena such as end effects, fissure effects, etc., the presence 
affect the results of laboratory 


of which adversely other 


methods. 

The results obtained by theoretical calculations, and experi 
mentally, are discussed. In view of the overall precision that 
can presently be obtained in reservoir calculations, the agree 
ment between the calculated and measured relative perme 
ability to liquid data can be considered satisfactory. 


In conclusion, for reasons of economy and simplicity, the 
procedure of calculating limiting relative permeability to 
liquid curves from capillary pressure data is indicated for 
general engineering purposes. It is shown that the above 
procedure can easily be extended to the cases where connate 


‘References given at end of paper. 
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water is present. Its use for reservoir studies is particularly 
recommended in conjunction with the method for measuring 
relative permeability to gas' which simultaneously yields the 
capillary pressure data necessary for the calculations. 


THEORETICAL 
Definition of Relative Permeabilities — 
Basic Equations for Heterogeneous Flow 
The equations by which the relative permeability concept 
is defined and upon which the formulation of all of the gas- 


oil flow problems rests at the present time are expressed as: 


_kK, 


By 


k 
J Grad P Grad P, (la) 


Grad P, (1b) 


kK 
Grad Py = - 
“ 


and ,, refer to liquid and gas: V is the volumetric 
rate of flow per unit gross area, « the viscosity, Grad P the 
potential gradient. and & the specific permeability of the 
porous medium.* (For horizontal flow, Grad P becomes the 
pressure gradient: otherwise, gravity must be included.) Ac- 
cording to these expressions, each of the constituent phases 
s considered similar to a homogeneous system where the 
volumetric rate of flow is proportional to the pressure gradi- 
ent, and for each of which the constants of proportionality. 
k, and ky, are termed effective permeabilities, by analogy to 
the specific permeability as defined by Darcy’s law in its 


where 


original form. In order to obtain a convenient basis of com- 
parison, the effective permeabilities are referred to the spe- 
of the considered porous medium, with 
the help of the relations: 

k, =K, h 

Re A Be © os. sk a> nea ee ee 
K, and Ky are defined as the relative permeabilities to the 
liquid and to the gas phase respectively, and frequently ex- 
pressed in per cent of specific permeability. It may be seen 
that Equations (1), which appear to be a direct generalization 
of Darcy’s law, correspond to the assignment at any given time 
of a set of “local” permeabilities to each point of the porous 
medium, and represent in a differential form the two fluid 
flow system as a simple superposition of the individual single 


cific permeability, k, 


flow systems 


The above interpretation implies that the effective or rela- 
tive permeabilities are independent of pressure or rate of flow. 


*Nomenclature at the end of the paper 
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an assumption which appears to be substantiated by most of 
the experimental results. It seems reasonable, therefore, to 
consider at first the relative permeability concept on the 
basis of the conventional assumption (suggested by experi 
mental observations), that relative permeabilities are functions 
of saturation only. The theoretical and experimental methods 
presented in this paper are confined to liquid-gas systems, for 
which it can be shown that each phase taken separately satis- 
fies the conditions under which Darcy’s law applies. 


ANALYTICAL EXPRESSION FOR SPECIFIC 
PERMEABILITIES 
Until recently the specific permeability, A, such as defined 


by Darcy's original relation 


k Grad P 

had to be considered as a purely empirical characteristic. A 
theoretical expression for Darcy’s law and an analytical inter- 
pretation of the specific permeability were established by J. 
Kozeny in 1926. His theory is essentially based upon a statis- 
tically homogeneous and isotropic distribution of the pore 
voids. It excludes any other simplifying assumption such as 
representing the porous medium by bundles of capillary tubes 
or by layers of spheres, ete., and leads to the final equation 


l f 
‘ - Grad P_. (4) 

bu tA 
where A is the interstitial surface area per unit of bulk vol- 
ume, / the porosity, and ¢ a textural constant. By identification 
ot Equations (3) and (4), it is possible to express the total 
permeability as 

1 f 

t A 
\ccording to theoretical. as well as numerous direct experi- 
it appears possible to assign to the textural 


A 


mental results,’ 
constant, t, the average value of 5, so that Equation (5) can 


be rewritten as 


(6) 


Up to the present time. direct experimental checks of the 
Kozeny equation had been carried out only for unconsolidated 
granular beds. However, with help of the thermodynamical 
approach presented in Appendix A, a more complete verifica- 
tion of the theory has been made. (See Appendix B.) It was 


found that on the average the values of permeability computed 


according to Equation (6) were within +12 per cent of the 


corresponding experimental measurements. These results were 


AREA OF SOLID IN 
CONTACT WITH 

AS PHASE 
AREA OF SOLID IN 
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BETWEEN TWO FLUIOS 


FIG. 1 — HYPOTHETICAL REPRESENTATION OF FLUID DISTRIBUTION IN 
RESERVOIR 
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established over the widest possible range of homogeneous 
core material, with permeabilities varying from 15 md to 
17,000 md and porosities varying from 9.3 to 35 per cent. Its 
consistency can, therefore. be considered as a good confirma- 
tion of Kozeny’s theory, as well as of the constant value to 
be assigned to the textural constant. 

The important physical concept to be derived from Kozeny’s 
theory is that the porous medium must essentially be consid- 
ered as a “flow matrix” of irregular shape. All these irregu- 
larities do, however, statistically average so that the perme- 
ability is determined only by the inside area and volume as 
indicated by Equation (5). It is then possible to define the 
following conditions for the general application of Darcy’s 
law, and for its extension to multiphase systems: 

Continuity of the pressure gradient, which requires the 
continuity of the flowing phase throughout the system. 


(a) 


(b) Microscopically steady flow pattern, i.e., laminar flow 
as opposed to turbulent or oscillatory regimes. 
Existence of a flow matrix defined by its inside surface 
area and volume per unit of bulk volume. 


CORRELATION OF THE RELATIVE PERME- 
ABILITY FUNCTIONS AND FLUID 
DISTRIBUTIONS 


Fluid Distribution at Static Equilibrium 


Since the specific permeability depends exclusively upon 
the structural properties of the porous medium, the study of 
relative permeabilities must include consideration of the spa- 
tial distribution of each of,the phases. 

When two fluids in a given proportion, or under given exter- 
nally imposed conditions occupy a porous medium, their dis- 
tribution at static equilibrium is determined by surface energy 
relations. These surface energy relations are best visualized 
with the help of the phenomena pertaining to interfacial ten- 
contact angle. The contact angle determines the 
curvature of the interface between two immiscible fluids 
which, under the influence of the interfacial tension gives rise 
to the capillary pressure. Under static conditions, stable equi- 
librium is attained when the external forces such as gravity, 
or artificially imposed pressure, upon one of the phases are 
balanced by those of capillary pressure. For a given capillary 
pressure (i.e., to that which is required to balance the exter- 
there is a corresponding curvature, the same for 


sion and 


nal forces) 
all the existing interfaces. Considering then the geometry of 
the porous medium, it is seen that for a given curvature, there 
must correspond a certain position of the interfaces, i.e., a 
definite microscopic fluid distribution and a definite saturation. 
\s a consequence, it can be anticipated that inside the porous 
medium each constituent phase has a definite geometrical form 
determined by the saturations.* 

In order to obtain quantitative results based on the above 
descriptive picture, the thermodynamic aspects of the capillary 
equilibria must be considered. Stable equilibrium at a given 
saturation corresponds to the minimum level of free energy 
in this case free surface energy —that can be reached. The 
first condition for stable equilibrium is that each of the con- 
phases be continuous. Furthermore, in a two fluid 
system, the “wetting” fluid is the one which has the smaller 
tension of adhesion against the rock. Consequently, at stable 
equilibrium the wetting fluid has to be distributed in such a 


stituent 


an exist only under conditions of stable equilibrium, 
possibility of defining capillary 


elationship « 
phenomena exclude the 
manner 


*Such a 
hyst 


pressure in 


since 
unequivocal 
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way as to cover always the greatest possible surface of solid. 
On the other hand, the non-wetting fluid will be distributed in 
such a manner as to cover the smallest possible surface of 
the solid, and, at the same time, to expose the smallest pos- 
sible interfacial area to the wetting fluid. It may then be seen 
that for a gas-liquid system the liquid having a tendency to 
preserve the greatest contact area with the solid will always 
fill the smaller pore openings, whereas the gas remains in 
the larger ones. 
Fig. 1 illustrates in a schematic way the areas referred to 
above for which the following relationships can be established: 
A, = R,, +1 = Total surface area of the liquid (wetting 
fluid) system. 
4, = R, +1 = Total surface area of the gas (non-wetting 
fluid) system. 
A = R,, + R, = Total surface area of the solid system, i.e., 
total grain surfacen. . .... =. (7) 
where R,, represents the contact area between the solid and 
the liquid, Re the contact area between the solid and the gas 
and / the interfacial area between the two phases. 


Extension of Static Fluid Distribution to Flow 
Conditions; Definition of Parallel Flow 


At static equilibrium, the distribution of two fluids in a 
porous medium depends essentially on the position of the 
interfaces, and as long as the pressure drop across a given 
interface remains constant, neither its shape nor its position 
will be changed. If, starting from a given static capillary 
equilibrium, flow is produced, a pressure gradient giving rise 
to this flow is established in each phase. For each interface 
this “flowing” pressure gradient may be considered as a local 
disturbance that is superimposed on a previously existing 
static equilibrium of capillary forces. Hence, the pressure drop 
across each interface, instead of being constant, over the entire 
interface varies from one point to another (see Fig. 2). Each 
interface will, therefore, have a tendency to modify its position 
in such a way as to equilibrate the resultant effect of the 
previously existing static pressure and the superimposed pres 
sure variation corresponding to the flowing pressure gradients. 
Under conditions of steady state flow, it may be visualized 
that each interface will establish itself in a position of “dy 
namic equilibrium.” Theoretically, this position is different 
from that of static equilibrium; however, for the purpose of 
analysis it is admissible to treat that type of flow for which 
the “disturbing” action of the flowing pressure gradient is 
negligible. Such an approach is justified whenever the flowing 
pressure gradient is sufficiently small in comparison to the 
capillary pressure, which is actually the case under most of 
the practical field conditions. 

The flowing pressure gradients that usually prevail over 
the greatest part of natural reservoirs have a maximum of 
some pounds per foot, for instance 10 psi/ft. The dimensions 
of the interfaces can be considered as being of the same order 
of magnitude as those of the pore openings, i.e., 1/10 mm as 
a maximum (in homogeneous media). That means that the 
change in capillary pressure, dP., resulting from the flowing 
pressure gradient over the dimensions of a given interface 
will in the average not exceed 1/300 psi. The capillary pres- 
sure itself, however, that establishes across the interface, is 
usually of the order of one psi and certainly no less than 0.3 
psi. Consequently, the ratio dP./P. which indicates the rela 
tive importance of the local disturbance, can always be ex 
pected to be smaller than 0.01 and may, therefore, be neg 
lected as far as the micromechanism of flow is concerned. 

For reasons of simplicity, one of the fluids had been assumed 
to remain stationary. In most cases, however. both fluids are 
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FIG. 2— EFFECT OF FLOW ON THE EQUILIBRIUM OF AN INTERFACE. 


flowing in the same direction, so that the actual dynamic 
disturbance is due to the difference of the flowing pressure 
gradients and might be still smaller than indicated above. 

The above example shows that when the flowing pressure 
gradients are an order of magnitude smaller than the capillary 
pressure, and when the flow originates from a condition of 
stable (static) capillary equilibrium, the distribution of the 
flowing fluids remains stable and the same as at static capil- 
lary equilibrium. The movement of the . ases can then be 
visualized as if they were flowing in an independent manner 
through two interwoven but distinct networks of channels. 
These networks, bounded partly by rock and partly by the 
other fluid, can be considered as extending “in parellel” 
throughout the system. Such a situation will be referred to 
as “parallel flow.” It is a generalization of the concept of 
laminar flow in hydraulics. which defines a steady state stream- 
line pattern. 

In the reservoir (and in many laboratory experiments), the 
liquid saturation generally decreases in the direction of flow: 
furthermore, there is usually an overall (slow) decrease of 
saturation with time. However, once it has been recognized 
that the perturbation due to flow is small enough so as not to 
modify locally the nature of the capillary equilibrium phe- 
nomena, one can consider this perturbation as a differential 
element causing the saturation and fluid distribution to vary 
as continuous functions of space and time. It is then always 
possible to divide the total space into sufficiently small por- 
tions for each of which rigorously parallel flow can be ad- 
mitted, if sufficiently small time intervals are considered. It 
is concluded that any heterogeneous flow system can be treated 
as a succession in time and space of “elementary” parallel 
flows, if throughout the system the saturation varies in a 
manner (with time and space) and capillary 
equilibrium is preserved. 


continuous 


ANALYTICAL EXPRESSION FOR RELATIVE 
PERMEABILITY TO LIQUID 


Referring to the conditions of applicability of Darcy’s law, 
it is readily seen that in the case of parallel flow, such as 
described above, Darcy’s law can be extended to each of the 
constituent phases; hence, Equation (5) can be applied to 
the liquid phase in the following manner: 


(a) The replacement of the absolute porosity, f, by Sf—the 
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FIG. 3 — EXAMPLE OF CAPILLARY PRESSURE CURVE 


liquid filled pore volume per unit of bulk volume, i.e., by 
the “effective porosity” with respect to the liquid. 


The substitution of the inside pore area, or total grain 
surface, A by A,, which represents the total boundary 
area of the liquid phase “flow matrix.” 


By substitution of the above terms into Equation (5), the 
effective permeability to the liquid phase is obtained 
Ss | 
hk > teih be re wae fe 
tA 
and division of Equation (8) by Equation (5) gives the ex- 
pression for relative permeability to liquid 


Ss 


(|) 
{ 

In order to obtain a more usable expression for the relative 
permeability to liquid, the following considerations must be 
taken into account. It is generally recognized that toward the 
end of a capillary displacement process an irreducible mini- 
mum saturation, S,. is reached, and no more liquid can be 
forced out, except by diffusion (the effects of which will be 
considered as negiligible in the present treatment). At this 
saturation the effective permeability to liquid becomes zero. 
Therefore, it seems logical to treat in first approximation the 
irreducible liquid as an essentially stationary element which 
reduces the porosity of the porous medium as well as the vol- 
ume of flowing liquid. Once this liquid is considered as part of 
the solid framework, it is necessary to use an “effective” area, 
{,. representing the surface that separates the flowing fluid 
from the immobile flow matrix composed of solid and station- 
further 


K 


ary fluid, instead of the specific solid surface. A: 
more, instead of using the total liquid saturation, S. a “reduced 
saturation.” S’. has to be introduced. 


a 
(10) 
nn 
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On substitution into Equation (9) the following expression 
for the relative permeability to liquid is then obtained: 


S* ( S-S, )( Ay ) 
(4 ) 1-S, A, 
{ 
It will be noted that the derivation of Equation (11) im- 
plies the two following assumptions: 


(11) 


(a) The molecular layers at the interfaces between solid /li- 

quid and gas/liquid are stationary. If they were not, the 
basic Kozeny treatment should be slightly modified so as 
to take into account a certain “slippage” of the fluid on 
part of its boundaries. 
The textural constant of the wetting fluid “matrix” is the 
same as that of the total pore space. Such an approxima- 
tion seems to be justified since in general the possible 
variations of the textural constant have been recognized 
to be comparatively small 


EVALUATION OF LIMITING RELATIVE 
PERMEABILITY TO LIQUID FUNCTIONS 
FROM CAPILLARY PRESSURE DATA 


Che calculation of K,, the relative permeability to liquid. 
requires a knowledge of the surface areas, 4, and A,. With 
the help of a thermodynamic approach indicated in Appendix 
\, the following expressions can be derived from a liquid/gas 
capillary pressure curve: 

S 
f P.ds 
l 
S, 
f P ds 
1 
S. 


4 f P.dS + 2 


+ Ss 


f x, + 2 (14) 
in which the integrals indicate the areas under the P. vs. S 
curves measured between the indicated limits (see Fig. 3). It 
is seen that a system of three equations is obtained for the 
unknowns A,. Ay. A,. and /. Consequently, while the 
values for A, and A, are directly indicated. separate solutions 
cannot be defined for A, or /. 


it is necessary to introduce an addi- 


four 


In order to evaluate A,. 
tional assumption concerning the distribution of the fluids. 
Such an assumption, expressed in terms of Ry. R,. and J, will 
then represent one more independent equation, which together 
with the above relations will furnish a system from which all 
the unknowns can be solved. It might be recalled that on the 
basis of thermodynamic considerations only a more or less 
complex “trend” could be indicated for the liquid distribu- 
tion, namely maximum contact area with the solid compatible 
with a minimum amount of interfacial area. In order to ex- 
press such a trend in a more exact manner. it would be nec- 
essary to consider in detail the geometry of the porous medium. 
which depends on very many parameters such as grain size 
distribution, sphericity or shape factor of the grains, packing. 
Obviously. at the present time, it is impos- 
sible to treat such a geometrical complex in a general and 
precise way. Therefore, no exact additional relation can be 
established that would lead to unique solutions for A, and /. 
The only logical approach that can be attempted consisis of 
establishing the most general limiting conditions of fluid dis- 
tribution that could be expected to apply to any. or at least 


consolidation, ete. 
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to most, of the possible geometrical systems. According to 
such a principle, separate calculations for the minimum and 
for the maximum values of K, will be shown. 


Evaluation of Limiting K,, Minimum Function 


In accordance with the previously established relations. 
Equations (7), (12), and (14), the following system of equa 
tions can be considered: 

2, = A, - 2! . (15) 
2,=R +1 via ts “a & «6 (16) 
where 2, and =, have a well defined value at any saturation 
and can, therefore. be considered as “constants” in regard to 
the unknowns 4,, /, and Rg. It is readily seen that according 
to Equation (15) for any set of =, and =,, i.e., for any satu 
ration, a maximum value will be obtained for A, if / is maxi 
mum and that according to Equation (16) the maximum for / 
obtains if R, is assumed to be zero. Consequently, the maximum 
value that A, may take at any saturation is 
{, (Max.) = 2, +2 2, = 2s, + 2s 
According to Equation (11) it is clear that at any saturation 
the smallest possible value will be obtained for the relative 
permeability to liquid, K,. if a maximum area, 4,, is consid 


ered, so that 
1. P : Zs 
— . (18) 
1, (Max) Zs, + = 


Physically the derivation of K,(Min) means that the dis 
tribution of the fluids is such as not to permit the existence 
of any contact between the gas phase and the solid surface 
The gas can then be visualized as flowing inside a network of 
channels completely surrounded by liquid, as schematically 
represented in Fig. 4, and it might be noticed that such an 
assumption concerns only the disposition of the fluids, and 
does not correspond to any hypothesis about the geometry of 
the porous medium itself. 


(17) 


K, (Min) = S’ 


Evaluation of Limiting K, Maximum Function 


Referring to the system of Equations (15) and (16), it is 
seen that the theoretical absolute minimum value for 4,, would 
obtain under the assumption / = 0. i.e., A, = Zy. Such a 
situation would correspond to Rg = Ag, i.e., to the assumption 
of having the gas phase —and equally the liquid phase 
bounded exclusively by solid surfaces. In that case, the porous 
medium should be visualized as formed by non-interconnected 
capillaries. It is clear that such an assumption is by far too 
limiting and leads actually to K, values which are too high 
Thus, instead of considering — even for the purpose of limit 
ing conditions — a porous medium as equivalent to a bundle 
of capillary tubes, it is reasonable to represent it as an iso 
tropic random packing of grains. Under this more general 
assumption, interfaces must exist between the gas and liquid 
and none of the three surface areas, /, R,. R, is zero. The 
lowest possible value of the ratio //R,, can then be evaluate: 
The use of this value in Equation (16) leads to a minimum 
value of /, of A,. and consequently, to a maximum of K,.* 

With the help of general capillary pressure relationship 
it is possible to show that 


=>) 
a er. ee ob ek) 


where P. represents the capillary pressure at any saturation 
S, and P,, the mean value of P. corresponding to the satura 


(19) 


*The detailed treatment for the derivation of K, maximum is indicated 
in Appendix C. 
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tion range S to S,. The above inequality leads to the following 
minimum value of 4,: 


9 


1-f & 7 
l + 

Sf Pes 
and the final expression for the limiting maximum relative 
permeability to liquid is then obtained as 


1, (Min) (20) 


K, (Max) 


EXPERIMENTAL 
Principle of the Method 


The experimental method has been conceived so as to meas- 
ure the relative permeability to the liquid phase under the 
conditions of parallel flow defined previously. According to 
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these conditions, capillary equilibrium must be maintained be- 
tween the two phases throughout the entire core specimen. 
The main difficulty in obtaining such an equilibrium results 
mostly from different sorts of “end-effects,” the influence of 
which extends back into the core. In the experimental method 
presented, the capillary pressure between gas and liquid is 
established by means of semi-permeable barriers, and con- 
trolled by the pressure imposed upon the gas. The elimination 
of secondary phenomena as outlined above is obtained by 
cementing the capillary barriers, permeable to liquid only. to 
each end of the core sample. The experimental procedure can 
then be outlined as follows: 


(1) The ensemble of core and end plates is fully saturated 
with liquid. 


(2) Through a central orifice in one of the end plates, gas is 
forced into the core and liquid is displaced through both 
the end plates until the system comes to equilibrium, i.e., 
reaches a saturation which corresponds to the capillary 
pressure that has been imposed upon the gas. 


The liquid is then allowed to flow through the system, 
whereas the gas remains immobile, held back by the 
semi-permeable barriers. The flowing pressure gradient 
in the liquid is kept constant and sufficiently low so as not 
to disturb the capillary equilibrium and the previously 
established saturatiop. 


When perfectly steady state conditions are reached, the 
exact values of the pressure drop in the liquid over the 
length of the core sample and the rate of flow are deter- 
mined. The effective and, consequently, the relative per- 
meabilities are then obtained and plotted against the 
corresponding saturation. 


Description of Experimental Apparatus 
and Procedure 


The relative permeability to liquid cell is shown in Fig. 5. 
In general construction it resembles the relative permeability 
to gas cell’ and many of the parts are interchangeable. Not 
shown in the diagram are the “O” ring seal between tube T3 
and the porous plate, and the grooves in the face of the metal 
plates B] and B2 for the distribution of liquid to the porous 
plates. These grooves ensure that the liquid phase is free to 
move unidirectionally through the core, thereby eliminating 
any divergent or convergent flow except in the immediate 
vicinity of the gas inlet tube T3. The ratio of the area of the 
end of this tube to the cross sectional area of the core is 1/64 
so that any resulting perturbations in the flow pattern are neg- 
ligible. Gas is admitted to the core through the tube T3, which 
is connected to a gas pressure source and to a manometer 
with which the pressure of the gas phase is measured. Liquid 
enters the cell through tube T1 and leaves through tube T2. 
These tubes are connected through a commutating valve to 
the flowmeter. 


The flowmeter consists of two precision-bore calibrated glass 
tubes which are side by side in the same horizontal plane. 
The right end of one tube and the left end of the other are 
connected to a commutating valve. Liquid fills only a por 
tion of each tube and the application of gas pressure to the 
gas-liquid interface in one of the tubes causes the liquid in 
that tube to flow through the core into the other tube. When 
the first tube is nearly empty of liquid, the gas pressure to 
that tube is released, the commutating valve is turned, and 
the gas pressure is applied to the second tube. which is now 


88 PETROLEUM TRANSACTIONS, AIME 


RELATIVE PERMEABILITY TO LIQUID IN LIQUID-GAS SYSTEMS 


nearly full of liquid. It, in turn, is emptied and the liquid 
reappears in the first tube. This cycle can be repeated indefi- 
nitely and though sizable volumes of liquid can be passed 
through the core, only a few milliliters of liquid are required. 
The fact that the liquid-gas interfaces in the tubes are in the 
same horizontal plane eliminates any gravitational gradient in 
the liquid phase, and only the differential gas pressure acting 
on the interfaces needs to be measured to determine the pres- 
sure gradient causing flow in the liquid phase. 


The commutating valve permits the reversal of fluid flow in 
the flowmeter tubes while maintaining at all times the same 
direction of flow within the core. A bleeder valve in one of 
the flowmeter lines permits the adjustment of the volume of 
liquid within the system. Since the two flowmeter tubes are 
parallel and of the same internal diameter, the two interfaces 
move at the same rate and in the same direction when liquid 
is flowing through the core under steady state conditions. By 
measuring the velocity of the interfaces and the differential 
gas pressure acting upon them, the effective permeability of 
the core under test can be calculated. 


The flowmeter also indicates the amount of liquid expelled 
from the core. At the start of the relative permeability meas- 
urement the liquid saturation of the core is unity and the dif- 
ferential linear separation of the interfaces is noted. As long 
as the liquid saturation of the core remains fixed, the separa- 
tion of the interfaces remains constant. When the gas pres- 
sure applied through the gas inlet tube T3 exceeds the thres- 
hold pressure of the core, liquid is expelled from the core 
and appears in the flowmeter tubes as a change in the separa- 
tion of the interfaces. This permits the calculation of the 
liquid saturation of the core at any time. A scheme similar to 
this has recently been described.” 


Several factors must be considered in the selection of the 
plates which function as semi-permeable barriers. Firstly, the 
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pore size of the plates must be such that they remain fully 
saturated with liquid at all gas pressures to which they are 
exposed. Secondly, they must have as high a transmission co 
efficient to liquid as is possible (in conformity with the first 
requirement) so that changes in the relative permeability of 
the core will not be unduly suppressed. This can, in part, be 
accomplished by using thin plates, yet they must remain thick 
enough to withstand breakage. The elimination of undesirable 
“fissure effects” is accomplished by cementing the plates to 
the core. 


After the core to be tested has been cut to shape and ex 
tracted, its pore volume and permeability (to n-decane) ar 
measured. It is dried and the porous end plates are cemented 
in the following manner. The core and plates are fully satu 
rated with water. A plaster of Paris slurry is spread over one 
side of the outlet plate, which is pressed against the end of 
the core to expel the excess plaster of Paris.* The core with 
outlet plate is dried and its permeability to n-decane is meas 
ured. It is redried and the inlet plate is cemented in the sami 
manner as was the outlet plate. The ensemble is then put into 
the cell where its liquid permeability is determined at a liquid 
saturation of unity. 


Air, under pressure, is then admitted to tube T3 and a por 
tion of the liquid in the core is expelled into the flowmeter 
When equilibrium is reached, the differential separation of the 
interfaces in the flowmeter tubes becomes zero and the per 
meability of the ensemble at a liquid saturation less than unity 
is measured. The air pressure in T3 is again increased and 
the cycle is repeated until the permeability of the ensemble 
falls to a value too low to be measured. From these measure 
ments the effective permeability of the core and end plates is 
obtained, and the relative permeability can then be determined 
for each corresponding liquid saturation. 


*This procedure of saturating the core prior to affixing the plates mir 
mizes the possibil.ty of plugging the core through movement of ti 
plaster slurry into the core if it were applied to a dry core. The liquik 
permeability of a core whose porous plates were attached in this mann« 
was found to be unchanged after their removal 
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RESULTS AND DISCUSSION 


Results 

The results obtained with n-decane and air using three alun- 
dum cores and three natural cores are shown in Figs. 6 to 11. 
The solid curves are the limiting K,, values calculated from 
displacement data obtained on these cores using a single-core 
displacement cell” with decane displaced by air. The limiting 
curves for the three alundum cores are much the same. This 
is to be expected since their capillary pressure curves were 
similar. With the exception of two points on Core 59A-15, the 
data fall within the limiting values. The agreement between 
the calculated and experimental data for the alundum cores 
is good for each core. 

In considering the three natural cores, it is seen that the 
limiting K,, values vary widely both as to position and spread. 
Here again their position is determined by the shape of the 
capillary pressure curve. Generally, it has been observed that 
a capillary pressure curve of near zero slope which undergoes 
a sharp transition to a vertical slope at the minimum satura- 
tion results in a large spread between the calculated K, 
limits, whereas capillary pressure curves which have an ap- 
preciable slope even at high liquid saturations usually result 
in a very narrow spread of the limiting K, values. Consider- 
ing the narrow spread, the agreement is good for the Bohanon 
core. Three runs were made on the Elk Basin core and the 
results cluster about the lower limiting curve. No single run 
has points which lie either entirely above or below the lower 
curve. The Bartlesville core has narrow limits and the experi- 
mental data which seem to be slightly S-shaped fall partly 
outside the limits. On the whole, the agreement obtained with 
the alundum cores is better than that obtained with the nat- 
ural cores. Although the exact reasons for this are not known, 
it is possible that the higher degree of homogeneity of the 
alundum cores is a contributing factor. 

Although K obtained by the Botset" and “Penn 
State”” methods might be regarded as questionable,” their 
rather mention. The results ob 


results 


extensive use warrants their 
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tained by these methods on an unconsolidated sand and on a 
synthetic sandstone are shown in Figs. 12 and 13. The experi 
mental points show greater scattering, which might be ex 
pected since no difference in pressure exists between the 
phases, so that appreciable “end effects” can be anticipated 
The majority of the points, however. lie within the calculated 
limits. 


DISCUSSION OF THE RESULTS 
Validity of the Parallel Flow Concept 


The theoretical and experimental methods presented in this 
paper for the determination of relative permeability to liquid 
are confined to the strict consideration of parallel flow mech- 
anisms, i.e., to the case where the flowing fluids remain in 
capillary equilibrium with each other. In most practical in 
stances, it is recognized that parallel flow may be “perturbed” 
by secondary phenomena, mainly the release of gas bubbles 
Such bubbles do not correspond to capillary equilibrium, and 
furthermore, cannot be subject to laminar flow. It can be 
anticipated, however, that in any point of the reservoir the 
existence of separate gas bubbles corresponds to a relatively 
short transient stage. Due to rapid diffusion effects and pos 


sible mechanical transportation, the separate bubbles have the 
tendency to reach stable capillary equilibrium and to merge 


into a continuous gas body. Recent experiments performed 
in this laboratory have indicated that (under atmospheric con 
ditions) the time required in order to reach phase continuity 

and consequently capillary equilibrium 
of 24 hours. 


It is clear that the formation of a continuous gas phase be 


is of the order 


comes possible only after a certain minimum amount of 


bubbles, corresponding to the “equilibrium gas saturation” 


has been released from solution. The phenomenon of equi 
librium gas saturation must be considered as a determining 
factor for the flowing behavior of the gas phase, and it is for 
that reason that a theoretical derivation of the relative perme 
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ability to gas function has not yet been presented. For the 
liquid phase, the continuity of which can always be expected 
to be uninterrupted, it is however, a reasonable approximation 
to derive the relative permeability function under the assump- 
tion that capillary equilibrium exists throughout the entire 
range of saturation. Accordingly, for most practical purposes 
the use of the parallel flow concept and the application of the 
presented methods appears to be justified in the case of any 
liquid-gas flow system. 


Discussion of the Experimental Method 


In the presented experimental technique, undesirable fea- 
tures of other methods have been overcome. The use of ce- 
mented porous end plates permits the maintenance of a pres- 
sure difference between the gas phase and the liquid phase. 
and end effects and fissure effects. These latter 
correspond to a progressive invasion by the gas phase of the 
microscopically large voids which exist between the core and 


eliminates 


porous plates (or two adjacent core sections) as shown by 
direct-X-ray measurements.” Furthermore, there is no distor- 
tion of the flow pattern from linear flow since at the ends of 
the core there are no pressure taps or sizable areas for gas 
flow. The determination of the pressure gradient is subject 
to unequivocal measurement and is effective over the complete 
length of the core. 

In the tests run by this method, the pressure gradients did 
not exceed a few centimeters of water. Under such conditions 
the use of an immobile gas phase is justified insofar as the 
theoretical discussions have shown that the relative perme- 
ability to liquid so obtained can be expected to be the same 
as though both phases were flowing. 

It is recognized that due to the immobility of the gas. a 
saturation distribution exists throughout the core. The proce- 
dure of plotting directly the measured relative permeability 
of the core against its average saturation might, therefore, be 
systematic error. However, the small 


subject to a pressure 
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gradients and rates of flow used in the present experiments 
preclude a large correction. To verify this, the points x in 
Fig. 10 were obtained using a pressure gradient which was 
twice that used for the other points. These points therefor: 
correspond to a steeper saturation distribution. Since, however 
the agreement between all the runs is very close, it appears 
that at the rates of flow used. the effect of any 
distribution is very small. 


saturation 


Under identical experimental conditions, two separate run- 
were made using cores 59A-15, Fig. 8, and EB-42, Fig. 10. Th 
results of the first run are indicated by circles and the second 
by triangles. All points lie within a narrow band no wider 
than would be expected on the basis of normal experimental 
deviations, so that the reproducibility of the method can be 
considered as satisfactory. 


Validity of the Calculated Data 


Since the limiting K, functions are expressed essentially by 
the cube of the “reduced” saturation, S‘ = S—S,/1-—S.,, they 
are very sensitive to the variations of the parameter, S,. in 
the regions of low saturations, i.e., when S approaches S,. The 
determination of the proper irreducible minimum saturation, 
S,, to be used in the calculations is comparatively delicate in 
view of diffusion effects, which can be quite important in 
liquid/gas systems. Both calculated limiting K,, curves might 
therefore, be subject to a systematic error, namely too high 
in the regions of low saturation, since the diffusion effects 
cause an apparent decrease of S,. This effect is apparent in 
Figs. 8 and 10, and has been observed in several othe: 
instances. 

The limiting K,, curves presented in this report were calcu 
lated from decane/air P. curves. Although this is justified for 
the purpose of checking the calculated results against the di 
rect K, measurements (which were obtained on decane-air 
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systems), the procedure should be slightly modified for prac- 
tical application. It is necessary to take into account the 
possibility of swelling clay minerals frequently contained in 
natural cores. Strictly speaking, the calculations should, there- 
fore, be based upon P. curves obtained by the displacement 
of decane (or any other suitable hydrocarbon) by air in the 
presence of connate brine. However, it can be anticipated that 
sufficiently reliable results will be obtained by simply using 
the brine/air displacement curves that are usually available. 
Justification for such a procedure is indicated by recent ex- 
perimental determinations of relative permeability to gas, 
which yielded practically identical results with and without 
having a “connate” liquid phase present.’ 

The fact that instead of calculating a unique A, function, 
only limiting curves are determined, might be examined in 
the light of the overall degree of precision that can be obtained 
by the present techniques. In many cases it can be observed 
that the spread between the experimental data is of the same 
order of magnitude as that of the theoretical curves. Further- 
more, in view of the sampling problems associated with petro- 
leum reservoirs, the present need is actually for approximate 
data on a large number of reservoir cores rather than precise 
data on only a few. In many practical instances, the use of 
the theoretical limiting K, curves might, therefore, be desir- 
able. Engineering calculations can be carried out either by 
using separately the upper or the lower curve, or by con- 
and determining (by comparison 
maximum error introduced. 


structing. the average curve 
to the limiting curves) the 


CONCLUSIONS 


From the work the following conclusions are 


drawn: 


foregoing 


(1) By means of a theoretical analysis of fluid micro-behavior, 
it has been shown that the relative permeabilities needed 
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for studies of liquid-gas reservoirs are directly related to 
the distribution of fluids and must primarily be functions 
of saturation. 


A procedure for calculating limiting K, curves from 
capillary pressure data has been established. The pos 
sibility of applying the calculated curves to the case 
where connate water is present has been indicated. 


An experimental method based upon the use of cemented 
capillary barriers and upon the principle of an immobile 
gas phase has been developed for measuring the relative 
permeability to liquid. It has been shown that in the 
above method the undesirable phenomena such as end 
effects, fissure effects, etc.. common to other types of 
procedures, have been eliminated. 


The general validity and the comparative simplicity of 
the experimental method have been demonstrated. It has 
been found that good agreement exists between the experi 
mental and calculated data. 


For reasons of economy and simplicity, the calculation 
of limiting K, curves seems indicated for routine engi- 
neering purposes. It is particularly recommended for use 
in conjunction with a method developed previously for 
determining relative permeability to gas.’ which yields 
K,, and P.. data simultaneously. 


ACKNOWLEDGMENT 


The authors wish to thank the Carter Oil Co. for permission 
to publish this paper. They are particularly indebted to W. A. 
Bruce and Charles D. Russell for their valuable suggestions 


and 


helpful assistance during the course of this research and 


in the preparation of this report. and to Miss J. J. Hume for 
much of the experimental data and calculations. 





o 
° 


RELATIVE PERMEABILITY — & 
> 
° 


nN 
io} 








UNCONSOLIDATED SAND (BOTSET) 
K2I72¢ 82 035 

— CALCULATED LIMITING DATA 

Sse us 

















80 


40 60 
LIQUID SATURATION — & 


FIG. 12 


PETROLEUM TRANSACTIONS, AIME 


2. Kozeny, J.: 


. Schriever. Wm. 
. Traxler and Baum: Physics, (1936) 7, 9-14. 
. Carman, P. 


RELATIVE PERMEABILITY TO LIQUID IN LIQUID-GAS SYSTEMS 


NOMENCLATURE 


(Symbols by Order of Appearance } 
Usual Subscripts: LZ refers to liquid 
G refers to gas 
/ refers to interfaces 
S refers to saturation 
Volumetric Rate of Flow per Unit Cross Sectional 
(rea 
Specific or Effective Permeabilities 
Relative Permeability 
Pressures 
Viscosity 
Porosity 
lextural constant 
Total Boundary Areas (in cm’ per cm* of bulk 
volume) 
Interfacial Surface Area Between Liquid and Gas 
(in cm* per cm’ of bulk volume) 
Contact Area Between Liquid and Solid (in em’ 
per cm of bulk volume) 
Contact Area Between Gas and Solid (in cm’ per 
cm‘ of bulk volume) 
Capillary Pressure 
Liquid Saturation 
Irreducible Minimum Liquid Saturation as Ob- 
tained from Capillary Pressure Curve 
Reduced Liquid Saturation 
Areas Under the Capillary Pressure Curve Times 
the Quantity — (f//v) 
Interfacial Tension 
Free Surface Energy 
Contact Angle of Liquid 
ence of Gas 


Between Liquid and Gas 


Against Solid in pres- 
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APPENDIX A 
THERMODYNAMICAL ASPECTS OF DISPLACE- 
MENTS IN POROUS MEDIA 


The capillary pressure curve obtained by displacing the 
wetting phase from an initially fully saturated porous medium 
may be defined as representing an isothermal equilibrium 
process, or more exactly, a succession of equilibria. The pres 
sure on the non-wetting phase (gas) is increased in a step- 
wise fashion resulting in the gradual displacement of the wet 
ting phase (liquid). As this equilibrium process takes place 
it can be said that (a) the interfacial area between phases 
is varying in a continuous fashion, probably increasing until 
the irreducible minimum liquid saturation is reached, and 
(b) the solid/liquid contact area diminishes and is replaced 
by a solid/gas, or more exactly by a solid/liquid-film/gas 
area. For the total system formed by the liquid and the solid. 
the displacement process represents an_isothermal-isobari« 
transformation. During such a transformation it is known 
that the external work done against the system (by the non 
wetting gas phase) is equal to the increase of free energy of 
this system. The capillary pressure curve can thus be inter 
preted as representing the change of free surface energy that 
takes place in the porous medium during the displacement 
process. 

Referring to the definition of surface areas as given by 
Equation (7), the change in free surface energy due to an 
increase of the interfacial area can be stated as: 

dF, =v dl F ae” (A-1) 
If during the process of desaturation, solid area previously 
covered by the liquid phase is covered by gas, the increase 
of the “solid-gas” contact area is dR, and the corresponding 
change in free surface energy can be stated as 
@.=@ea tee . «ss (A) 
Addition of Equations (A-1) and (A-2) obtains: 
dF = dF, + dF, =v (dl + dR, cos @) . (A-3) 
where 7 represents the interfacial tension between the two 
fluids, and @ the contact angle. But dF is the total change in 
free (surface) energy corresponding to a change in saturation 
dS, and is, consequently, equal to the work Pdv done against 
the system, or if a core of unit bulk volume is considered: 
dF = —P..f.dS = y (dl + dRz cos 6) (A-4) 
In the case where the contact angle @ is zero, cos 6 = 1, and 
P..f.dS = -—v (dl + dR) (A-5) 
or 
P. f.dS = - yd Ay, 
and 


S 
Ag =- f J P, ds . 
7 


In explanation of the above equations, if, during the dis 
placement process, equilibrium is not attained. or if an un 


(A-6) 


(A-7) 


*This relation implies that the receding liquid leaves a fi'm over the solid 
surface. It is, however, understood that the liquid particles of this molec 
ular film have no actual bearing upon the distribution and the flow 
mechanism of the fluids 
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stable film is left on the rock surface, the work of displace- 
ment may have a considerably different value than that indi- 
cated in the general Equation (A-4). Equilibrium is obtained, 
however, if the contact angle @ is zero, which is usually the 
case in liquid gas systems. The exact value of Ag can be 
derived only if @= 0, since up to the present no relationship 
exists which would permit the separate identification of d/ 
and dRg. 

Under the above circumstances, the areas under the capil- 
lary pressure curve can be assigned the following meanings 
(see Fig. 3) 4 

Pe 
/ 
Ze : J P. dS represents the interfacial area at 
7 
the saturation S together with the rock area exposed to the 
non-wetting phase, less the amount, /,, of interfaces that ex- 
isted at the bulk surface of the fully saturated core; /, is 
comparatively very small and can be neglected so that 


S 
j E 
: J P. dS =1 + Ro-Iy™~1 + Re = Ae . (AB) 
The maximum limiting value of the above integral can be 


represented by 
s 
j F 
= J P ds = Ry LIM 


In this expression J.) 1. 


T Tsay (A-9) 


represents the interfacial area be- 
tween the gas and the “connate” liquid, and consequently 
Reas, the (maximum) gas/solid contact area that obtains 
at the irreducible minimum liquid saturation. On the basis 
of considering the connate liquid as being immobile, it is 
seen that 2s, which corresponds to the envelope of the con- 
nate liquid and solid area. represents the surface area A, of 
the total effective flow matrix.* so that 


(A-10) 


S 
/ 
| P.dS 1, R, t R, 


*Such an interpretation is equiva'ent to treating /,,,, as part of the total 
solid” area, which always remains covered by the liquid phase and is 
at any saturation 


therefore, automatically included in R 
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Is is then finally seen that 


S 


f P.dS = Ag-R,-1 
S 


/ 
¥ ¥ 7. — ~ 


or, combining Equations (A-10) and (A-11) 
te =R,-1=A,-21 
The above derived relations (A-8), 
used for the derivation of K, (Max) and K, (Min). 


APPENDIX B 
CALCULATION OF SPECIFIC PERMEABILI- 
TIES FROM CAPILLARY PRESSURE DATA 


According to the derivations presented in the Main Report 
and in Appendix A, the specific permeability of a porous 
medium can be expressed as 


1 (1-S,)’ f 
k= a ae (B-1) 
5 (Ay) 


or using Equation (A 


(B-2) 


10) 
S 

P. ds (1-S,) 
] 


imtly 


where f, S and S, are fractional: P. is expressed in dynes per 
sq cm, ¥ in dynes per cm. and & in sq em. 

The above relations imply that not all of the fluid filling 
a porous body participates in the flow. and that the fraction 
held immobile is approximately equal to the irreducible mini 
mum saturation as indicated by capillary displacement experi- 
ments.* Equation (B-2) may be used to calculate the specific 
permeability from cap‘llary pressure data. 

Table I shows the comparison between calculated and meas 
ured dry air permeability for 27 samples which varv widely 
in permeability. porosity, and textural characteristics. It might 
be noted that in 21 cases the calculated permeability is smaller 
than that actually measured. Such a systematic trend may be 
explained by the fact that the measured dry air permeabili 
ties have not been corrected for gas slippage effects. 

It is seen that the average discrepancy between the meas- 
ured and the calculated results is greater for the natural cores 
than for the artificial cores. It must, however, be indicated 
that the degree of precision of the experimental data on most 
of the natural core samples is not very high. The core material 
itself seemed to be extremely sensitive to chemical agents so 
as to be subject to definite changes (in surface and even tex 
tural characteristics) with time. 

In the light of these remarks the overall agreement between 
calculated and measured values of permeability can be con- 
sidered as satisfactory and may be interpreted as a good veri 
fication of the presented theoretical approach. 


APPENDIX C 
DERIVATION OF THE LIMITING K, MAXIMUM 


Under any circumstances. it can be considered that the gas 
inside the porous medium is flowing through a certain chan- 
nel system surrounding clusters of rock grains containing 
liquid in their interstices. The assembly of all these clusters 


*It has already been recognized by King,” and later by Kozeny,"* as we'l 
as by other workers,” that in a single phase system some part of the 
fluid seems to be held back in internal angles and fissures, and in zones 
of grain contacts where molecular forces of adhesion become prevalent 
Therefore, in the original Kozeny equation.’ the porosity, f, should be 
assigned a slightly lower value, and the surface area, A, should actually 
correspond to some “effective flow matrix” and not to the specific solid 
surface. 
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may be treated like an entirely water saturated “reduced” 
porous medium, the porosity of which would be equal to 
Sf 
1-f + Sf 
sv that under the assumption of isotropic random packing. 
a total cross sectional liquid area 


a total cross sectional solid area 
S} water volume 


1—f 
\ schematic illustration is given in Fig. 14, where one pore 
opening between two rock grains is represented. 
The total outside surface area of the “reduced” porous 
medium is equal to 4, = 1 + Ry. The following interpretation 
can then always be given to the ratio 


I 
R, 


= ne (C-1) 
rock volume 


namely, 
/ liquid surface exposed to gas : 
= — ——— (C-2) 
R solid surface exposed to gas 
Under the previous assumption of isotropic packing, the in- 
equality 
/ a, S} va 
>—= ee 
R, a. l-j 
can be established with help of the following reasoning. The 
liquid is the wetting phase so that it actually ‘ 
the rock, creating a total liquid-gas contact area, /, that is 
greater than the (average) cross sectional area, @,, of the pore 
space containing the water. As a consequence, the total solid 
surface that remains exposed to the gas, R., will be smaller 
than the (average) cross sectional area. a.. of the solid grains 


“spreads” over 


so that one has simultaneously: 

CJ =i>'4, P (C-3") 

CMeSree@: < sli; (C-3") 
where i and r represent the projections of the curved surfaces 
/ and R 
C is the factor of projection or shape factor of these surfaces, 
and is assumed to be the same for the rock grains and the 
interfaces. By combination of (C-3") and (C-3"), inequality 
(C-3) is obtained. In fact. at low liquid saturations C is 
smaller for the interfaces than for the rock which still in- 
creases the resultant inequality (C-3). 


rhe 


respectively on the considered cross-sectional plane. 


ratio 
i 
a; 
can be approximated in the following way: the surface 7 can 





PROJECTIONS ON CROSS SECTION 


FIG. 14— SCHEMATIC COMPARISON BETWEEN ROCK SURFACES AND 
INTERFACES 
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be considered as corresponding to a mean hydraulic radius 
given by 


(C-4) 


whereas the mean hydraulic radius of the pore space, Le., 
that corresponding to the cross sectional area. a, can be shown 
to be given by 
27 
A, = aS arin Sieh ae ar’ ed 
P.. 
where /,, represents the mean capillary pressure over the 
saturation range S to S,, Le., 


P,= vas fs new Ee 


5. -S, 
The ratio of " surtaces i ba a, can then be expressed as 


ieee 


Combining Equation (C-7) with the relations (C-3), (C-3") 
and (C-3') gives 


l S/ -. ) 

Ry 1-f a 
ees 
The combination of (C-9) and a relation 


Tra 


trom whence 


/ minimum = R, - 


(See Appendix A 


Equation A-8) 


aes) 


Substitution of the above expression in the equation 


De 
f P.dS + 21 
S 


/ Minimum = A, - C-10) 


(See Appendix A 
Equation A-12) 
gives 


S 


Ss. 2f 
I. P.dS + 


4, Minimum = 


aT rd 
Sf - 
In accordance with Equations aly of the Main Nie pew 
(A-10) of Appendix A, the final expression for the maximum 
of the relative permeability is then given by 

ie \* 

S* | 
| f Pads | 
U4) J 


K, Maximum = es 


| 
S 2s, P.dS 


frase ri(RY 


or using the symbols of the Main Report and Appendix A 


K, Maximum at - : ay (C-13) 


“i(F y| 
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Member 
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By Charles G. 
Bartlesville, Okla., 


The authors’ work is a most stimulating contribution to the 
study of the complex problems involved in the flow of fluids 
through porous media. This reviewer has profited greatly from 
a study of the paper which undoubtedly represents the end 
careful analysis. The following comments 
are presented in an effort to implement rather than detract 
from the value of the ideas presented. 


product of much 


(Among the fundamental assumptions, stated and implied, 
upon which the paper is based are the following: (1) In 
single-phase fluid flow through a porous medium the fluid 
filling all minute crevices, cracks and pores is stagnant and 
does not participate in flow. (2) The total volume of stagnant 
liquid may be identified with the irreducible minimum liquid 
saturation measured by a capillary pressure experiment. (3) 
f a flowing liquid is under consideration, the volume of each 
pore which is filled with stagnant liquid is determined by the 
e., by the solid geom- 
etry and the surface tension of the liquid. (4) There is no 
slippage at the boundary between flowing liquid and stagnant 
liquid. (5) In relative permeability to liquid measurements 
when the gas phase is immobile there is no slippage at the 
liquid-gas interfaces within the voids of the porous medium. (6) 
In analyzing the free energy change involved in the displace- 
ment of liquid from a porous medium by gas it is assumed 
that a solid-liquid interface is replaced by a solid surface cov- 
ered by a liquid film having the surface tension of the bulk 
liquid or by a surface of equal energy. It would be profitable 
assumptions in turn. 


curvature of the liquid-gas interface, i 


to discuss these 


The Kozeny equation’ has been applied by Carman’** and 
others to determination of the specific surface areas of finely 
divided materials by measuring the fluid permeabilities of 
unconsolidated beds. The method has enjoyed wide acceptance 
in recent years, and the surface areas measured thereby agree 
well with area measurements made by independent methods. 
although it is recognized’ that the internal areas of porous 
silica gels, ete., are not evaluated by the method. For perme- 
ability experiments to measure a solid surface area, the solid 
surface must offer a frictional resistance to flow, and the flow- 
ing fluid must everywhere contact the surface. In recent work 
at the Bureau of Mines laboratories’ we have obtained agree- 
ment to better than three per cent between the specific surface 
areas of a sample of 70- to 100-mesh glass beads as calculated 
from liquid permeability experiments to water and isoéctane. 
and independent miscroscopic measurements. These 
glass beads were of the average size of grains in many sand- 
stone rocks. Packed beds of the beads probably 
would retain 10 to 15 per cent water at completion of a dis- 
placement experiment as an irreducible minimum water satu- 
ration, but no correction was made for this in our work such 
as the authors suggest in their Equation (B-1). On the basis 
of the above remarks, it is felt that assumptions (1) and (2) 
do not have general applicability to all porous systems. How- 
the sandstone reservoir rock material, for 
example, the presence of cementing material, clay minerals, 
and very fine silt between coarser sand grains alters the pic- 
ture. Undoubtedly the fluid in not 
participate in flow. 


precise, 


reservoir 


ever, in case of 


such fine crevices does 


In the case of assumption (3) it is difficult to understand 
how the curvature of a hypothetical liquid-gas interface can 
determine the boundary between flowing and stagnant liquid 
in a continuous liquid phase. It would seem more appropriate 


to attribute such a phenomenon to attraction of the solid for 
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the stagnant liquid, but no analysis of this situation has been 
developed. On the other hand the concept serves as a working 
hypothesis to stimulate further research. 

Assumption (4) is a difficult one to accept, but it may sit 
better on further reflection. Assumption (5), on the other 
hand, seems to be untenable. It is involved in the authors’ 
definition of A, which is used throughout their development. 
This reviewer would suggest that the analysis be reconsid- 
ered on the basis of a more plausible definition of A,, elimi- 
nating the gas-liquid interface as one where friction is 
involved. 

In Appendix A the authors, apparently without realizing it, 
have made assumption (6), which seems a reasonable one to 
this reviewer if the saturating liquid wets the solid with a zero 
contact angle. The authors assume the work done in the dis- 
placement process is that necessary to form new liquid surface 
having an area Ag. It should be pointed out that this results 
in Rg always being zero, a condition the authors use in evalu- 
ating their minimum K, function. However, if Rq is always 
zero, the ratio //R, used in evaluating the maximum K,, func- 
tion loses its significance. 
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AUTHOR’S REPLY TO MR. DODD 


Dodd’s very clear analysis of the assumptions used in the 
theoretical developments of this paper contributes toward the 
recognition of some fundamental aspects of fluid mechanics in 
porous media and represents a most valuable basis for con- 
structive discussions. The authors want to take this opportu- 
nity for indicating that one of their main purposes has been 
to attempt the definition of a series of “working hypotheses” 
and to outline thereby a logical basis for further research. 

With respect to assumption (5) the authors can hardly 
follow the views expressed by Dodd. It is not excluded that 
there may be some slippage at the liquid-gas interfaces. Ac- 
cordingly, a more elaborate derivation could possibly be 
worked out, provided there were any tangible indication as to 
the exact amount of slippage to be accounted for. It can be 
anticipated that the consideration of slippage on part of the 
liquid boundaries would lead to a situation that is somewhat 
different from what is usually admitted, namely applicability 
of Darcy’s law to each of the flowing phases. It is believed 
that in any event neither the principle of the theoretical treat- 
ment nor the nature of the boundaries would be altered. At 
this time, therefore, it does not appear that the definition of 
the liquid boundary 4, ought to be modified or even that the 
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assumption of zero slippage at the liquid-gas interfaces is 
untenable. 

It is felt that some confusion might arise in connection 
with the comments concerning assmption (6). As Dodd points 
out, there is actually no formation of true molecular contact 
between the gas and the solid surface since it is implied that 
the latter remains at all times covered by a liquid film. This 
film must be visualized merely as an immobile coating layer 
of infinitesimal thickness. With respect to the spatial distribu- 
tion of fluids, it is necessary to distinguish between the inter- 
facial area, /, separating the gas from the flowing bulk liquid, 
and the area of triple contact “gas/liquid - film/solid” which 
for reasons of simplicity has been termed “solid-gas” contact 
area, R,. From the foregoing it is then clear that R, must 
always be assigned a definite value and that it becomes zero 
only in the case of complete gas channeling such as may be 
seen from Figs. 1 and 4 


DISCUSSION 


By M. R. J. Wyllie, Gulf Research and Development Co., 
Pittsburgh, Pa.. Member AIME 


The authors are to be congratulated on a most interesting 
and provocative paper which will, I am sure, give rise to much 
discussion. For this reason it is my intention to discuss here 
only those aspects of the authors’ work which appear to be 
fundamentally at variance with certain published works of 
Rose and myself.’* My hope is either to reconcile the dif- 
ferences or failing that, at least to underline them so that some 
later worker may be in a better position ultimately to resolve 
them. 

Firstly, | am most happy to see that more work applying 
the Kozeny equation to consolidated porous media is now in 
print; nevertheless, I feel that a word of caution is in order 
since the application of the Kozeny equation to such media 
has not yet received wide support; indeed, Carman, one of 
the principal protagonists of the theory, has outlined at length 
the rather severe theoretical requirements which must be met 
before its use is justified.’ Nevertheless, I agree that the valid- 
ity of this equation when applied to consolidated porous 
media can only be tested by experiment and the authors data 
make further such tests possible. 

Where I differ from the authors is (a) in the employment 
of the “effective porosity” concept in the type of porous media 
used by them and (b) in the magnitude of the Kozeny con- 
stant, t, applicable to consolidated porous media. 

With regard to (a) I believe that this concept, as suggested 
by Kozeny,’ Zunker® and Carman’ for clays, is valid; I have 
found it useful myself’ for the theoretical treatment of fluid 
flow through mud filter cakes. Nevertheless, clays are unique 
in that they undoubtedly absorb strongly an appreciable water 
film. Thus, the concept of the effective flow porosity of a clay 
particle system being less than the porosity determined solely 
by the amount of “interstitial” water present in the system is 
not unreasonable; indeed some such postulate is essential if 
the well-known fact is to be rationalized that the permeability 
of a clay particle system may be essentially zero when its 
porosity is appreciable. A similar concept may apply to all 
porous systems, as pointed out by the authors; here the dis- 
cussion by Carman* quoting the results of King,” Erikson” 
and Darapsky,” all of whom tended to a view similar to that 
held by the authors, is relevant. Nevertheless, it is, I believe. 
important to note that all direct substantiations of the Kozeny 
equation, including those cited by the authors, involve the 
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use of the total surface areas and total porosities of the sys 
tems investigated; i.e., in Equation (5) of the authors, k, / 
and A were directly measured and thus t found. It is, of course, 
possible that, by the use of appropriate effective porosities 
and surface areas, the same constants would have been found 
by the workers quoted in the authors’ references;*""" never 
theless, this was not done and hence it does not appear wholly 
appropriate for the authors to cite these data as a justification 
of their choice of a Kozeny constant of 5. 

It may be emphasized that, as noted by the authors, the 
Kozeny constant of 5 is experimentally proved to be only an 
approximate figure for flow through isotropic unconsolidated 
porous media. In non-isotropic unconsolidated porous media 
the constant varies with the direction of flow being, respec- 
tively, approximately 3 and 6, for flow parallel to and per 
pendicular to the long axes of a bed of fibers.” It is the view 
of Sullivan and Hertel* and of Rose and myself that the 
Kozeny constant is a function of what has been called the 
tortuosity of the porous system through which flow occurs. In 
the case of the fiber bed quoted above the tortuosities are ap 
proximately 1 and 2 respectively and, if a reasonably con 
stant shape factor of about 3 is assumed (Carman* adduces 
evidence for this), the directional Kozeny constants of 3 and 6 
are obtained as the product of the tortuosity and shape factor. 
In consolidated porous media, Rose and I’ have noted from 
electrical evidence that tortuosities much higher than those 
applicable to unconsolidated systems seem applicable, and 
hence we assume that Kozeny constants much exceeding 5 
are applicable to such unconsolidated systems. We have ad 


Table I 
NATURAL CORES 
K 1 K,,- K, 
ris Ky, 

‘ md in ¢ 
15.7 30.0 + 90.0 
17.0 24.5 + 44.0 
25.0 22.8 8.8 
26.0 39.5 +29.0 

117.0 86.5 26.6 

127.5 104.0 18.5 

133.0 103.0 22.5 

139.0 108.0 

154.0 112.5 

159.0 131.0 

174.5 149.0 

182.5 199.0 

213.0 150.0 

226.0 158.0 

243.0 192.0 

301.0 244.0 

470.0 271.0 

21.5 660.0 590.0 

18.3 2,543.0 2,710.0 

35.0 17,200.0 12,000.0 


Average Error on Natural Cores: 26.0 Per Cent 


Designation ye 
B14 9.3 
B 3 15.4 
21.2 
11.3 
14.2 
13.6 
16.3 
15.8 
14.8 
16.2 
12.6 
15.3 
16.6 
17.0 
17.3 
16.8 
17.8 


Formation 
Third Deese Sand 
Third Deese Sand 
Bartlesville Sand 
Third Deese Sand 
Third Deese Sand 
Third Deese Sand 
Third Deese Sand 
Third Deese Sand 
Third Deese Sand 
Third Deese Sand 
Third Deese Sand 
Third Deese Sand 
Third Deese Sand 
Third Deese Sand 
Third Deese Sand 
Third Deese Sand 
Third Deese Sand 
Nichols Buff Sand 
Bromide Sand 
Unconsolidated Sand 


B13 
B10 
B25 
B5 
B37 
327 
B33 
B19 
B55 
B41 
B 1 
B48 
B 8 
B4 


22.0 
27.0 
17.5 
14.5 

9.0 
29.5 
32.0 
21.0 
18.6 
42.5 
10.6 

6.6 
30.0 


Average Error: 


f% 
Kx 
Ko = 
Ko meal Ky 
Ky 
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duced some experimental evidence to support this contention. 
In addition we have noted’ that tortuosity appears to vary 
with saturation in all porous media and give as the relation- 
ship between 7, the tortuosity at 100 per cent saturation, and 
T., the tortuosity at some saturation, S, of wetting phase, the 
relationship, 
t./7 ST 

where / is the electrical resistivity index pertaining to the 
saturation, S. It may be noted that a variation of tortuosity 
with saturation has recently been suggested also by Fatt and 
Dykstra.” 

It is important then to point out that Rapoport and Leas 
use a constant of 5 for consolidated porous media and specifi- 
cally state that to a good approximation they believe there is 
no variation of the constant with saturation. Hence their posi- 
tion is quite different from that of Rose and myself and 
clearly the divergence in outlook can only be settled by refer- 
ence to the experimental facts. Table I in the authors’ paper 
gives experimental results which the authors cite as a general 
confirmation of their theoretical reasoning. They note that 
although the agreement between measured and computed per- 
meabilities is not precise, there are considerable experimental 
difficulties attached to work of the kind carried out and thus 
the over-all average agreement of +12 per cent between com- 
puted and measured permeabilities supports the validity of 
the Kozeny constant of 5 which they employed in all cases. 
Certainly I do not question the experimental difficulties in- 
volved in compiling Table I and the overall agreement noted 
must, therefore, be considered as rather good support of the 

-_ 


Comparison Between Calculated and Measured Specific Permeabilities 


ARTIFICIAL CORES 


Core 


Materia! 


Designation I% 


Alundum 
Alundum 


Alundum 
Alundum 


27.4 
30.0 
33.2 


567.0 
1,385.0 
2,710.0 


566.0 
982.0 
2,700.0 


Alundum 
Alundum 
Penn State 
Average Error on Artificial Cores: 13.1 Per Cent 
22.6 Per Cent 

Porosity 

Measured Permeability in md 

Calculated Permeability in md 


Error in Per Cent 
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authors’ contentions; alternatively the data may be indicative 
of the fact that there are at least two plausible methods of 
skinning this particular cat. 
The computed permeabilities in Table I are derived essen- 
tially from the relationship: 
(1-S.)° f 


1’, 


(1) 


(2) 


with ¢ > 5 in all cases. (Actually ¢ = F’f't, where F = forma- 
tion factor and ft, = shape factor and is approximately 
2.5 — 3.0.) 

Now in all cases (1 
1, < A. (It seems pertinent to note that the use by the authors 


S.)'°f<f and it seems probable that 


S, : ‘ 

of —f/y { “P.dS = Ay, is not wholly in accordance with the 
l 

recommendations of Leverett in Reference’ of the authors.) 


Depending upon the relative values of (1—S,)’f. Ay and ¢ 


it would appear that both Equations (1) and (2) above could 
give approximately correct values of k. At this stage I do not 
believe it possible to decide which of the two equations is the 


more fundamental. 
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AUTHOR'S REPLY TO MR. WYLLIE 


The authors fully appreciate the importance of Wyllie’s 
comments. They feel, however, that a sufficiently detailed an- 
swer would require considerable divergence from the issue at 
hand. In the near future it is planned that a note containing 
additional experimental data will be written concerning the 
significance of the textural constant. 
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DISCUSSION 
By J. S. Osoba, Humble Oil and Refining Co., Houston, Tex. 


The authors have presented a useful method of calculating 
relative permeability to the wetting phase under drainage con- 
ditions when the capillary pressure versus saturation relations 
are known for reservoir rock. The agreement between calcu- 
lated and experimental relative permeabilities is proof of the 
validity of the calculated method provided the experimental 
method is free of errors. 


Several methods of measuring relative permeability have 
been investigated in our labortory. The relative permeability 
to the wetting phase measured by three of these methods, the 
Penn State, single core dynamic, and gas drive technique, were 
nearly identical to the curves calculated from the capillary 
pressure using Rapoport and Leas’ equation for 
K, (min). One method, the Hassler method, which is most 
nearly like the method used by Rapoport and Leas, gave a 
relative permeability to the wetting phase that was lower than 
the calculated curve. using the equation for K, (min). 


curves 


Since there are laboratory methods whereby the relative per- 
meability to the wetting phase can be measured and the results 
coincide with the calculated relative permeability to the wet- 
ting phase, and since these methods are many times faster 
and less difficult to operate than the method presented in this 
paper, I do not believe the authors are justified in concluding 
that their method of calculating relative permeability is supe- 
rior to experimental methods. 

With reference to another point in the paper, in calculating 
specific permeability by the Kozeny equation [Equation (6) }. 
the result is very sensitive to the value chosen for the irredu- 
cible minimum saturation on the capillary pressure curve. This 
is due to the fact that the irreducible minimum saturation 
determines the factor A [in Equation (6) ]. and a small varia- 
tion in this irreducible minimum causes a large change in A. 
Since this factor appears in the equation as A’, the accurate 
determination of this factor is very important. It would be 
of some help in evaluating the data presented in Table I if 
the procedure used to determine the factor A were explained. 


AUTHOR'S REPLY TO MR. OSOBA 


Osoba’s statement concerning the agreement that he found 
between and direct 
liquid relative permeability is most gratifying and provides a 
further substantiation of the theory advanced by the authors. 


the calculated curves measurements of 


It is readily admitted that there is a certain resemblance 
between the Hassler method and the experimental technique 
presented in this paper insofar as both require the use of 
semipermeable capillary barriers. It may, however, be stated 
that the results obtained by the Hassler method are subject to 
due to the distortion of the flow pattern 
from linear whereas the procedure used by the authors has 


a systematic error 


been conceived so as to eliminate any such anomalies. 


As regards the superiority. from the viewpoints of cost and 
time, of experimental methods over that of calculating. the 
latter has been recommended merely in preference to that 
experimental method which is presented in this paper. The 
authors certainly do not deny the possibility of more econom- 
ical methods. They would, on the contrary, advocate the use 
of more advantageous methods, as soon as their reliability and 
economy could be sufficiently demonstrated. x * * 
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The pressing duties of arranging a 
successful AIME Annual Meeting in 
St. Louis, Mo., last month prohibited 
Robie’s writing his column for this is- 
sue. With apologies to the patrons of 
this page, Jess E. Adkins, editor of the 
JOURNAL OF PETROLEUM TECHNOLOGY, 
is substituting for Robie. 


Annual Meeting 


The nearly 2,500 members who at- 
tended the Annual Meeting last month 
in St. Louis came away from that city 
with the feeling that those who did not 
attend missed a good thing. The 17Ist 
General Meeting of the AIME was in- 
deed a success. The amount of effort 
expended in preparing some 400 papers 
for presentation is immeasurable, but 
evident in the manner in which tech- 
nical sessions of all divisions and 
branches were attended. 

Dr. Donald H. McLaughlin, immedi- 
ate past-president of the Institute, was 
moved to remark at the Welcoming 
Luncheon that the last time the AIME 
met in general session in St. Louis was 
in 1884. The records of that meeting 
are not at hand, but it is wondered if 
the city in those days could have 
handled as gracefully the many simul- 
taneous meetings which took place in 
1951. No less than four hotels were 
involved in this most recent gathering. 
The Jefferson Hotel was almost taken 
over by the Institute while the Statler 
Hotel was predominantly occupied by 
the Petroleum and Metals branches. 

The ten petroleum papers and the 
production review session were all well 
received, evoking considerable discus- 
sion from the floor. Several of the pa- 
pers will become a part of the Trans- 
actions while a summary of domestic 
production for 1950 is carried else- 
where in this issue of PETROLEUM 
TecHNoLocy. A similar summary of 
foreign production for 1950 will appear 
in the April issue. 

At the Welcoming Luncheon, the In- 
stitute heard the eminent physicist. 
Arthur Holly Compton, chancellor of 
Washington University, warn the world 
against another earth-wide conflagra- 
tion. He expressed his belief, though 
not in an optimistic tone, that by the 
very nature of the destructiveness of 
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. . as followed by EDWARD H. ROBIE 


Secretary AIME 








the atomic bomb (and subsequent de 
velopments along the same lines) such 
a catastrophe would be avoided at least 
for a time. He also defended, from the 
physicist’s standpoint, the development 
of such a destructive force. The role of 
the scientist is not, Dr. 
pointed out, to bestow happiness upon 


Compton 


a people; rather it is to bestow great 
ness, 

Texas Railroad Commissioner Wil 
liam J. Murray made an urgent plea 
for the submission of accurate reserve 
estimates by petroleum engineers. As 
guest speaker at the Petroleum Branch 
banquet, Murray was concerned about 
the apparent chasm between the painted 
picture of petroleum reserves and pro- 
duction possibilities and the actual pic- 
ture. He reported that allowables have 
jumped in recent months due to mili- 
tary needs but that “pre-Korea” esti- 
mates of production capacities have 
fallen short of their claims. He posed 
these questions: What amounts should 
import volume reach? What allowables 
should be set? To what extent should 
we depend upon Middle East produc- 
tion? The answers, Murray observed, 
could be found only through the scru 
tiny of more accurate engineering re 
ports of estimated reserves 

W. M. Peirce. AIME president, at- 
tended the Petroleum Branch banquet 
and urged that the Institute’s goal for 
1951 be a substantial increase in mem- 
bership. 


Random Notes 


There was no deficiency of “favors” 
at the various functions in St 
At the Welcoming Luncheon members 
were given palladium tie clips while 


Louis. 


the ladies received ear clips of the 
same metal. At the Dinner-Smoker 
clothes brushes of aluminum (except 
for the bristles) were donated and 
leather bound desk sets were given out 
at the Annual Banquet 
offered a sight seldom seen anywhere 
except in the United States. Tickets 
were marked with the words “Dress 


. The banquet 


optional.” As it turned out, the manner 
in which a man was dressed could al 
most be matched with his table location 
in the room. White ties were popular 
at the speakers’ dais while the front 
center section was dotted with dinnet 
jackets. Those in the rear were prob 
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ably the most comfortable of all in con- 
servative street clothes. The ladies, of 
course, were formal to a man, as it 
No such problems of choice 
existed at the Petroleum Branch Ban- 
quet which was informal. 


were. 


Technical sessions of the Branch were 
all held in the St. Louis Room of the 
Statler Hotel. An overflow crowd was 
the rule but no one turned away be- 
cause of a lack of seating facilities. . . 
They merely chose a good location on 
the floor and listened Indian style. 

\ word of appreciation should be ex- 
tended the projectionist on duty during 
the petroleum sessions. He knew what 
he was doing every minute. . . The 
Statler lived up to its advertising. Cour- 
tesy was the rule, and service (in prep- 
aration of meetings) was carried out 
unobtrusively. A microphone was pro- 
vided on five minutes notice which is a 
. The Dinner- 
Smoker might also have established a 
record. Five male voices stole the show. 

.. The maturity of St. Louis was evi- 
dent everywhere. Some, who might view 
business district traffic signals as a 
necessity, might even think it quaint. 
Traffic policemen seemed to have things 
under control at all times, however. At 
least, no pedestrians were run down 
during the four-day meeting. 


record of some sort. 


(merican ingenuity was manifested 
in the acts of some members who waited 
until the last minute before flying to 
Weather conditions 
throughout the Midwest had planes on 
the ground over the week end of Feb. 
17-18. Stranded in Springfield, Mo.. 
Bob Gilmore of DeGolyer and Mac- 
Naughton rallied several plane loads of 
folks and chartered an earthbound bus. 
They all arrived in St. Louis heavy-eyed 
but present at 4 a.m. . . Bob had a 
Petroleum Branch treasurer's report to 
make. 

All St. Louis cabs bore “Welcome 
AIME” signs on windows. So did 
.. The AIME Section 


Delegates report to the Board of Direc- 


the meeting. 


some of the bars. 
tors received many compliments for 
constructiveness. Some action was taken 
. Alaska has a new local 
. . Economy was the keynote at 
the Board meeting. Measures requiring 


on all counts. . 


section 


expenditures received heavy considera- 
tion. * = 9 
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LPG C 
Injected Recovered 
Ga Gal 


63,452 60,867 
62,997 62,801 
13,487 13,416 
13,075 12,951 


LPG 
Recovered 
% In 
95.9 205 
99.7 
99.5 
99.1 


LPG Storage Cavities 


Continued from Page 11, Section | 


example, a second string of tubing with a bottom hole pump 
may be employed to remove salt water from the reservoir as 
LPG is injected into the annular space, thus avoiding exces 
sive pressures on such formations. Fig. 4 shows the necessary 
formation depth to permit a desired surface pressure with 
propane storage. 

Also in formations that will not support excessive gradients, 
LPG may be injected through the tubing under gas at sufficient 
pressure to maintain the desired surface pressure. Storage in 
contact with gas, however, will alter the vapor pressure of 
the stored product which may be objectional in some cases 


TESTING 


The Upton County project has been tested only once. On 
this test, 1,500 bbl of propane were injected, of which ap 
proximately 96 per cent was recovered. This was deemed ade 
quate proof by the operator of the feasibility of the method 
and the reservoir was placed in commercial service. Enlarge 
ment of the Winkler County storage system was shut down at 
approximately 8,000 bbl capacity and four tests were con 
ducted to determine the practicality of the method, and to 
allow public demonstrations of the process. The results of 
testing are shown in Table I. 

In all tests LPG was injected and withdrawn at a rate of 
approximately 300 bbl per hour, the limiting factor apparently 
being only the size of the tubing. The effluent LPG was found 
to have increased in vapor pressure, which is believed to be 
primarily due to absorption of residue gas used over the LPG 
in the surface tanks to increase the pressure on the suction 
end of the injection pump. Also, cobalt bromide tests on the 
recovered LPG indicated a very slight trace of moisture. The 
primary question to be answered by this testing was what 
proportion of injected products could be recovered. It is be 
lieved that the moisture content resulted from the withdrawal 
of all recoverable products in each test, and that under normal 
operation, withdrawn products probably will be dry. 


APPLICABILITY 


To date, no attempts have been made to apply this prin 
ciple to formations other than salt, but as many shales readily 
disintegrate in water, it is believed that in many localities 
shale beds exist suitable to the process. Also in some cases 
it may be feasible to dissolve space out of massive limestone 
with hydrochloric acid. 

Effort is being made to locate all salt deposits suitable to 
ihe process, but the study is by no means complete. The Per 
mian salt of West Texas extends over into Eastern New Mex 
ico, up through West Central Oklahoma and Central Kansas 
almost to the Nebraska line. A sizable deposit exists in 
Western Colorado and Eastern Utah and some small deposits 


March, 1951 


Vapor psi 
@ 105°F 


JOURNAL OF PETROLEUM TECHNOLOGY 


Recovery Tests of Winkler County Storage System 


Corrosion Cobalt 
F Test Test 
Out In Out In Out 


225 No No Dry 
No No Dry 
No No Dry 
No No Dry 


Trace 
Trace 
Trace 
Trace 


are in Southern California. Also extensive salt beds underlie 
most of Michigan and extend into Eastern Ohio, Northern 
West Virginia, Western and North Central Pennsylvania, and 
Southwest New York. The salt domes of the Gulf Coast are 
believed applicable to the process. 

In consideration of the probable, though untried, applic- 
ability to shales and lime stones and proven economy of the 
system in salt beds, it is believed the formation of underground 
storage space by artificial means will alleviate a considerable 
portion of the LPG storage problem. 
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With Geolograph, ‘you 
actually watch the dritt- 
ing change as the bit 
bites into contrasting 
strata — foot by foot! 
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Personals 





W. T. Nicutincace of Rock Springs. 
Wyo., has been elected president of 
the Mountain Fuel 
Supply Co. of Salt \ 
Lake City, Utah. 

He has formerly 

been with White- 

hall Petroleum 

Corp., Ltd., the 

New York Oil Co., 

and the Ohio Oi) 

Co. In 1929, he 

was made chief ge- 

ologist for Mountain Fuel Supply and 
in 1941 he became a vice-president. In 
1945, he was made vice-president in 
charge of exploration, production and 
transmission. He will move soon to Salt 
Lake City from Rock Springs. Wyo. 


oe 


FrRaNK LinpeMaN, Jr., has been 
named general superintendent of the 
producing department of Stanolind Oil 
and Gas Co. in Tulsa, Okla. Prior to 
his present position, he was in charge 
of Stanolind’s marine operations in the 


Gulf of Mexico. 


+ 


Wituiam B. BLeaKLey is now with 
the University of Tulsa as assistant pro- 
fessor in petroleum engineering. He 
was formerly with Stanolind Oil and 
Gas Co.. Alvin. Tex. 


+ 


J. A. Cricu Ton is now vice-president 
in charge of operations for the San 
Juan Oil Co. with 
headquarters in 
Dallas, Tex. He 
was formerly as- 
sociated with the 
firm of DeGolyer 
and MacNaugh- 
ton, Dallas, as 
vice-president and 
director. With the 
San Juan Oil Co., 
he will be in charge of the company’s 
operations in leasing, drilling, develop- 
ment and production. Presently, he is 
chairman of the North Texas Section 
of the Petroleum Branch. AIME 
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Wituiam E. Bernuarp resigned from 
the Seaboard Oil Co. of Delaware. Cas- 
per, Wyo 


— 


OrMAN PEARSON Brown resigned 
from the Shell Oil Co., Denver, Colo.. 
after nine and one-half years as geolo- 
gist in the Midcontinent and Rocky 
Mountain areas. He has accepted a 
position with the Wilshire Oil Co., Inc., 
Los Angeles, Calif... to open Rocky 
Mountain Division offices in Denver 


A 


Ricuarp C. Harris is now resident 
manager, Venezuelan Atlantic Refining 
Co., Caracas, Venezuela. He was for- 
merly located in Philadelphia for the 
same firm 


+ 


Avpert Lazar has accepted a posi 
tion with the Michigan Consolidated 
Gas Co., Dearborn, Mich. 


+ 


Georce L. Lockett, Mene Grande 
Oil Co., has been transferred from Car 


acas, Venezuela. to Maracaibo. 
af. 


Joun Curtis McCartuy has been 
transferred from the Tulsa office of the 
Stanolind Oil and Gas Co. to the West 
Texas district office in Midland 


— 


GLENN STALEY ts now connected with 
the New Mexico Oil Co.. Hobbs. N 
Mex. 


+ 


Rex W. Woops has resigned his posi 
tion as chief reservoir engineer with 
the Creole Petroleum Corp., Caracas. 
Venezuela, and plans to open a con- 
sulting office in Tulsa, Okla.. specializ 
ing in reservoir engineering 
estimation and evaluation 


reserve 
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MarsHae L. Stirtinc, formerly lo- 
cated in Dowell Incorporated’s Tulsa, 
Okla., office, has 
been made head of 
the firm’s Fort 
Worth, Tex., office. 

He has been with 

Dowell since June, 

1939. He will act 

as sales manager 

in the Fort Worth 

office. Prior to his 

employment at 

Dewell, Stirling was sales manager with 
the Gas Lift Corp. of Houston, Tex. 


— 


Ira J. MEADOR is production foreman 
for the Skelly Oil Co., Tulsa, Okla. 


— 


R. B. Save is now with the Seaboard 
Oil Co. of Delaware, Corpus Christi, 
Tex 


— 


Owen F. THorNTON is now associated 
with The Louisiana Land and Explora- 
tion Co. and will be in charge of the 
Petroleum Engineering Department. For 
the past seven years Thornton has been 
with The Texas Co. as senior petroleum 
engineer where he was concerned with 
reservoir studies, valuations and other 
engineering problems throughout the 
United States, Canada and South Amer- 
ica. He was employed by The Tide 
Water Associated Oil Co. in the East 
Texas area after graduation and until 
1944, 


+ 


Wituiam E. Strives has resigned as 
director of sales for Core Laboratories, 
Inc., to become as- 
sociated with Ken 
Corp.. Long Beach, 

Calif. Stiles, with 
headquarters in 
Dallas, will man- 
age the manufac- 
turing and sales of 
“Ken-Oil” oil-base 
drilling fluid out- 
side of California 
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Proposed for Membership, Petroleum Branch 





Totel A 31, 1950, 
was 17, 082. in sddition: P4534 Soden Associates 
were enroll; 

ADMISSIONS COMMITTEE 

A. J. Phillips, Chairman; George B. Corless, 
John T. Sherman, Lloyd C. Gibson, R. D. Molli- 
son, Plate Malozemoff, Ivan Given. Alternates, 
H. P. Croft, F. W. Hanson, T. D. Jones, Frank A. 
Ayer, H. W. Hitzrot, H. K. Masters. 

Institute members are urged to review this list 
as soon as the issue is received and immediately 
to wire the Secretary's office, ht message 
collect, if objection is offered to “the admission 
of any applicant. Details of the objection should 
follow by air mail. The Institute desires to extend 
its privileges but does not desire to admit persons 
unless they are quolified. 

In the following list C/S means change of 
status; R, reinstatement; M, member; J, Junior 
Member; AM, Associate Member; S, Student Asso- 
ciate; F, Junior Foreign Affiliate. 


CALIFORNIA 

, ee Dietzel, Norman E. (J) (C/S- 
La Bouff, John J., Jr. (J) 
Hamann, Paul M. (J) (C/S- 





Angeles 
S-J). 
San Gabriel -- Christensen, Robert D. 
(C/S-S-J) ; Griffin, Will C., Jr. (J) (c/S8-5). 
Ventura — Bolstad, John H. (M) (R,C/S- 
J-M); Johnson, Bernard R. (J) (C/S-S-J). 


ILLINOIS 


Mattoon 
Olney 
S8-J). 


KANSAS 
Bunker Hill 
S-A). 
English, Homer C. (J) (C/S-S-J). 


Bossler, David P. (J) (C/S-S-J). 
Patrick, Roy E. (M). 


Hammond, Ralph E. (M). 
Loveless, Ernest T., Jr. (J) (C/S- 


Cassell, Robert F. (A) (C/S- 


Wichita 


KENTUCKY 
Prestonsburg 
(C/8-8-J) 


LOUISIANA 


New Orleans Peterson, Albert F., Jr. (J) 
8-S-J); Pyeatt. Woodrow W. (M) (R,C/S- 
S-M). 


MISSISSIPPI 


Laurel — Smith, Francis H. (A) 


MISSOURI 
: Brentwood 


S-J). 


OKLAHOMA 
Bartlesville — Dodd, Charles G. (M). 
Locker, George R. (J) (C/S-S-J). 
-Flesher, Cranston W. (J) (C/S- 


“Healdton Goad, Wilmer E. (J) (C/S-S-J). 

Oklahoma City — Johnson, Grady E. (M). 

Tulsa — Brainerd, Harold W., Jr. (J); Rig- 
sins, William H. (M); Shaw, George R. (J) 


(C/8-8-J). 
TEXAS 


Addicks — 
Beaumont - 
3.M). 


Huzzey, Thomas Edwin (J) 


Remmert, William (J) (C/S- 


Gibson, Robert M. (J) (C/S-S-J). 
Junod, William B. (M) (R,C/S- 


Johnston, William R. (J) (R, 


J). 
mo Spring Miller, Robert N. (J) (C/S-S- 


"lee rger — Nease, Charley W. (J) (C/S-S-J). 

:  eromnane Rodelander, Clovis E. (J) (C/S- 
S-J) 

Corpus Christi Armstrong, John N., Jr. 
£2); _ Charles L. (M); Dornell, Grant 
- Dallas Hudson, Eddie J. (J) ; Stuart, Stan- 

y ). 

ae James M. (J); Thorn- 

’. (M). 

Jackson, James J. (J) (C/S-S-J) 
Levelland — King, Glenn W. (J) (C/S-S-J) 
Midland — Moore, Joe L. (M) (C/S-J-M). 
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Palacios —- Smith, William P. (J) (C/S-S- 
« premont Turner, Willie J. (M) (R,C/S- 


). 
«Randolph Field Kozik, Eugene (J) (C/S 


VIRGINIA 


Arlington White, Foster 
M). 


WEST VIRGINIA 


Charleston Detamore, Waher A. (J) 


WYOMING 
Casper — Churchwell, Robert M. (M) (C/S- 
J-M); Korn, Norman (J) (C/S8-S-J). 
ALASKA 
Juneau 
CANADA 
Alberta Glenn, 
M). 


Glover, Norman J. (J) (C/S-S-J). 
William D. (M) (C/S-J- 


VENEZUELA 
Barcelona Siegele, Harold L. (J). 
Caracas — Hall, Merle Leroy, Jr. (M) (C/S- 
J-M) xx 





Employmen 


t Notices 





The Journat will post notices of 
men and jobs available. Companies and 
AIME members are invited to use this 
space, for which there is no charge. 
Except as noted below, address replies 
to: Code (appropriate number), Jour 
NAL OF PETROLEUM TECHNOLOGY, 601 
Continental Bldg., Dallas 1. Show re- 
turn address on envelope. These replies 
will be forwarded unopened and no 
fees are involved. 

Replies to the positions coded Y4735. 
Y4658(a). Y4658(b) and Y5642 below 
should be addressed to: Engineering 
Societies Personnel Service. 8 West 40th 

t., New York 18, N. Y. The ESPS, on 
whose behalf these 
lished here, collects a fee from appli 
cants actually placed 


notices are pub- 


PERSONNEI 


@ February 
wishes position related to geology. Has 
MS in geology and is available im- 
mediately for employment. Salary re- 
quired, $325. A veteran of World War 
II and a Lieutenant (jg) in USNR. 
Code 142. 


graduate, 26, married. 


POSITIONS 


% Petroleum engineer with at least five 
years’ experience and capable of test- 
ing and completing wells to 12,000-ft 
depth. Also should be able to handle 
entire course of production. Headquar- 
ters Dallas, Tex. Independent producer 
employer. Salary open. Code 528. 

@ Engineers. (a) Petroleum Enginee: 
Trainee (CNY-1279-W), 21-24, B.S. in 
petroleum or mechanical engineering. 
No experience 
work in oil fields preferred. Will per- 
form general field and office petroleum 
engineering work including bottom-hole 
checking well depths, per- 


required, but summer 


pressures, 
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forming routine analysis work in the 
laboratory, checking and posting pool 
maps, gas oil ratios, etc. Single status. 
Salary, $5,700 a year. (b) Petroleum 
Engineer Trainee (CNY-1280-W) 

24, B.S. in petroleum engineering, pref- 
erably with some experience in petro- 
leum industry. Will be required to learn 
Spanish. Will be trained to obtain a 
general knowledge of the industry, and 
obtain a good idea of a petroleum or 
chemical engineer’s job. Salary if sin- 
gle, $5.700 a year: if married, $6,600 a 
year. Single status for six months. Lo- 
cation, Venezuela. Y4735. 

@ Instrumentation engineer, graduate 
chemical, with a minimum of eight to 
ten years’ experience in chemical ot 
oil refinery work. Must know electric 
and pneumatic instrumentation. Salary, 
$6,000-$7,200 a year. Y4658(a). 

@ Buyer with experience in purchasing 
materials for the chemical industry 
from large petroleum refiners. Salary, 
$4,800-$6,500 a year. Location, Néw 
York, N. Y. Y4658(b). 

@ Engineers. (a) Chief Reservoir En- 
gineer, 25-40, with degree in engineer- 
ing and about five to eight years’ ex- 
perience in petroleum engineering. Lo- 
cation, Peru. (b) Senior Petroleum En- 
gineer, 25-40, with degree in petroleum 
engineering or equivalent. and about 
four to eight years’ experience in res- 
ervoir engineering. Location, Peru. (c) 
Petroleum Engineering Laboratory An- 
alyst. 25-40, with an engineering degree. 
Location, Colombia. Wellsite Mud En- 
gineer, 25-35, with three to five years’ 
experience. Location, Colombia. Y4642. 
@ Petroleum reservoir engineer with 
four years’ experience to work on Tulsa 
engineering staff large independent. 
Supply usual application data including 
salary. Replies will be confiidential. 
P. O. Box 2039, Tulsa. x 8's 
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Rock Bit Design 


Continued from Page 21, Section | 
footage. Several tables were shown to 
indicate that rate of penetration is not 
always desirable, especially in very 
deep holes where round trip time is 
considerable. A cross section of the jet- 
type rock bit was shown. Because of 
the extreme mud velocity in bits of this 
type, the nozzles must be well stream- 
lined and made of special abrasive re- 
sistant materials. 
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BIBLIOGRAPHIES 
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1948. $4.00. 


Precision Investment 
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1949. $2.50. 
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Materials Handling. 114 
references. 1949. $2.00. 


Machinery Foundations; 
Design, Construction, Vi- 
bration Elimination. 120 


references. 1950. $2.00. 
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$2.00. 
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1950. 
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Applications of Solar 
Heating. 149 references. 
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Management of Construc- 
tion Jobs. 53 references. 
1950. $2.00. 
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Memorial | 


4n appreciation by D. R. McKeithan 
and W. A. Schaeffer, Jr 


Alvin D. Whitman (Member 1924) 
died of a heart attack on January 3 
1951. in Denison, Tex. He was born in 
1892, in Pulaski, Va.. but his earlier 
years were spent in lola, Kans., and 
Bartlesville, Okla. When 16 and for the 
next four years he 
worked with his 
uncle, the late T. 

J. Jones, well 

known mining and 

metallurgical en- 

gineer, and was 

engaged in dia- 

mond drilling and 

other activity in 

various western 

states and in Alaska and South 
ica. 


\mer- 


Returning to Oklahoma in 1912, he 
became a tool dresser and cable tool 
driller, and in 1915 went to Burma as 
drilling superintendent for Nathsingh 
Oil Co. He served in France for two 
years as first lieutenant in the Corps 
of Engineers during World I, and be 
fore returning home married Miss Car 
olyn Kingston, of Manchester, England 
During the period 1920-29 he was drill 
ing superintendent for Twin States Oil 
Co. and was with Shell as superintend 
ent in West Texas 1929-31. With discoy 
ery of the East Texas Field in 1931, he 
entered the drilling contracting busi- 
ness, in which he remained as contrac- 
tor and independent operator until his 
death. His headquarters and home had 
been in Houston, Tex., since 1934. Dur 
ing World War II, he 
the Corps of Engineers, this time as a 
major. 


again served in 


Well read. a keen student. and pos 
sessed of broad engineering and me 
chanical ingenuity, he could claim a 
share in the development of many im- 
provements in drilling and producing 
equipment and methods. x * * 





Meetings of Interest 

Donald L. Katz will present a se- 
ries of lectures on Natural Gas En- 
gineering at the University of Kan- 
sas, March 26-27. Advance registra- 
tions should be sent to E. A. Me- 
Farland, Mgr., Lawrence Center, 
University of Kansas, 
Kansas. The enrollment fee is $7.50. 
Accommodations are available in 
fraternity houses at $2.50 up and at 
local hotels from $3.00 up 


I awrence, 
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CHONS.... 
to the editor 





To the editor: 

Now that criticism is being asked 
for with respect to our publication, I 
will give you the opinion of an ordi- 
nary AIME member. 

Part of the available material for 
publication seems to consist of papers 
previously given in large meetings. 
Naturally these will be old stuff to 
some but this, per se, should not be 
the sole objection. The main objec- 
tion to some papers applies . . . to 
their being given originally. 

For instance, there have been those 
that were by rights suitable for organ- 
izations of high class physicists or 
members of the Faraday Society. Not 
one per cent of our membership could 
follow or understand them. The sum- 
mary ... of these papers may be of 
benefit and can be understood. I am 
not belittling the fine work of these 
extraordinary men but do believe that 
papers could be given without 
vast amounts of Greek letter formulae 
that are impressive but sleep-produc- 
ing when not understood. Anything 
worthwhile can be expressed so that 
the organization as a whole can assim- 
ilate and make use of it... 

Neither am I saying that a talk 
should be limited to high school 
mathematics. There is a happy me- 
dium, however, and very few members 
have continued their mathematical 
studies since leaving school; hence. 
while they appreciate the work and 
its necessity, the actual feeling is relief 
when the paper is over. 

We have too much data on oil and 
gas developments in the various areas. 
A journal should pertain more to what 
occurs in (its) own organization, with 
selected papers that may be by non- 
members, which will fit the picture. 
Irving Fatt has some excellent sug- 
gestion in the final paragraph of 
his letter (PETROLEUM TECHNOLOGY. 
Dec., 1950). But papers on the sub- 
jects he suggests should be presented 
in a manner that will not be overpow- 
ering and only read by a minority 
group. 

I am not trying to reduce our or- 
ganization to a lower plane in making 
these suggestions. We can have erudite 
papers but the manner and form 
should be changed. A lecture does not 
have to be in Latin to have class . . . 

Cc, E. Clason 
Halliburton Oil Well 
Cementing Co. 
Duncan, Okla. 
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Technical Note 67 
Continued from Page 17, Section | 

the conversion factors do not depend significantly on the con 
centration, the 0.5 factor for the SO; is actually the mean 
determined for three solutions of Na.SO,, using SO™ concen 
trations of 676 ppm, 2,030 ppm, and 3.030 ppm. The resistiy 
ity conversion factors found for these quantities of SO, were 
0.553, 0.532, and 0.412 respectively. It is therefore recom 
mended that a factor of 0.5 for SO; be used for the majority 
of oil field brine samples, in which the SO- concentration is 
relatively small compared to the C/- concentration. However: 
in some regions, such as the Rocky Mountains, the SO, con- 
centration is large enough to be relatively important, and in 
these areas the factor corresponding to the existing SO; con- 
centration should be used. 

In order to check the accuracy of this method of calculating 
formation water resistivities, the resistivities of formation 
waters from 26 wells in the Gulf Coast region were measured 
in a four pole resistivity cell, using 60 cycle current and a 
vacuum tube voltmeter of 1 megohm input impedence. The 
water resistivities ranged from 0.0972M at 68°F to 3.522M at 
68°F. The chemical analyses of the waters were also available. 
and the resistivities of these samples were calculated from 
their analyses according to the methods outlined by Jones 
and by Tixier, as well as by the method described above. 
using resistivity vs. concentration data taken from the Inter 
national Critical Tables. 

For the method presented here, the average of the absolute 
values of the per cent deviation between the measured and 
calculated resistivities was 3.3 per cent, the calculated resis- 
tivity being on the average 1.1 per cent too great. For the 
Park Jones method, which uses a series of conversion factors 
to change the primary and secondary salinity and alkalinity 
to equivalent quantities of NaCl, the average of the absolute 
values of the per cent deviation between the measured and 
calculated resistivities was 4.6 per cent, with the calculated 
resistivity being on the average 3.2 per cent too great. For 
the Tixier method, which adds the reciprocals of the resistivi 
ties of the salts formed in the Palmer recombination scheme 
to obtain the reciprocal of the water resistivity, the averag¢ 
of the absolute value of the per cent deviation was 7.8 per 
cent with the calculated resistivity being on the average 6.3 
per cent too low. Tixier also describes a simplified method of 
calculation which is identical with that of Park Jones, except 
that a factor of 0.32 rather than 0.37 is used to convert from 
secondary salinity to equivalent NaCl. Incidentally, in both 
systems, the 0.4 conversion factor used to go from primary 
salinity as Na,SO, to equivalent NaCl is almost certainly in 
error. In laboratory prepared solutions of high SO, and low 
Cl- content, use of this factor led to resistivity errors as large 
as 60 per cent. As noted before, high SO; concentrations are 
encountered mainly in the Rocky Mountain regions. 

The principal advantage of the method presented here is 
its speed and convenience. It is felt that for most purposes 
and in particular where the SO; concentration is low, any of 
the three methods of obtaining water resistivity will give re 
liable results. 
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For Those Who 
Wissed Out 


Petroleum Branch Members who failed 
to attend the AIME Annual Meeting in 
St. Louis, Mo., Feb. 18-22, may now 
receive the 10 papers presented there 


These Papers Are Now 
Available in Preprint 
Form For Only 
$2.00 Per Set 


And In Addition — 


Production Review Reports 
For All Domestic Areas 
Are Included 
At No Extra Cost! 


You Get This — 
3 Papers on Waterflooding 


3 Drilling, Logging and 
Acidizing Papers 


¢ Reservoir and Production 
Engineering Papers 
And — 


Production Review Reports of 
All Domestic Areas 


ALL FOR ONLY $2.00 


Petroleum Branch, AIME 
601 Continental Bldg. 
Dallas 1, Texas 


Please send me sets of papers presented at 
AIME Annual Meeting in St. Louis,. Enclosed is my 
check (or money order) in the amount of $ 


Name 
Address 


City, State 
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Darwin H. Clark Co., Los Angeles, Calif. 


Lufkin Foundry & Machine Co. . . . . Sec. 1-24 
Cudlipp Advertising Agency, Lufkin, ten 


Magnet Cove Barium Corp. ; Third Cover 
Rives, Dyke and Co., Houston, Tones 


B and W Incorporated . 
Gene fees, Gere... 3 ee 
Harold Warner Co., Buffalo, N. Y. 


Core Laboratories, Inc. , : 
Duvall Williams, Dallas, eo 


Oil Base, Inc. ar oe! ge 
Dozier-Graham- Susmen, bes Aageie, Colit. 


Otis Pressure Control, Inc. 
Dan Goodrich, Dallas, Tex 


Professional Services «... . « «» « « « « See 2 
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Dowell Incorporated 
MacManus, John & Adame, hee Detroit, Mich. 


Eastman Oil Well Survey Co. . . ‘ 
Ted Workman Advertising, Dollas, en. 


G. M. Giannini & Co., Inc. . 
Western Advertising Ageney, los ‘Angeles, Calif. 


Grove Regulator Co. . . 
Rives, Dyke and Co., Menten, Sex, 


Halliburton Oil Well Cementing Co. . 
Lannan & Sanders, Dallas, Tex. 


Harry W. Dietert Co. 


Johnston Oil Field Service Corp. ‘ 
Wallace Davis & Co., Houston, Tex. 


The Second National Bank, Houston, Tex. 
Gano, Bachrodt, Edwards, Inc 


Schlumberger Well Surveying Corp. Second Cover 


Rives, Dyke & Co., Houston, Tex. 


The Geolograph Co., Inc. Sec. 2—3 


Lowe Runkle Co., Oklchoma City, Okla. 
United Geophysical Co., Inc. 
Dozier-Graham-Eastman, Los Angeles, Colif. 
Welex Jet Services, Inc. . . 
Albert Evans & Assoc., Ft. Worth, Tex 
* In previous Issues 





PETROLEUM BRANCH MEETING CALENDAR 
MARCH Through AUGUST 


MAY 


South Plains Sub-Section 3 
8 Delta Section 5 


Illinois Basin Chapter 3 Delta Section 


Oklahoma City Section Nlinois Basin Chapter 


9 Southwest Texas Section 9 ee 
14 Mid-Continent Section 


Pacific Petroleum Chapter 18 
16 Kansas Section 


17 South Plains Sub-Section 
17 Oklahoma City Section 19 
28 Pacific Petroleum Chapter 23 


Kansas Section 
Kansas Section 

19 South Plains Sub-Section 
Oklahoma City Section 


Pacific Petroleum Chapter 


Illinois Basin Chapter AUGUST 
Mid-Continent Section 5 Delta Section 2 


6 Southwest Texas Section 7 


Illinois Basin Chapter 


Delta Section Delta Section 


7 Illinois Basin Chapter 
1 Mid-Continent Section 


‘ : 15 
ansas Section 20 Kansas Section ’ 


South Plains Sub-Section 21 
Oklahoma City Section 21 


Southwest Texas Section 13 Mid-Continent Section 


Kansas Section 


South Plains Sub-Section 16 South Plains Sub-Section 


Oklahoma City Section 16 
Pacific Petroleum Chapter 27 


Oklahoma City Section 


Pacific Petroleum Chapter 25 Pacific Petroleum Chapter 
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CEMENT CONTAISIWATIOM 


Several times during the drilling of every well cement contamination 


of mud becomes a problem. These Magcobar Products have been proved 
*ractive for pre-treating the mud prior to drilling the cement and restoring the desirable 
drilling qualities to cement-cut mud. When you encounter cement contamination, 
use these Magcobar materials . . . separately or in combinations 
. to save rig time and money and eliminate the costly practice of 
d scarding cement-cut mud. Your Magcobar field engineer will gladly give assistance 
at the rig and there are nearly 300 strategically located Magcobar dealers 
with stocks of these products, ready to serve you day and night 


MAGNET COVE BARIUM CORP., MALVERN ARK., HOUSTON, TEX. 
ONE OF THE DRESSER INDUSTRIES 


i COI ay, : 


Complete 
DRILLING MUD SERVICE 
DEALER! 





YOU GET ALL THREE WITH 


LANE-WELLS 
PERFORATING 


“And with Lane-Wells, they all come 
first—because they’re all equally important. 





Penetration you've got to have—sure—and 
Lane-Wells delivers it. But you've got to 
have it in the right place, and you've got to 


get it there without risking your well or 


eewemme your men. And Lane-Wells does just that. I 


ye don't believe you can equal them in 
/ 


/ 
th 


measurement accuracy, equipment or oil- 
country knowledge. They've done a lot oi 
jobs for me, and I know they're good! 


h Yor 
LANE@ WELLS 


LOS ANGELES + HOUSTON + OKLAHOMA CITY 
General Offices, Export Office and Plant + 5610 SO. SOTO STREET, LOS ANGELES 58, CALIFORNIA 
LANE-WELLS CANADIAN CO. IN CANADA © PETRO-TECH SERVICE CO. IN VENEZUELA 





